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To  quantify  the  performance  abilities  of  existing  or  proposed  narigation  systems,  the  U.S.  Air  Force  has,  for 
the  last  several  y»rs,  compared  the  performance  of  the  system  under  test  to  the  p:erformance  of  a  baseline 
navigation  system  known  as  the  Completely  Integrated  Reference  Instrumentation  System  (CIRIS).  CIRIS 
obtains  a  highly  accurate  navigation  solution  by  combining  information  from  three  major  subsystems:  inertial 
navigation  system  (INS)  information,  barometric  altitude  information,  and  range  and  range-rate  information 
from  ground  transponders  which  have  b-en  precisely  surveyed.  Although  the  navigation  solution  produced  by 
CIRIS  is  highly  accurate,  it  will  soon  be  inadequate  as  the  standard  against  which  future  navigation  systems  can 
be  tested.  This  research  proposes  an  alternative  to  CIRIS  -  a  hybrid  Enhanced  Navigation  Reference  System 
(ENRS)  which  is  designed  to  take  advantage  of  a  newer  INS  (the  LN-B3),  certain  features  of  the  current  CIRIS, 
and  certain  features  of  dilTerential  corrections  to  Global  Positioning  System  (DGPS)  measurements.  Analysis  is 
conduct^  using  a  Kalman  Alter  development  package  known  as  the  Multimode  Simulation  for  Optimal  Filter 
Evaluation  (MSOFE).  Both  a  large  order  imth  model  for  the  ENRS  (in  which  a  full  24  satellite  constellation  is 
modeled)  and  full-  and  reduced-order  Kalman  filters  are  developed.  Results  suggest  (hat  the  proposed  ENRS 
.  (with  DGPS  aiding)  provides  a  navigntion  solution  one  ordes  of  magnitude  better  than  CIRIS.  ^ 
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Preface 


This  thesis  is  the  fourth  in  a  series  devoted  to  increasing  the  accuracy  of  the 
Completely  Integrated  Reference  Instrumentation  System  (CIRIS)  used  by  the  Cen¬ 
tral  Inertial  Guidance  Test  Facility  (CIGTF),  Holloman  AFB,  NM  to  validate  the 
accuracy  of  inertial  navigation  systems  (INSs).  However,  the  material  contained 
within  these  pages  is  not  limited  only  to  this  application.  The  incorporation  of 
Global  Positioning  System  (GPS)  and  Differential  GPS  (DGPS)  mecisurements  ex¬ 
tends  to  any  INS  in  production  or  development  through  the  use  of  Kalman  filtering. 
As  GPS  usage  expands  throughout  the  military  and  civilian  communities,  I  hope 
this  thesis  provides  a  small  contribution  in  this  area. 

As  the  fourth  in  a  series  of  theses,  I  would  be  delinquent  if  I  did  not  mention 
the  people  responsible  for  the  excellent  work  which  proceeded  this  thesis.  Captains 
Solomon,  Snodgrciss,  and  Staeey  all  deserve  more  thanks  than  I  can  give  them  in  this 
format,  so  I  will  just  say  “well  done”.  Their  previous  theses  provided  the  backbone 
from  which  my  DGPS  work  grew,  and  I  would  not  have  completed  one-fourth  of  the 
work  I  accomplished  were  it  not  for  their  legjicy. 

I  would  also  like  to  thank  Mr  Darwin  Abbey  and  Mr  Scott  Dance  of  Intermetrics 
for  taking  time  out  of  their  busy  schedule  to  explain  DGPS  to  an  apprentice.  The 
DGPS  filters  within  this  thesis  can  largely  be  attributed  to  their  recommendations, 
and  I  hope  performance  is  up  to  their  expectations.  I  can  honestly  say  I  am  looking 
forward  to  working  with  these  gentlemen  when  I  get  to  CIGTF  to  further  refine  the 
ba.sic  filters  contained  in  this  thesis. 

Next,  a  few  kind  words  to  the  gentlemen  on  my  thesis  committee  are  in  order. 
Captain  Randall  Paschall,  my  advisor,  kept  me  “on  course"  throughout  this  work 
and  his  guidance  allowed  me  to  complete  as  much  work  as  I  did.  Thattks,  Randy, 
for  showing  me  the  light  at  the  end  of  the  tunnel;  Dr  Peter  Maybeck,  a  visionary  in 


the  field  of  stochastic  estimation,  gave  me  the  clear  insights  on  filter  order  reduction 
and  tuning  used  throughout  this  thesis.  Without  his  advice,  this  thesis  would  not 
be  worth  the  paper  it  is  printed  on.  Lt  Col  David  Meer  and  Maj  Robert  Riggins 
provided  their  time  and  enthusiasm  over  my  work  when  it  wa.s  needed,  and  I  thank 
them  for  their  efforts. 

A  special  word  of  thanks  must  go  to  Mr  Dan  Zambon,  a  man  greatly  responsible 
for  my  success  in  completing  this  thesis.  Dan  gave  me  the  opportunity  to  utilize  six 
different  computer  systems  and  over  300  mega-bytes  of  storage  space  for  all  my 
simulations  and  results.  Were  it  hot  for  his  dedicated  efforts  keeping  the  computer 
systems  working  at  high  efficiency,  I  would  never  have  completed  the  first  Monte 
Carlo  simulation  or  tuned  the  first  filter  let  alone  the  three,  thousand  Monte  Carlo 
runs  this  thesis  encompasses. 

Also,  my  parents  must  also  be  recognized  for  instilling  in  me  the  need  for  higher 
education.  Their  support  throughout  the  years  made  it  possible  for  me  to  complete 
this  thesis  and  earn  yet  another  degree.  I  am  only  now  beginning  to  realize  how 
much  their  guidance  throughout  my  life  has  contributed  to  my  current  success. 

Finally,  a  word  to  my  wonderful  wife  Annette.  I  will  never  be  able  to  repay 
you  for  the  many  months  you  were  neglected  so  that  I  could  complete  this  thesis. 
Your  endless  encouragement  is  the  driving  force  making  all  my  successes  possible.  I 
thank  God  for  giving  me  someone  as  special  as  you  to  spend  the  rest  of  my  life  with. 

William  Joseph  Negast 
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AFIT/GE/ENG/91D-41 

Abstract 

To  quantify  the  performance  abilities  of  existing  or  proposed  navigation  sys¬ 
tems,  the  U.S.  Air  Force  has,  for  the  last  several  years,  compared  the  performance 
of  the  system  under  test  to  the  performance  of  a  baseline  navigation  system  known 
as  the  Completely  Integrated  Reference  Instrumentation  System  (CIRIS).  CIRIS 
obtains  a  highly  accurate  x^avigation  solution  by  combining  information  from  three 
major  subsystems:  inertial  navigation  system  (INS)  information,  barometric  altitude 
inforination,  and  range  and  range-rate  information  from  ground  transponders  which 
have  been  precisely  surveyed.  Although  the  navigation  solution  produced  by  CIRIS 
is  highly  accurate,  it  will  soon  be  inadequate  as  the  standard  against  which  future 
navigation  systems  can  be  tested.  This  research  proposes  an  alternative  to  CIRIS 
-  a  b,.  brid  Enhanced  Navigation  Reference  System  (ENRS)  which,  is  designed  to 
take  advantage  of  a  newer  INS  (the  LN-93),  certain  features  of  the  current  CIRIS, 
and  certain  features  of  differential  corrections  to  Global  Positioning  System  (DGPS) 

,  meeisurements.  Analysis  is  conducted  using  a  Kalman  filter  development  package 
known  as  the  Multimode  Simulation  for  Optimal  Filter  Evaluation  (MSOFE).  Both 
a  large  order  truth  model  for  the  ENRS  (in  which  a  full  24  satellite  constellation  is 
modeled)  arid  full-  and  reduced-order  Kalman  filters  are  developed.  Results  suggest 
that  the' proposed  ENRS  (with  DGPS  aiding)  provides  a  navigation  solution  one 
order  of  magnitude  better  than  CIRIS. 
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INCORPORATION  OF  DIFFERENTIAL  GLOBAL 
POSITIONING  SYSTEM  MEASUREMENTS 
USING  AN  EXTENDED  KALMAN  FILTER 
FOR  IMPROVED  REFERENCE 
SYSTEM  PERFORMANCE 


I.  Introduction 

The  Completely  Integrated  Reference  Instrumentation  System  (CIRIS)  is  a 
transponder  aided  Inertial  Navigation  System  (INS)  test  reference  currently  used  by 
the  Air  Force  for  the  development  and  testing  of  new  aircraft  navigation  systems 
(20:1-1).  The  Office  of  Primary  Responsibility  (OPR)  for  the  CIRIS  system  is  the 
Central  Inertial  Guidance  Test  Facility  (CIGTF),  6585th  Test  Group,  Air  Force 
Systems  Command  (AFSC),  Hollomail  AFB,  NM.  The  operation  of  CIRIS  involves 
flying  the  INS  to  be  tested,  referred  to  as  the  test  article,  and  CIRIS  through  an 
aircraft  trajectory  of  interest  across  a  CIRIS  transponder  range.  The  data  from' 
■  each  system,  is  recorded  during  the  flight  and  compared  to  analyze  hov/  well  the  test 
article  performed.  Up  to  this  point  in  time,  CIRIS  was  considered  mdre  accurate 
than  the  test  articles  and  formed  the  baseline  for  determining  the  performance  of 
aircraft  IN.Ss. 

1,1  Background 

CIRIS  determines  the  aircraft’s  latitude  and  longitude  with  a  1«7  accuraicy  of 
13  feet(ft)  horizontal  and  40  ft  vertical;  the  north  and  west  velocity  to  0.1  ft/sec 
(fps)  l<r;  and  the  vertical  velocity  to  0.4  fps  l<r  (20:1-1).  This  accuracy  is  due  to  the 
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Range/ Range- Rate  System  (RRS)  transponder  aiding  (8,  9).  Recently,  there  have 
been  a  few  state-of-the-art  aircraft  INSs  developed  (and  many  more  in  the  design 
stage)  approaching  the  accuracy  of  CIRIS.  Interestingly,  many  of  thsse  new  INSs 
use  the  Global  Positioning  System  (GPS)  to  increase  their  accuracy.  In  order  to  use 
CIRIS  as  a  baseline  against  these  new  INSs,  GIRIS  must  be  enhanced  to  provide  an 
order  of  magnitude  more  accurate  navigation  solution. 


1.2  Problem  Definition 


By  using  Differential  Global  Positioning  System  (DGPS)  measurements  to  aug¬ 
ment  the  navigation  solution  of  CIRIS,  it  should  be  possible  to  increase  the  accuracy 
of  CIRIS  to  produce  an  order  of  magnitude  better  estimate  of  the  navigation  solu¬ 
tion.  Specifically,  this  thesis  concentrates  on  developing  a  post-processing  Extended 
Kalman  Filter  (EKF)  of  70  states  or  less  to  augment  the  CIRIS  navigation  solution 
with  DGPS  pseudorange  and  delta-range  measurements.  An  EKF  is  developed  in¬ 
stead  of  a  smoother  due  to  the  limited  computer  storage  capacity  available.  The  EKF 


is  limited  to  70  states  or  less  to  ensure  24  hour  turn-around  time  for  post-processing 
real  measurements  during  INS  testing  at  CIGTF.  Pseudorange  and  delta-range  mea¬ 


surements  are  explained  in  detail  in  the  DGPS  sections  of  Chapter  IV. 

1  2.1  Summary  of  Effort.  The  Navigation  Reference  System  (NRS)  filter 
Stacey  designed  (22)  is  tuned  by  Monte  Carlo  analysis  and  its  performance  relative 
to  CIR[S  is  analyzed.  A  41  error-state  reduced-order  INS  filter  model  is  developed 
based  on  the  recommendations  of  Lewantowicz  and  Keen  (10).  97-  and  69-state  NRS 
filter  models  are  evaluated  and  compared  to  Stacey’s  NRS  filter  and  CIRIS.  A  DGPS 
error  model  is  developed  and  augmented  to  the  INS  ano  RRS  error  models  to  form 
the  89-iitate  Enhanced  Navigation  Reference  System  (ENRS)  truth  and  filter  mod¬ 
els.  Finally,  a  reduced-order  ENRS  filter  of  48  error-states  incorporating  delta-range 
and  range-rate  measurements  is  designed,  developed,  and  compared  to  the  NRS  and 
CIRIS. 
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1.3  Summary  of  Previous  Research 

Previous  research  at  AFIT  by  Captains  Joe  Solomon  (GE-89M),  Britt  Snod¬ 
grass  (GE-90M),  and  Richard  Stacey  (GE-91M)  focused  on  improving  CIRIS  per¬ 
formance  (20,  21,  19,  22).  Solomon  modified  the  CIRIS  error  model,  produced  an 
extended  Kalman  filter  designed  to  enhance  CIRIS  performance  (20),  and  performed 
,  preliminary  work  on  the  GPS  error-state  model  (21).  Snodgrass  evaluated  and  modi¬ 
fied  Solomon’s  filter  design  (19).  Stacey  augmented  the  LN-93  baro-altimeter  model, 
impleihented  a  reduced-order  Ih(S  error  model,  designed  high  order  CIRIS/GPS  truth 
and  filter  models,  and  developed  a  Kalman  filter  to  integrate  GPS  pseudorange  mea¬ 
surements  into  CIRIS  (22). 

Capt  Solomon  conceptrated  his  effort  on  designing  an  improved  Kalman  filter 
for  CIRIS  (20).  He  used  an  85-state  truth  model  for  the  LN-39  INS  and  42  states  to 
model  lb  RRS  transponders  for  a  total  of  127  states  in  his  full-state  truth  model.  He 
designed  a  127-state  Kalmcin  filter  using  this  design  and  verified  its  performance  using 
actual  CIRIS  flight  data.  He  then  developed  a  reduced-order  70-state  Kalman  filter 
which  performed  as  well  eis  the  full-order  filter  but  required  much  less  computational 
time  to  implement.  In  a  special  study,  Solomon  designed  a  28-state  error  model  for 
the  GPS  system  assuming  stationary  satellites  (21).  Solomon  never  developed  an 
integrated  CIRIS/GPS  filter. 

Capt  Snodgreiss  continued  the  development  of  Solomon’s  70-state  reduced- 
order  filter  (19).  He  improved  the  filter  software  by  modifying  the  structure  of 
the  Source  code,  making  it  easier  to  change  filter  parameters,  thus  in'"'*asing  its 
efficiency  and  capabilities.  He  then  tuned  various  filter  parameters  which  improved 
filter  performance.  Finally,  he  designed  a  fixed-interval  smoothing  algorithm  to  use 
as  a  post-processor  to  generate  the  optimal  reference  trajectory.  Snodgrass  was  never 
able  to  evaluate  his  smoother  due  to  limitations  in  time,  and  computer  memory. 

Capt  Stacey  returned  to  a  93-state  truth  model  for  the  INS  in  his  thesis  when 
he  incorporated  the  error-model  for  a  Litton  LN-93  INS  (11)  instead  of  the  LN-39 


INS  used  by  Solomon  and  Snodgrass  (20,  19,  22).  He  increased  this  to  a  96-stai,e 
model  through  improvement  of  the  baro-altimeter  model  by  using  a  four-state  model 
instead  of  a  one-state  model.  He  then  designed  a  72-state  reduced-order  INS  model 
and  filter  to  reduce  computational  times.  He  augmented  this  reduced-order  model 
with  26  RRS  transponder  error-states  and  30  GPS  error-states  for  a  128-state  truth 
and  filter  model.  He  completed  one  simulation  run  with  this  model,  referred  to  as 
the  Navigation  Reference  System.  It  is  important  to  note  Stacey  never  implemented 
GPS  delta- range  or  RRS  range-rate  measurement  equations  in  his  research. 

1.4  Research  Objectives 

This  research  finishes  developing  the  NRS  filter  started  by  Stacey.  The  full- 
order  Kalman  filter  is  analyzed  and  the  performance  increa.se  over  CIRIS  is  charac¬ 
terized.  Then  both  the  truth  and  filter  models  are  reduced  by  decreasing  the  72-state 
INS  model  to  41  states,  thus  decreasing  the  full-order  truth  and  filter  models  to  97 
states  from  Stacey’s  128  states.  This  decrease  is  necessary  for  a  timely  completion 
of  this  thesis  and  matches  the  results  of  the  higher  order  truth  and  filter  models 
developed  by  Stacey  (22). 

An  Enhanced  Navigation  Reference  System  (ENRS)  is  then  developed.  The 
ENRS  augments  standard  CIRIS  with  OGPS  iiseudorange  and  delta-range  measure¬ 
ments.  The  ENRS  is  evaluated  against  the  NRS  as  well  as  CIRIS  to  show  the 
performance  increase  available  from  augmenting  CIRIS  with  DGPS  measurements. 
By  increasing  accuracy  using  the  ENRS,  the  position  and  velocity  estimates  are  more 
accurate  than  the  current  CIRIS’s  navigation  solution. 

The  following  i^bjectives  are  based  on  the  statement  of  the  problem  as  pre¬ 
sented  in  the  previous  sections. 

1.  Perform  10-run  Monte  Carlo  analysis  on  the  128-state  NRS  filter  designed  by 
Stacey  (22).  Compare  these  results  to  CIRIS  performance  and  show  the  increment  in 
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perform^ce  possible  with  the  NRS.  This  analysis  completes  the  work  Stacy  started 
on  this  full-order  filter  and  completes  all  of  the  stated  objectives  in  his  thesis. 

2.  Decrease  both  the  truth  and  filter  models  of  the  INS  down  to  4  T  states  from 
the  72  states  Stacey  implemented.  Show  this  reduced  state  truth  model  calculates 
position,  velocity,  and  attitude  error-states  with  no  discernible  difference  from  the 
93  error-state,  INS  model  (one  baro-state)  but  reduces  computational  times  greatly. 

3.  Develop  a  new  full-order  NRS  EKF  incorporating  the  above  reduced-state  INS 
truth  model.  This  97-error-state  filter’s  performance  is  compared  to  CIRIS  as  well  as 
Stacey’s  128-error-state  NRS  filter.  This  proves  that  the  navigation  solution  accuracy 
of  CIRIS  can  be  improved  with  the  incorporation  of  GPS  pseudorange  measurements. 
A  69-error-state  reduced-order  NRS  filter  is  implemented  and  its  performance  com¬ 
pared  to  CIRIS,  Stacey’s  NRS  filter,  and  the  97-state  full-order  NRS  filter. 

4.  Develop  a  DGPS  error-state  truth  model.  This  model  is  similar  to  the  GPS 
truth  model  but  produces  less  variance  in  the  position  estimates  due  to  removal  of 
most  of  the  atmospheric  and  ephemeris  errors  and  elimination  of  the  space  vehicle 
(SV)  clock  errors  (1,  5,  18). 

5.  Design  and  implement  an  ENRS  EKF  using  the  above  DGPS  truth  model  and 
the  CIRIS  truth  model  which  has  been  previously  developed  (20,  19,  22).  This  89- 
state  full-order  filter  differs  from  the  97-state  NRS  filter  only  by  a  change  in  the  GPS 
states  due  to  implerpentalion  of  a  DGPS  nmdel.  The  performance  of  the  ENRS  filter 
is  compared  to  CIRIS  and  to  the  NRS  filters  to  show  the  position  accuracy  improve¬ 
ment  possible  with  the  incorporation  of  DGPS  ps«*udorange  .measurejiients.  Then  a 
48-error-state  reduced-order  ENRS  filter  i.s  «levelo|>e<l  to  r«Hluce  computational  tinM*s 
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while  maintaining  accuracy.  Its  performance  is  compared  to  CIRIS,  Stacey’s 
NRS  filter,  and  the  97-  and  69-error- state  NRS  filters  described  above. 

6.  Design  and  implement  DGPS  delta-range  and  RRS  range-rate  measurements  into 
the  89-state  ENRS  truth  model  and  the  48-state  EN'RS  filter  to  improve  the  velocity 
error  estimation.  Show  the  velocity  estimates  with  this  filter  are  more  accurate  than 
CIRIS  and  the  48-state  filter  above  having  no  velocity  measurements. 

7.  Compare  the  accuracy  of  the  48-state  ENRS  filter  to  all  previous  filters.  This 
comparison  shows  the  reduced-order  ENRS  filter  with  delta-range  and  range-rate 
measurements  outperforms  all  the  NRS  filters  and  has  an  order  of  magnitude  in¬ 
crease  in  position-error  accuracy  and  an  increa.se  in  velocity -error  accuracy  over  the 
CIRIS  navigation  solution. 

1.5  Research  Approach 

Before  any  filter  analysis  can  be  accomplished,  the  Multimorle  Simulation  for 
Optima]  Filter  PN-aluation  (MSOFE)  software  package  must  be  thoroughly  under¬ 
stood  (3).  MSOFE,  as  u.sed  in  this  thesis,  propagates  and  updates  the  EKF'  and 
also  propagates  the  truth  model  of  the  .system.  As  in  Stacey’s  thesis,  a  two-hour 
fighter  flight  profile  is  used  to  generate  the  triith  model  for  each  analysis  (22).  Once 
MSOFE  i*  mastered,  filter  evaluation  and  tuning  is  accomplished.  Chapter  II  con¬ 
tains  an  overview  of  the  EKF  and  also  describes  filter  order  retluction  and  tuning. 
The  remainder  of  this  section  gives  a  brief  overview  of  the  approach  to  completing 
each  of  the  research  objectives  outlined  in  the  previous  section. 

1,  Stacey’s  truth  and  filter  model  is  evaluated  using  lO-rtin  .Monle-(’arlo  analy¬ 
sis  and  the  NRS  filter  performance  is  comparer!  to  ('IRIS  to  show  the  inrrea.s«*  in 
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navigation  solution  accuracy  p)ossible  by  incorporating  GPS  pseudorange  measure¬ 
ments  into  CIRIS.  Stacey’s  NRS  truth  and  filter  models  are  described  in  Chapter  III 
and  the  results  are  analyzed  in  Chapter  V. 

2.  Stacey’s  truth  and  filter  models  are  based  on  1.28  error-states.  It  presently  takes 
over  18  hours  for  a  single  MSOFE  one-run  analysis  to  be  completed  for  a  2-hour 
simulated  fighter  flight  trajectory.  A  10-run  Monte  Carlo  analysis  is  considered  the 
minimum  which  should  be  completed  before  any  reasonable  estimate  of  filler  perfor¬ 
mance  can  be  evaluated.  Decreasing  the  INS  truth  model  from  72  states  to  41  states 
adlows  a  10-run  analysis  to  be  completed  in  only  4  days.  However,  for  the  41-state 
INS  model  to  be  considered  adequate,  it  must  provide  a  irue  navigation  solution  as 
accurate  as  the  full  LN-93  OS-state  INS  truth  model.  Using  MSOFE  and  the  rec¬ 
ommendations  of  Lewantowicz  and  Keen,  the  41 -state  INS  truth  model  is  evaluated 
against  the  93-state  INS  truth  model  showing  almost  identical!  performance  (10). 
The  reduction  of  INS  states  is  explained  in  Chapter  III  with  results  in  Chapter  V. 

3.  Once  the  INS  truth  model  is  reduced,  new  full-order  97  error-state  NRS  truth 
auid  filter  models  are  developied.  A  lO^run  Monte  Carlo  analysis  is  performed  on  the 
filter  and  the  results  are  comp2U’ed  to  CIRIS  and  Stacey’s  NRS  filter.  Since  four  days 
is  still  a  long  time  to  wait  for  results  and  the  objective  of  this  thesis  is  to  develop  a 
filter  of  less  than  70  states,  a  69-state  reduced-order  filter  is  implemented  by  elimi-' 
nating  all  but  two  GPS  error-states  (the  receiver  clock  states).  This  filter  completes 
a  10-run  Monte-Carlo  analysis  in  36  hours,  but  with  proper  tuning  (see  Chapter  II) 
produces  a  navigation  solution  just  as  accurate  a.s  Stacey’s  .NRS.  filter  or  the  97-state 
NRS  filter  described  above.  The  work  in  this  section  provides  a  baseline  to  analyze 
how  well  the  DGPS  model  performs.  The  filter  reduction  is  explained  i.T  detail  in 
Chapter  III  with  the  analysis  of  filter  performance  presented  in  Chapter  V. 
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4.  The  DCPS  model  is  very  similar  to  \he  GPS  model  developed  by  Stacey.  The 
differences  lie  in  the  fact  that  the  SV’s  clock  error  is  totally  eliminated  while  the 
atmospheric  and  SV  ephemeris  errors  are  greatly  reduced.  A  complete  discussion  of 
DGPS  as  well  as  a  detailed  analysis  of  the  22  error-states  which  make  up  the  DGPS 
error  model  are  given  in  Chapter  I\'. 

5.  The  89-state  truth  and  filter  models  in  the  full-order  EXRS  are  only  different 
from  the  97-state  NRS  models  discussed  above  by  the  change  in  the  number  of  GPS 
error-states  when  using  DGPS.  A  10-run  Monte  Carlo  analysis  is  performed  and  the 
results  discussed  in  Chapter  IV.  The;n,  in  order  to  meet  the  70-filter-state  require¬ 
ment  and  reduce  the  10-run  time  to  less  than  one  day,  a  48-state  filter  model  is 
implemented  by  again  eliminating  all  but  two  of  the  DGPS  error-states  (the  two 
clock  states)  and  21  more  INS  slates  according  to  the  recommendation  of  Lewan- 
towicz  and  Keen  (10).  This  48-state  filter  is  discus.sed  in  Chapter  IV  with  the  filter 
results  presented  in  Chapter  V. 

6.  DGPS  delta-range  and  RRS  range-rate  measurements  are  both  needed  if  the  ve¬ 
locity  error  estimates  of  the  ENRS  filter  are  to  be  more  accurate  than  CIRIS  which 
uses  range-rate  measurements  alone  for  velocity  aiding.  The  theory  presented  behind 
delta-range  and  range-rate  measurements  is  given  in  Chapter  IV.  A  DGPS-INS  filter 
of  22  error-states  (no  RRS)  and  RRS-INS  filter  of  46  error-stales  (no  DGPS)  are 
used  to  tune  the  velocity  aiding  measurements  separately  for  optimal  performance 
using  a  10-run  Monte  Carlo  analysis.  The  results  of  these  two  filters  are  presenteil 
in  Chapter  V. 

7.  Finally,  delta-range  and  range-rate  measurements  are  incorporated  into  the  48- 
state  ENRS  filter  to  complete  the  ENRS.  A  25-run  Monte  Carlo  analysis  i.s  perforttied 
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for  a  final  compcirison  against  all  previously  discussed  filters.  The  results  of  this  com¬ 
parison  are  presenti^  in  Chapter  V. 

1.6  As.<itiTnptions 

Because  of  the  lack  of  performance  of  Stacey’s  4-state  baro-altimeter  model, 
the  Litton  93-state  error  model  using  one  baro-altimeter  state  is  considered  the  true 
error  model  of  the  INS  subsystem  of  CIRIS.  The  41-state  error  model  implemented 
as  truth  in  this  thesis  is  shown  to  be  a  very  close  approximation  to  the  93-state 
LN-93  error  model.  Snodgrass’s  RRS  error  model  is  used  as  the  truth  model  for  the 
RRS  subsystem  as  it  was  implemented  in  Stacey’s  thesis  (19,  22).  Finally,  Stacey’s 
GPS  error  model  derived  from  articles  by  Cox  (4),  Martin  (12),  and  Milliken  and 
Zoller  (17)  is  considered  adequate  to  use  as  a  starting  point  for  the  DGPS  subsystem 
for  the  ENRS. 

1.7  Overview  of  Thesis 

Chapter  II  provides  a  brief  description  of  the  EKF  theory  used  in  this  thesis. 
Filter  order  reduction  and  tuning  are  dso  presented.  Chapter  III  overviews  the  NRS 
developed  by  Stacey.  An  explanation  of  the  procedure  used  in  reducing  the  INS 
error-state  model  from  72  to  41  states  is  given.  The  97-state  NRS  truth  and  filter 
models  and  the  69-state  reduced-order  filter  model  are.explained.  Chapter  IV  begins 
with  an  explanation  of  the  ENRS.  DGPS  principles  and  error-states  are  presented. 
Detailed  analysis  of  the  DGPS  full-state  truth  and  filter  tnodels  are  shown  along  with 
the  development  of  the  reduced-order  ENRS  EKF  and  the  delta-range  and  range- 
rate  measurement  equations.  Chapter  V  discusses  the  results  of  all  the  Monte  Carlo 
analyses  of  the  above  EKFs.  Chapter  VI  contains  the  conclusions  drawn  from  the 
work  presented  in  this  thesis.  A  summation  of  the  overall  performance  of  the  NRS 
and  ENRS  filters  is  given.  Recommendations  for  further  research  in  this  area  are 
also  addressed. 
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II.  Extended  Kalman  Filtering 


2. 1  Overview 

This  chapter  refamiliarizes  the  informed  reader  with  the  extended  Kalman  filter 
(EKF)  dynamics,  measurement,  propagation  and  update  equations  as  they  pertain 
to  this  thesis,  and  discusses  filter  state  reduction  and  tuning.  The  reader  who  is 
not  familiar  with  Kalman  filtering  is  directed  to  Maybeck’s  textbooks  on  stocheistic 
control  in  which  robust  derivations  of  the  EKF  equations  and  discussions  of  state 
reduction  and  filter  tuning  are  presented  (14,  15,  16). 

2.2  The  Extended  Kalman  Filler 

When  both  the  dynamics  and  measurement  equations  of  a  system  of  interest 
are  constructed  using  a  linear  combination  of  the  states  modeled,  a  linear  Kalman 
filter  can  be  developed.  However,  the  error-state  model  for  this  research  is  not  a 
linear  function  of  the  error-states  in  either  the  dynamics  or  measurement  equations. 
Because  of  this,  a  linearized  Kalman  filter  of  some  form  must  be  used.  The  lin¬ 
earized  Kalman  filter  implemented  here  is  the  extended  Kalman  filter,  which  allows 
relinearization  at  each  measurement  update  time.  The  equations  describing  the  EKF 
used  in  this  research  are  derived  from  a  nonlinear  continuous-time  dynamics  equation 

(14): 

X(0  =  flx(0.f]  +  G(0w(0  (2.1) 

In  this  case,  the  state  dynamics  matrix  f[x(t),t]  is  a  nonlinear  function  of  the  state 
vector  x(t)  and  time  i  .  For  this  thesis,  G(t)  =  I  and  the  white  Gau.ssian  noise 
w(0  is  defined  with  mean: 


=  E{w(f)}  =  0 


(2.2) 


and  noise  strength  Q{t)  of: 


£{w(l)w’'((  +  r)}  =  Q(i)i(r)  (2.3) 

In  addition,  the  discrete-time  measurement  equations  z(t  j)  are  also  a  nonlinear  func¬ 
tion  of  the  state  vector  and  time  (15):  ■ 

z(t,)  =  h{x(t, ),<,]-!- v(<.)  (2.4) 

The  discrete-time  noise  vector  v(t,)  is  also  zero-mean  and  its  covariance  is  described 
by: 

£;{v(<0v^(t,)}  = 

A  system  must  be  linear  in  order  for  the  Kalman  filter  to  yield  the  “optimal”  state 
estimates,  so  the  nonlinear  Elquations  (2.1)  and  (2.4)  are  linearized  by  the  method 
given  in  Mayueck’s  Volume  II  (15).  A  nominal  state  trajectory,  Xn(0  is  generated 
satisfying  x„(to)  =  Xno 

Xn(t)  =  f[Xn(0,<]  (2.6) 

where  f[x„(t),  t]  is  specified  in  Equation  (2.1).  The  nominal  measurements  are 
defined  as: 

- 2n(t.)  =  h[x„(t,),M  (2.7) 

The  “perturbation”  of  the  states  are  obtained  by  subtracting  the  nomin^d  trajectory 
(Equation  (2.6))  from  the  original  nonlinear  dynamics  equation  (Equation  (2.1)):, 

[x(<)-Xn(O]  =  *Ix(t),<l-fIXn(O,0  +  G(Ow(O  (2.8) 

The  equation  above  is  approximated  to  first  order  by  a, Taylor  series  expansion: 

«x„(0  =  F((;Xn(0j+G(0w(0  (2.9) 


R(*f )  for  ti  =  tj 
0  for  ti  ^  tj 
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where  ^x(t)  represents  a  first-order  approximation  of  the  process  [x(t)  —  x„{t)], 
and  F[<;Xn(t)]  is  a  matrix  of  partial  derivatives  of  f{x(<),<]  with  respect  to  x(<) 
evaluated  along  the  nominal  trajectory  x„(t)  (15): 


F[<;x„(0] 


af[x(<),<]| 


(2.10) 


IX=Xa(0 


The  perturbed  discrete-time  measurement  equation  is  derived  in  like  fashion  and  is 
expressed  as  (15): 

^z„(<,)  =  H[<,;x„(<,)] -I- v(li)  (2.11) 


where 


H[<i;x„(t,)]  = 


dh(x(<.),<.]| 


(2.12) 


lx=x„(«,) 


The  nonlinear  continuous-time  dynamics  and  discrete-time  measurement  equations 
have  now  been  successfully  converted  to  linearized  “error-state”  equations.  Estimates 
of  the  whole-valued  quantities  of  interest  are  calculated  from  (15): 


x(<)  =  x„(t)  -h  Sx(t) 


(2.13) 


The  expression  above  for  the  linearized  Kalman  filter  is  useful  provided  that  the 
linearization  assumption  is  not  violated.  However,  if  the  nominal  and  “true”  trajec¬ 
tories  differ  by  too  large  an  amount,  unacceptable  errors  may  result  (15).  It  is  for  this 
reason  that  extended  Kalman  filtering  is  ,  useful  in  many  cases  where  perturbation 
techniques  .done  do  not  suffice.  Extended  Kalman  filtering  allows  for  nlinearizing 
about  newly  declared  nominals  at  each  sample  time,  to  enhance  the  adequacy  of  the 
linearization  process,  and  thus  of  the  resulting  filter  performance  as  well.  (22) 

Discrete-time  measurements  are  incorporated  into  the  EKF  by  the  following 
state  equations  (15): 

K[U)  =  P((-)H’'|((0;x(i-)|  {h[((, );*((, -)|P(l,-)H’'(((i);x(<-)l  +  (2.H) 
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x((+)  =x(li  )  +  K((4{z,-h|x(ij  );(ii))}  (2.15) 

p(in  =  p(*r)  -  K('OH((y;*(ir))P(ir)  (2-i6) 

where 

H[U-,x{t-)]  = 

The  state  estimate  and  cov^iance  are  propagated  from  t,-  to  by  integrating  the 
following  equations  (15): 

x(t/ti)  =  f  [x(t/t.),  t]  (2. 18) 

P{t/U]  =  F[t;x{t/U)]Pii/ti)  +  P{t/U)F'^[mm)]  +  G(0Q(0G^(0  (2T9) 

where 

and  the  initial  conditions  are: 

x(t./<.)  =  x(tt)  (2.21) 

Piti/ti)  =  P{tt)  (2.22) 

The  EKF  equations  given  are  programmed  in  MSOFE  for  optinial  error-state  de¬ 
termination.  It  is  important  to  point  out  MSOFE  uses  the  U-D  Factored  form  of 
the  update  and  propagation  equations  (3).  A  complete  development  of  linear  and 
linearized /extended  Kalman  filters  is  given  by  Maybeck’s  Volumes  1  and  2,  respec¬ 
tively  (14,  15).  An  extensive  explanation  of  Kalman  filtering  as  it  pertains  to  this 
research  is  given  in  Chapter  II  of  Stacey’s  thesis,  and  the  information  presented  in 
this  section  is  paraphrased  from  it  (22). 


F(t;x(t/tO]  = 


af[x(t),  t]  I 


lx=x(«Ai) 


(2.20) 


2.3  Filter  Order  Reduction 


When  an  engineer  is  faced  with  the  task  of  designing  a  Kalman  filter  for  a 
system  with  a  large  state  vector,  it  sometimes  becomes  desirable  or  necessary  to 
decrease  the  number  of  states  modeled  in  the  filter,  often  referred  to  as  “filter  order 
reducf  on.”  How  and  why  the  filter  is  to  be  utilized  are  important  considerations 
in  setting  the  final  number  of  states  modeled.  As  an  example,  a  filter  which  is 
implemented  off-line  on  ,  a  large  mainframe  or  supercomputer  can  include  a  larger 
number  of  filter  states  than  a  filter  implemented  on  an  8-bit  aircraft  navigation 
computer  operating  real-time  to  estimate  INS  position  and  velocity  errors.  Even  in 
the  case  of  this  thesis,  where  a  post-processing  filter  is  developed  to  run  off-line  on 
a  Hewlett-Packard  minicomputer,  time  considerations  limit  the  filter  to  less  than  70 
states  (a  10-run  Monte  Carlo  analysis  of  a  2-hour  flight  trajectory  must  be  completed 
in  24  hours  or  less).  Proper  choice  of  states  from  the  truth  model  is  critical  to  ensure 
accurate  position  and  velocity  error  estimates.  The  next  few  paragraphs  offer  some 
general  recommendations  for  filter  order  reduction  with  a  bias  toward  the  work 
performed  in  this  thesis. 

A  good  starting  point  for  filter  order  reduction  is  determining  which  states  can¬ 
not  be  eliminated.  For  example:  the  position,  velocity,  attitude,  auid  baro-aUimeter 
states  are  the  10  states  of  interest  for  an  INS  and  cannot  be  considered  for  elim¬ 
ination  (assuming  they  are  observable)  (14)’.  Those  readers  familiar  with  aircraft 
Kalman  filters  notes  that  these  states  compose  the  bulk  of  the  states  used  in  aircraft 
navigation  filters.  However,  the  other  83  states  which  compose  the  full-order  LN-93 
INS  error  model  are  states  which  are  not  directly  of  interest,  but  contribute  to  the 
position,  velocity,  and  attitude  errors  (see  the  discussion  of  Stacey’s  INS  model  in 
Chapter  III).  Which  of  the  83  states  can  be  eliminated?  A  smart  engineer  starts  by 
conducting  a  literatilre  search  on  the  type  of  model  he  is  implementing  and  analyzes 
the  v/ork  previously  accomplished  with  good  results.  The  paper  by  Lewaiitowicz  and 
Keen  (10)  is  a  good  example  for  this  thesis.  It  presents  four  different  reduced-order 
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filters  for  the  INS  ranging  from  41  states  down  to  17  states,  all  performing  reason¬ 
ably  well.  This  paper  forms  the  basis  for  both  the  truth  and  filter  state  reductions 
executed  in  this  thesis. 

What  happens  when  there  is  no  literature  on  the  subject?  The  engineer  must 
then  analyze  the  states  mathematically  or  use  software  tools  such  as  MSOFE  (3). 
Normally,  the  time  and  frequency  response  of  the  states  are  determined.  Then,  states 
with  negligible  magnitudes  can  be  eliminated,  or  they  can  be  combined  with  those 
of  larger  magnitudes  and  similar  frequency  responses.  Note  that  filter  “tuning”  is 
often  necessary  when  states  are  combined  or  eliminated,  and  this  subject  is  covered 
in  the  next  section.  The  engineer  must  be  extremely  careful  when  performing  state 
reduction  to  ensure  that  accurate  state  estimates  of  the  quantities  of  interest  are 
mdntained. 

A  filter  analysis  tool  such  as  MSOFE  also  helps  the  engineer  in  other  ways  than 
those  expressed  above  (3).  When  a  full-order  filter’s  state  estimates  are  compared  to 
the  truth  model  of  the  system,  sometimes  even  a  well-tuned  filter  provides  estimates 
of  states  that  do  not  adequately  follow  the  true  states  dynamics.  Normally,  this  lack 
of  filter  performance  in  estimating  states  can  be  attributed  to  a  lack  of  observability 
of  the  states  (14).  Maybeck  states  that  it  is  improper  to  include  unobservable  truth 
model  states  in  the  filter  design  model  (13).  Examples  of  this  phenomena  in  this 
thesis  are  the  GPS  space  vehicle  (.SV)  position  errors  (see  Appendices  D,  E,  and 
G).  The  filter’s  covariance  values  for  these  states  show  an  adequately  tuned  filter. 
However,  the  filter’s  state  estimates  never  follow  the  true  error  states  due  to  a  lack 
of  observability  of  the  states.  MSOFE  allows  the  engineer  to  perform  Monte  Carlo 
or  covariance  analysis  on  the  filter  and  accurately  assess  its  state  estimation  per¬ 
formance  under  a  wide  range  of  conditions.  Maybeck  recommends  using  covariance 
analysis  on  linear  or  linearized  Kalman  filters  because  only  one  simulation  run  is  re¬ 
quired  to  determine  the  true  states  covariance  (14).  However,  this  one  simulation  run 
tadces  a  much  longer  time  to  complete  when  large  dimension  filters  are  implemented 
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and  is  not  as  useful  for  an  EKF  implementation,  both  satisfied  in  this  thesis.  Monte 
Carlo  analysis  provides  quicker  results  if  many-run  analysis  is  not  used,  though  not 
it  is  not  as  accurate  as  the  covariance  analysis.  Often  filter  designers  use  1-  to  5-run 
Monte  Carlo  analysis  at  first  to  get  quickly  “in  the  ballpark”  and  then  use  covariance 
ajialysis  to  determine  true  filter  performance.  Only  Monte  Carlo  analysis  is  used  in 
this  thesis  due  to  time  constraints  and  the  implementation  of  an  EKF. 

A  brief  discussion  of  filter  order  reduction  has  been  given.  The  general  proce¬ 
dure  used  to  eliminate  filter  states  successfully  in  this  thesis  has  been  outlined.  This 
procedure  is  by  no  means  all  inclusive,  but  should  give  the  inexperienced  engineer  a 
foundation  for  performing  filter  order  reduction. 

2.4  Tuning  Kalman  Filters 

When  the  truth  model  of  a  system  of  interest  is  composed  of  linear  dynamics 
and  measurement  equations,  a  linear  Kalman  filter  based  upon  this  system  proves 
to  be  the  “optimal”  estimator  by  almost  any  criteria  (14).  When  a  full-order  filter 
of  the  truth  model  is  implemented,  the  dynamics  and  measurement  noises  used  to 
construct  the  filter  (similar  to  the  noises  used  in  the  nonlinear  filter  derivation  given 
previously)  would  normally  match  the  truth  model  noises.  However,  when  filter 
order  reduction  is  performed  on  this  linear  filter,  the  noise  strengths  may  have  to  be 
adjusted  to  compensate  for  the  missing/combined  states.  Along  similar  lines ,  when  a 
linearized/extended  full-order  Kalman  filter  is  implemented  to  simulate  a  nonlinear 
system,  the  noise  strengths  again  may  have  to  be  adjusted  in  both  the  full-order  and 
reduced-order  filters  to  compensate  for  the  approximations  used  in  imphmenting 
linearized  equations  to  model  a  nonlinear  system  and  for  the  reduction/combination 
of  states.  This  filter  rioise  strength  adjustment,  or  “tuning”,  i.s  overviewed  in  the 
next  few  paragraphs  as  it  applies  to  this  thesis. 

The  dynamics  and  measurement  equations  used  in  this  research  hnve  been 
shown  previously  (Section  2.2).  The  closer  the  full-order  filter  approximates  the 


truth  model,  the  less  the  noise  strength  Q(<)  and  covariance  R^t,)  are  adjusted 
(normally  increased)  to  compensate  for  nonlinearities  in  the  true  dynamics.  For  the 
full-order  filters  implemented  in  this  research,  there  are  three  main  reasons  why  the 
dynamics  noise  strength  Q(i)  on  some  of  the  states  needed  increasing. 

The  first  reason  occurs  when  states  which  have  noise  (or  have  no  noise)  associ¬ 
ated  with  them  in  the  truth  model  need  increased  noise  (or  some  noise)  in  the  filter 
model  to  track  the  dynamics  of  the  true  state  adequately.  Adding  noise  increases 
the  uncertainty  of  the  filter-assumed  state’s  dynamics  to  compensate  for  nonlinear 
behavior  or  to  keep  the  Kalman  filter  gains  in  that  channel  from  going  to  zero.  Usu¬ 
ally  the  states  in  the  filter  needing  this  type  of  noise  adjustment  are  the  states  which 
are  not  part  of  the  measurement  equation.  In  this  thesis,  some  of  the  INS  states  not 
part  of  a  measurement  equation  need  a  small  amount  of  added  noise  to  compensate 
for  nonlinear  dynamics.  If  the  states  are  not  directly  part  of  the  measurement  equa¬ 
tion,  and  the  dynamics  equations  do  not  totally  describe  the  true  states  behavior, 
the  EKF  has  difficulty  estimating  these  states.  The  noise  increase/addition  is  not 
large,  just  enough  to  increase  the  uncertainty  in  the  dynamics  equation  to  account 
for  using  a  linearized/extended  Kalman  filter  (14). 

The  second  reason  that  dynamics  noise  needs  to  b"*  added  into  a  filter  model 
occurs  when  the  filter  variance  of  a  state  goes  negative.  Maybeck  states  that  this 
numerical  difficulty  is  a  normal  occurrence  when  the  range  of  numbers  the  filter 
estimates  is  large  as  in  this  filter  when  some  of  the  INS  states  are  on  the  order  of  l.Oe-6 
and  the  GPS  clock  bias  state  is  often  l.OelO  (14).  The  problem  of  negative  variances 
arises  due  to  the  limited  numerical  precision  of  computers  which  must  multiply  high 
order  matrices  together.  True  state  variances  which  showed  this  tendency  often  have 
a  Value  of  0.0  in  the  truth  model,  and  the  filter  does  a  good  job  of  estimating  the 
state,  but  the  variance  often  goes  negative.  A  sma.Il  amount  of  added  noise  in  the 
filter  keeps  the  variance  of  the  state  positive  and  does  not  degrade  the  filter’s  slate 
estimate. 


Tne  third  and  final  reason  why  dynamics  noise  is  added  involves  filter  order 
reduction.  States  are  eliminated  in  the  INS  model:  states  which  do  not  show  up 
directly  in  the  measurement  equations,  but  which  impact  the  states  which  are  part 
of  the  measurement  equations.  To  compensate  for  the  eliminated  states,  the  noise  in 
,the  states  they  affect  is  increased.  The  noise  increase  is  small,  since  the  eliminated 
states  all  have  small  magnitudes,  but  is  necessary  to  ensure  a  well-tuned  filter.  In 
this  thesis,  the  first  two  INS  states  in  the  reduced-order  filters  relating  to  latitude 
and  longitude  need  a  small  amount  of  noise  added  into  the  filter  to  compensate  for 
the  elimination  of  INS  states. 

There  are  two  reasons  why  the  mea.surement  noise  covariance  R(i,)  is  adjusted 
in  the  filters.  The  first  reason  is  the  same  ^  why  the  dynamics  equation  noise  is  ad¬ 
justed  upward  in  the  full-order  filters,  to  increase  the  uncertainty  in  the  measurement 
equations  due  to  using  linearized  equations  to  model  nonlinear  systems.  Remember, 
the  EKF  only  has  access  to  mecisurements  which  are  a  combination  of  states,  not  to 
the  states  themselves. 

The  second  reaison  for  increasing  the  measurement  noise  R(t,),  and  the  most 
important  and  time  consuming  reason  for  increasing  noise  in  these  filter  models,  is 
due  to  filter  order  reduction.  When  states  are  eliminated  in  the  filter  model  which 
are  part  of  the  measurement  equation,  it  is  necessary  to  increa.se  R(f,).  This  increase 
in  noise  is  often  not  small,  as  some  eliminated  states  have  large  covariance  values. 
A  good  starting  point  for  this  noise  increase  is  to  add  the  variance  of  the  eliminated 
states  from  the  full-order  model  to  the  original  variance  of  the  measurement  noise 
(13).  If  this  is  done  progressively  for  each  eliminated  state,  the  resitu;>,g  R(<i)  is 
often  very  close  to  what  is  needed  for  a  well-tuned  filter. 

An  overview  of  how  and  why  filter  tuning  is  performed  in  this  thesis  is  given. 
Once  again,  these  procedures  are  not  universal  to  all  Kalman  filter  tuning.  They 
are  basic  reasons  why  Kalman  filters  need  to  be  tuned,  and  any  engineer  involved  in 
Kalman  filter  tuning  needs  to  be  familiar  with,  them. 
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2.5  S'ATnmary 

Thi ;  chapter  gives  a  surnmary  of  the  Kalman  filter  techniques  utilized  in  this 
thesis.  The  general  EKF  equations,  filter  order  reduction,  and  filter  tuning  are 
examined,  with  a  bias  toward  how  they  were  accomplished  in  this  thesis.  A  more 
detailed  derivation  of  the  EKF  equations  as  they  pertain  to  this  research  can  be 
found  in  Stacey’s  thesis  (22).  Now  that  the  basic  theory  is  covered,  the  NRS  truth 
and  filter  models  are  described  in  the  next  chapter. 
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III.  The  Navigation  Reference  System  (NRS) 

3.1  Overview 

This  chapter  overviews  the  Litton  LN-93  error  model  and  the  NRS  developed 
by  Stacey  (22)  along  with  the  97-  and  69-stafe  i\RS  models  implemented  in  this 
thesis.  First,  the  original  Litton  LN-93  error-state  model  is  overviewed.  Second, 
Stacey’s  INS,  RRS,  and  GPS  models  are  briefly  described  since  they  form  the  basis 
for  all  the  other  models  and  filters  in  this  thesis.  Third,  the  steps  used  to  reduce 
the  93-state  INS  error-model  to  41  states  while  maintaining  the  accuracy  of  the  full 
93-8tate  model  are  explained.  Fourth,  a  97-slate  NRS  model  incorporating  the  new 
41-state  INS  truth  model  is  described.  Finally,  a  G9-state  reduced-order  NRS  filter 
is  described.  For  a  detailed  explanation  of  Stacey’s  NRS,  refer  to  Chapters  3,  4,  and 
5  of  his  thesis  (22). 

3.2  The  93‘State  LS-93  Error  Model 

This  section  overviews  the  original  93-state  Litton  INS  error  model  (11).  Six 
categories  of  errors  are  associated  with  this  error  model.  The  first  category  represents 
the  “general”  error-states,  states  which  are  combinations  of  several  other  states  in 
the  error  model.  Position,  velocity,  platform  tilt,  and  vertical  channel  errors  are 
in  this  category.  The  second  category  of  errors  consists  of  states  described  by  first- 
order  Markov  processes.  Both  the  gyros  and  accelerometers  have  exfwnentially  time- 
correlated  errors,  and  the  baro- altimeter  state  is  also  included  in  this  group.  The 
third  and  fourth  categories  of  error  states  are  the  gyro  and  accelerometer  bias  errors, 
respectively.  These  states  are  modeled  as  r^mdom  constants  in  the  truth  model.  The 
fifth  category  of  error  states  are  again  first-order  Markov  processes  mo<leling  the 
gyro  and  accelerometer  thermal  transients.  Finally,  the  last  category  of  error  stat<*8 
are  gyro  compliance  errors  mo<lel«l  as  random  const^U1ts.  The  following  equation 


describes  the  overall  93-state  error  model  using  the  six  categories  of  states  (22):' 

Sx  =  [  6x3^  Sx^^  Sx^^  Sx^  (3.1) 

where  is  a  93  x  1  column  vector  and: 

8Xi  represents  the  “general”  error  vector  containing  13  position,  velocity,  attitude, 
and  vertical  channel  errors. 

consists  of  16  gyro,  accelerometer,  and  baro-altimeter  exponentially  time-correlated 
errors,  and  “trend”  states.  These  states  are  modeled  as  first  order  Markov  pro¬ 
cesses  in  both  the  truth  (system)  model  and  in  the  Kalman  filter. 

^Xj  represents  gyro  bias  errors.  These  18  states  are  modeled  as  random  constants 
in  the  truth  model  and  are  modeled  as  random  walks  (with  small  magnitude 
pseudo- noises)  in  the  Kalman  hlter. 

^x^.  is  composed  of  the  accelerometer  bias  error  states.  These  22  states  are  modeled 
in  exactly  the  same  manner  as  the  gyro  bias  states. 

6x^  depicts  accelerometer  and  gyro  initial  thermal  transients.  The  6  thermal  tran¬ 
sient  states  are  first  order  Markov  processes  in  the  system  and  Kalman  filter. 

6x^  models  the  gyro  compliamce  errors.  These  18  error  states  are  modeled  as  biases 
in  the  system  model  and  as  random  walks  in  the  Kalman  filter. 

The  truth  model  system  state  space  differential  equation  is  of  the  form: 

F||  Fi2  F|;i  Fi4  F|5  Fie  w, 

0  Fn  0  0  0  0  #Xa  w, 

0  0  0  0  0  0  6x3  0 

0  0  0  0  0  0  8x^  0 

0  0  0  0  Fss  0  6X5  0 

0  0  0  0  0  0  J  [  «x«  J  [0 
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Appendix  A,  Tables  A.l  through  A. 4,  contains  a  full  listing  of  the  93  INS  error  states 
as  well  as  the  RRS  (Table  A. 5)  and  GPS  (Table  A. 6)  states  which  are  explained 
in  later  sections  of  this  chapter.  Appendix  B,  Tables  B.l  through  B.IO  contain  the 
non-zero  entries  of  the  LN-93  error-state  dynamics  matrix  F{t)  and  the  the  dynamics 
equation  noise  strength  matrix  Q(t).  Stacey’s  thesis  describes  the  LN-93  error  states 
in  more  detail  including  the  mathematical  models  of  first  order  Markov  processes, 
random  constants  and  random  walks  used  in  all  the  truth  models  and  filters  (22). 

3.3  Stacey’s  128~State  NRS  Afodel 

This  section  overviews  the  NRS  model  Stacey  developed  in  his  thesis  (22). 
Stacey  only  performed  a  single  Monte  Carlo  run  on  this  filter,  and  its  performance 
was  inconclusive.  This  thesis  analyzes  Stacey’s  128-state  filter  using  lO-run  Monte 
Carlo  analysis.  The  results  of  this  analysis  can  be  found  in  Chapter  V  with  plots  of 
filter  performance  in  Appendix  D. 

3.3.1  The  72-Siate  LN-93  Error  Model.  Stacey  used  the  Litton  LN-93  error 
state  model  as  the  basis  for  the  work  performed  in  his  thesis  (22).  However,  Stacey 
recognized  the  inadequacy  of  the  a  baro- altimeter  stale  in  this  model  and  designed  a 
4-state  baro-altimeter  error  model  (22)  to  implement  as  the  truth  and  filter  models  in 
his  thesis.  This  increased  the  INS  error  model  to  96  total  error  states.  Unfortunately, 
Stacey’s  four-state  baro-altimeter  model  never  performed  up  to  expectations,  and  the 
ENRS  designed  in  this  thesis  returns  to  a  single  baro-dtimeter  state  (22).  Tables 
A.l  through  A,4  in  Appendix  A  list  the  93  INS  error  states. 

When  Stacey  merged  his  96  state  INS  model  with  the  26  side  RRS  model 
and  the  30  State  GPS  model,  he  found  the  VAX  III  workstations  he  was  u.sing  did 
not  contain  enough  memory  to  compile  and  link  the  152-state  Fortran  code.  He 
knew  Solomon  had  successfully  executed  code  with  127  states,  so  he  reduced  the 
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INS  model  by  24  states,  maintaining  the  93-state  error  performance,  bringing  the 
total  error  model  to  128  states  (20,  22).  Stacey  eliminated  the  last  24  states  in  the 
INS  model  and  then  ran  simulations  to  prove  the  position,  velocity,  and  attitude 
errors  were  essentially  the  same  as  in  the  full  96-state  error  model.  At  this  point, 
Stacey  declared  this  72-state  INS  model  to  be  the  truth  and  filter  models  in  his 
analysis  (22). 

3.3.2  The  26-State  R RS  Error  Model.  The  Range/Range-Rate  System  (RRS) 
is  a  navigation  aiding  system  which  comprises  a  significant  part  of  CIRIS.  Conse¬ 
quently,  the  RRS  is  used  in  both  the  NRS  and  ENRS.  Navigation  information  is 
obtained  by  “interrogating”  ground  transponders  and  subsequently  processing  the 
electromagnetic  (EM)  signals  which  the  transponders  emit.  The  information  ob¬ 
tained  allows  high  quality  range  and  range-rate  measurements  to  be  calculated  by 
the  RRS  interrogating  hardware  (19).  Using  these  range  and  range-rate  measure¬ 
ments,  updates  to  the  NRS  position  (and  ENRS  position  and  velocity  estimates)  are 
then  possible. 

Solomon,  Snodgrass,  and  Stacey  each  have  given  detailed  explanations  of  the 
RRS  error-states  and  the  range  measurement  equation  (20,  19,  22).  This  section 
reviews  these  topics  to  acquaint  the  reader  to  the  dynamics  and  measurenient  equa¬ 
tions.  The  range-rate  measurement  is  explained  in  the  next  chapter  when  the 
ENRS  is  developed.  The  following  sections  havfe  been  excerpted  from  Chapter  4 
of  Stacey’s  thesis,  and  the  reader  should  review  his  thesis  if  a  detailed  explanation 
of  the  Range/Range-Rate  System  (RRS)  is  needed. 

3.3.2. 1  RRS  Model  Equations.  The  RRS  error  state  v«:tor  is  composed 
of  26  elements  (shown  in  Table  A.5,  Appendix  A).  The  RR.S  states  occupy  numbers 
through  -1-26  in  the  complete  error  state  model,  where  S, ^5  represents 
the  total  number  of  states  used  to  model  the  INS  subsystem  in  the  truth  model. 
The  first  two  RRS  states  are  simple  random  constant  (bias)  states  which  model 


the  effects  of  user  hardware  (RRS  iriterrogator)  range  and  range-rate  calibration 
errors,  respectively.  Coupling  of  these  states  into  the  RRS  measurement  equation  is 
presented  in  the  next  subsection.  The  error  state  model  equation  for  these  states  is 
(19): 


where 


Xir  =  range  equivalent  of  interrogator  bias 

xin  =  velocity  equivalent  of  interrogator  bias 


The  initial  state  estimates  and  covariances  for  these  states  are  (19): 


and 


Pfrr,6«(to)  — 


1/1*  0 
0  10“^/<*/sec* 


(3.4) 

(3.5) 


While  the  two  states  discussed  above  are  common  to  all  RRS  measurements, 
there  exist  two  sources  of  errors  which  are  unique  to  each  individual  transponder. 
First  is  the  error  due  to  transponder  surveyed  pasition  uncertainty  (x,  y,  z  compo¬ 
nents  in  earth-centered  earth  fixed  (ECEF)  frame  (2)),  and  second  is  the  error  due  to 
atmospheric  propagation  delays  between  the  user  and  each  individual  transponder. 
The  three  position  error  sources  are  well  modeled  by  tandom  bias  states  and  the 
atmospheric  error  states  are  represented  by  first  order  Markov  proces.ses.  Then  for 
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each  transponder,  four  states  axe  used  to  define  error  sources  (19): 


0  0  0 
0  0  0 
0  0  0 
0  0  0 


(:  6) 


The  subscript  i  in  the  equation  above  denotes  the  transponder  and  t  =  300  sec 
is  the  transponder  atmospheric  error  state  time  constant.  The  initial  conditions  for 
these  states  are  (19): 

^x,yrt,atm{^o)  —  0  (^•'^) 

25ft^  0  0  0 

0  25ft*  0  0 

P,.y.,.a<m(fo)  =  '  (3.8) 

0  0  25/t*  0 

0  0  0  lOO(PPM)* 

and 

E  =  0  (3.9) 

0  0  0  0 

^  <»<"»(  v.,.atm(H-r)}=r  ®  °  °  U(r)  (.-^.lO) 

0  0  0  0 

0  0  0  %^) 

with  =  10"*°/^*'  Equations  (3.3)  through  (3.10)  were  developed  by  the  Cubic 
Corporation,  designers  of  the  RRS  system  for  CIGTF,  and  are  based  on  real  static 
and  dynamic  measurement  analysis  of  the  RRS  (20).  Once  again,  the  set  of  equations 
above  (Equations  (3.6)  through  (3.10))  apply  to  a  single  transponder.  There  are  six 
such  sets  of  equations  fcr  RRS  transponders  which  are  used  in  this  thesis.  The 
complete  RRS  error  vector  is  specified  in  Appendix  A. 
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3.3. 2.2  RRS  Range  Measurement  Equation.  In  CIRIS,  RRS  range  mea¬ 
surements  aid  in  estimating  position  errors  of  the  reference  INS,  ^he  RRS  range 
measurement  is  derived  from  the  time  delay  detected  between  the  time  at  which  the 
reference  hardware  (CIRIS  or  NRS)  sends  an  interrogation  signal  and  the  time  at 
which  a  reply  from  the  transponder  is  received.  This  temporal  difference  is  multi¬ 
plied  by  the  speed  of  light  (and  divided  by  two,  to  account  for  the  “round  trip”  of 
the  EM  signal)  to  obtain  an  uncorrected  range  measurement.  Correction  factors  are 
then  applied  in  order  to  compensate  for  delays  introduced  by  the  propagation  of  EM 
signals  through  the  atmosphere  and  to  correct  for  errors  introduced  by  equipment 
calibration  biases  (19).  The  uncorrected  range  measurement  as  obtained  from  RRS 
is  modeled  as:  . 


Runs  ~  Rt  "¥  V  (3.11) 

where 


^HRS  ~ 

RRS  range  measurement,  from  transponder  to  user 

Rt 

true  range,  from  transpo’^de;  to  user 

^Ratm  ~ 

range  error  due  to  atmospheric  delay 

SFhr  = 

range  error  due  to  equipment  calibration 

V  = 

zero-mean  white  Gaussian  measurement  noise 

Equation  (3.11)  is  a  model  for  the  range  as  determined  by  the  RRS  truth 
model.  It  includes  the  true  range  (which  is  never  precisely  known)  along,  with  terms 
which  reflect  the  sources  of  error  discussed  in  the  previous  subsection. 

3. 3.2. 3  Range  Calculation  From  INS  Data.  In  order  to  formulate  a  dif¬ 
ference  measurement  as  shown  in  Chapter  II,  two  sources  of  range  information  must 
be  obtained.  The  first  is  the  RRS  range  measurement  which  is  modeled  by  Equa¬ 
tion  (3.11).  Another  range  indication  is  computed  from  the  INS  indicated  position 
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and  RRS  (surveyed)  positions.  (Both  the  indicated  INS  position  and  transponder 
surveyed  position  contain  uncertainties  which  must  be  considered.)  In  this  approach, 
the  user  (INS)  indicated  position  is  represented  by  an  vector  expressed  in  the 
Litton  ECEF  coordinates  as; 


(3.12) 


while  the  true  RRS  transponder  position  is  represented  in  the  Litton  ECEF  frame 
by: 


(3.13) 


Then  the  calculated  range  from  the  user  (INS)  to  the  transponder  is  given  by: 


R 


ms 


I  ^  j/r  I 

I  2t,  ]  ^  Zt  \ 


(3.14) 


Equation  (3.14)  may  be  rewritten  as: 


^ms  =  -  Xrf  HVv  “  Vrf  +  i^u  “ 


(3.15) 


If  the  nonlinear  equation  above  is  “perturbed”  to  reflect  uncertainties  in  user  and 
transponder  positions,  then  a  first-order  Taylor  series  may  be  written  to  approximate 
the  range  (2,  14).  The  tiuncated  first-order  series  is  of  the  form  (20); 
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After  substituting  Equation  (3.15)  into  Equation  (3.16)  and  evaluating  the  partial 
derivatives,  the  INS-derived  range  approximation  becomes  (19,  20): 
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At  this  point,  the  difference  measurement  may  be  formed  as: 
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(3.18) 


The  range  mecisurement  noise  v  has  been  determined  to  have  a  variance  of  4  in 
the  truth  model  when  RR5  range  measurements  occur  at  6  second  intervals  (20, 19). 
Note  that  the  true  whole- valued  range  (Rt)  formerly  nresent  in  both  individual  range 
me^urements  is  cancelled  in  the  differencing  operation  (2,  14).  Also  note  that  the 
bracketed  coefficients  in  the  equation  above  appear  in  the  H  matrix  presented  in  the 
EKF  update  equations  overviewed  in  Chapter  II.  Also  note  that  ihe  measurement 
equations  must  be  implemented  in  a  single  coordinate  frame,  the  Litton  ECEF  frame 
being  the  frame  of  choice  in  this  thesis  (11).  Although  not  presented  in  this  thesis, 
Capt  Stacey  does  an  outstanding  job  of  deriving  the  necessary  coordinate  frame 
transformations  to  implement  all  the  measurement  equations  in  his  thesis  (22). 

S.S.3  The  30-State  GPS  Error  Model.  GPS  is  designed  to  be  a  highly  accu¬ 
rate,  stand-alone  navigation  system.  However,  for  this  research,  GPS  is  used  as  a 
subsystem  to  improve  the  navigation  solution  of  the  LN-93  based  NRS.  In  a  manner 
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somewhat  reminiscent  of  the  RRS  transponder  system  discussed  in  the  previous  sec¬ 
tion,  GPS  navigation  information  is  obtained  from  EM  signal  propagation  through 
the  media  (space  and  atmosphere)  between  the  user  (NRS)  and  each  of  the  space 
vehicles  (SVs)  which  the  user  “locks”  into  a  reception  channel  of  the  GPS  receiver. 
In  a  stand-alone  GPS  receiver,  navigation  information  is  obtained  by  receiving  GPS 
SV  ephemeris  data  which  are  broadcast  continuously  from  each  active  ( “locked-on” ) 
SV,  correlating  the  phase  of  the  signal  with  a  matching  signal  in  the  GPS  receiver, 
and  correcting  for  known  error  sources  to  produce  a  highly  accurate  ratine  estimate 
,  between  the  user  and  each  SV  which  is  monitored.  In  this  thesis,  uncorrected  range 
measurements  (known  as  pseudorange  measurements)  are  channeled  to  a  Kalman 
filter  which  prov  des  estimates'of  the  error  sources.  Although  not  used  by  the  NPS, 
range-rate  information  may  be  obtained  from  GPS  ephemeris  in  a  similar  manner. 
As  in  the  RRS,  GPS  range  measurements  make  refinements  to  the  NRS  navigation 
solution  possible. 

A  basic  dynamic  error  model  for  the  GPS  system  was  developed  by  Solomon 
(20)  and  revised  in  (6).  However,  substantial  changes  to  the  basic  GPS  model  have 
been  made  by  Stacey  (22).  In  the  references  cited  (6,  20),  a  simplified  GPS  model 
was  assumed.  It  consisted  of  four  stationary  space  vehicles  (SVs)  and  did  not  per¬ 
form  geometry  optimization  calculations.  In  St^ey’s  model,  a  24-SV  “optimal” 
constellation  based  on  a  paper  by  Gteen  (7)  is  modeled.  Stacey’s  model  includes 
orbital  calculations  for  all  SVs,  and  simulates  GPS  receiver  operation  as  well.  These 
enhancements  are  discussed  in  detail  in  Chapter  V  of  Stacey’s  thesis  (22).  The  mea¬ 
surement  model  equations  for  the  GPS  system  follow  a  parallel  development  to  that 
of  the  RRS  measurement  model  shown  previously. 

3.3.3. 1  GPS  Error-Model  Equations.  The  GPS  error  state  vector  is 
.  composed  of  30  elements  (shown  in  Table  A.6,  Appendix  A).  The  GPS  states  oc¬ 
cupy  the  thirty  “uppermost”  states  in  the  NRS  error  state  model.  The  first  two 
GPS  states  model  the  user  set  (receiver)  clock  bias  and  drift  errors,  respectively. 
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xuclki,  =  range  equivalent  of  user  set  clock  bicis 

xucikjr  —  velocity  equivalent  of  user  set  clock  drift 


(3.20) 


(3.21) 


Note  the  large  uncertainties  associated  with  the  user  clock  states.  These  values 
are  representative  of  the  normal  receiver  clock  bias  and  drifts  associated  with  GPS 
receivers  used  in  military  and  civilian  airborne  applications.  Until  the  user  clock  error 
is  determined,  it  is  the  single  largest  source  of  error  in  GPS  range  measurements. 
While  the  two  states  discussed  above  are  common  to  all  GPS  measurements,  there 
exist  five  sources  of  errors  which  are  unique  to  each  individual  SV. 

One  error  source  specific  to  each  SV  is  code  range  quantization  error.  At  the 
heart  of  any  GPS  receiver  exists  a  pair  of  interacting  tracking  loops  (12).  One  of 
these  loops,  the  “code  tracking  loop”  is  the  source  of  pseudorange  error  which  is 
modeled  as  a  first  order  Markov  process  (12)  with  an  exponential  autocorrelation 
function.  Other  significant  error  sources  include  the  tropospheric  and  ionospheric 
propagation  delays.  Both  of  these  error  sources  are  modeled  and  corrected  to  a  large 
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degree  by  the  GPS  receiver.  However  the  uncompensated  error  contribution  of  these 
error  sources  is  still  significant.  Both  of  these  are  also  modeled  as  first  order  Markov 
processes  (with  different  time  constants).  Still  other  sources  of  error  which  must 
be  included  in  the  GPS  rnodel  are  SV  clock  error  and  SV  position  error,  which  are 
treated  as  random  bias  states.  The  reader  should  note  the  strong  possibility  for 
observability  problems  in  a  model  such  as  this.  In  this  case,  three  position  error 
states  are  used,  but  the  measurements  provide  new  information  only  along  the  line- 
of-sight  vector  between  the  user  and  SV.  All  error  sources  discussed  above  which  are 
unique  to  each  SV  are  included  in  an  error  state  vector  shown  below: 

SRci  ]  -1  0  0  0  0  0  0  [  1  Wci 

SRtrop  0  '  -555  0  0  ,  0  0  0  SRtrop  Wtrop 

SRion  0  O-j^OOOO  6Rion  Wian 


0  0 
0  0 
0  0 


0  0  0  0  0 
0  0  0  0  0 
0  0  0  0  0 
0  0  0  0  0 


0.25/t2  0  0 

0  1.0/f2  0 


^apsi^o)  — 


l.O  ft-  0  0  0 

0  25/t*  0  0 

0  6  25/<2  0 


0  25/t*  0 

0  0  25/f* 


(3.22) 


(3.23) 


E  {^apsi*)}  ^  0 


(3.24) 
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Hpwever,  two  major  problems  exist  in  using  this  simple  definition  for  range. 
First,  the  EM  signal  is  not  propagated  entirely  in  vacuum.  The  signal  originates 
in  space  (where  the  assumption  of  vacuum  is  acceptable),  but  must  subsequently 
propagate  some  distance  through  the  earth’s  atmosphere  as  well.  Naturally,  the 
signal  delay  introduced  by  atmospheric  propagation  must  be  taken  into  account. 
Second,  in  order  to  preserve  any  hope  of  accurately  determining  range,  it  is  critical 
to  determine  the  ii,M  signal  transit  time  with  an  extremely  high  degree  of  accuracy. 
Recalling  that  light  (or  any  EM  signal)  propagates  on  the  order  of  3  x  10*m/s,  it  is 
readily  apparent  that  even  a  very  small  error  in  determining  the  EM  signal  propaga¬ 
tion  time  can  wreak  havoc  on  attempts  to  use  uncorrected  pseudorange  information 
to  improve  the  navigation  system  solution  in  the  NRS  (or  any  other  such  system). 

As  a  consequence  of  the  concerns  above,  it  is  imperative  to  develop  a  much 
higher  fidelity  model  for  range  estimation.  A  typical  model  for  the  pseudorange 
measurement  between  the  user  and  space  vehicles  is  given  by  the  following  equation 
(20): 

^ars  =  /?«  +  6R,rop  +  Jr  SR^^,^  -p  SR^,^,^  +  v  (3.27) 

where 


^nps 

Rt 

SR,, 

S  Rfrep 

SR,am 


SRsdk 

SRi’dk 

V 


GPS  pseudorange  measurement,  from  SV  to  user 

true  range,  from  SV  to  user 

range  error  due  to  code  loop  error 

range  error  due  to  tropospheric  delay 

range  error  due  to  ionospheric  delay 

range  error  due  to  SV  clock  error 

range  error  due  to  User  clock  error 

zero-mean  white  Caiussian  measiirrrnent  noise 


As  in  the  comparable  e(|uation  for  RKS,  the  GPS  psetniorange  (*<|uation  alM)ve 
Includes  the  true  range  (which  can  never  be  known  exactly)  along  with  terms  which 
reflect  sources  of  error  and  uncertainty  inherent  to  GPS  range  measurements. 
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3.3. 3.3  GPS  Pseudorange  Calculation  From  L\'S  Data.  As  in  the  RRS 
subsystem,  it  is  desirable  to  formulate  a  difference  measurement  in  the  CPS  model 
as  well.  Once  again,  two  sources  of  range  information  must  be  obtained.  Like  the 
RRS  case,  the  first  source  is  the  range  measurement  which  comes  from  the  GPS 
subsystem  and  which  is  modeled  by  Equation  (3.27).  The  secoi.d  range  estimate  is 
constructed  by  differencing  IN.S-indicated  position  and  SV'  (broadcast)  positions  to 
calculate  the  range.  (Note  that  the  indicated  INS  position  and  SV  position  contain 
uncertainties  which  must  be  considered.)  The  user  (INS)  position  is  represented  by 
an  vector  expressed  in  the  Litton  ECEF  which  is  repeated  below  for  convenience. 
SV  position,  ,  is  represented  in  like  fashion. 


'  ' 

t 

II 

X 

Vv 

.  .  X,  =  . 

“V 

\  J 

Then  user-to-SV  pseudorange  may  be  calculated  as: 


Equation  (3.29)  may  also  be  written  in  the  equivalent  forip: 


(3.28) 


(3.29) 


-  \l(Xv  -  >*■5)*  +  (Vt.  -  y.^)*  +  (i,,  -  -v)*  (3.30) 

Invoking  perturbation  theory  (2,  14),  the  nonlinear  equation  above  is  written  as  a 
first-order  Taylor  series  to  approximate  the  (INS  derivetl)  calculation  of  user  to  SV 
range.  The  truncate<I  (to  first-order)  series  is  expres-sed  a.s: 
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When  Equation  (3.30)  is  substituted  into  Equation  (3.31)  and  the  partial  derivatives 
evaluated,  the  INS-derived  pseudorange  approximation  becomes: 
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Now  the  GPS  pseudorange  differtnce  measurement  is  formed  as: 
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The  pseudorange  measurement  noise  w  Stacey  used  has  a  variance  of  2  /<*  when 
GPS  pseudorange  measurements  occur  every  10  seconds.  It  should  also  be  noted 
that  the  97-state  NRS  truth  model  dso  has  slightly  higher  noise  modeled  in  its  GPS 
pseudorange  measurement  equation.  Stacey  had  implemented  a  measurement  noise 
variance  of  only  2  jO,  while  further  resewch  reveals  the  actual  noise  variance  to  be 


J  fi^  fee  this  tyi:.''  of  GPS  rP'-del  (12).  Of  course,  the  filter  measurement  noise  was 
correspondingly  intr-'ased  Xr,  noted  in  the  RRS  ceise,  the  true  whole-valued  range 
(Rt)  fcrmerly  prese  u  the  individual  pseudorange  representations  {R[f,s  ^.nd  R^ps) 
is  cancelled  •!;  t^e  t  ifF'^rencing  operation.  The  lu  rV  ted  coefficients  in  the  equation 
adove  are  u.'^ed  in  tne  H  !.iatri>  appear,  in  the  /iM’  update  equations  shown  in 
Chapter  Ih 

3.4  R^  heiion  to  a  LX  03  Error  ^*odel 

The  NRS  developed  by  Stacey  >  'S  states  in  both  iio  nth  and  filter 

modela  (22).  A  one-run  Monte  Carlo  a  .„h  sis  a  2-bour  flight  trajectory  of  Stacey’s 
NRS  using  MSOFE  on  a  G  MIP  VAX  worh.siation  takes  over  18  hours  to  complete. 
In  order  to  reduce  the  cc  i^cutationaJ  I  »’)den,  the  INS  truth  and  filter  model  in  the 
ENRS  is  further  reduced  torn  'he  72  .;tates  Stacey  used  down  to  the  41  states  fol¬ 
lowing  the  jyrcommendstioiis  of  Lewantowicz  and  Keen  (10).  This  31  state  reduction 
significantly  decea.ses  the  computatioiial  time  needed  for  analysis  of  the  NRS  and 
ENRS  filters,  but  does  not  significantly  affect  the  estimation  accuracy  of  the  errors 
of  interest  (position,  velocity,  and  attitude).  Tables  A. 7  and  A.8  in  Appendix  A  list 
the  41  states  used  in  this  new  INS  truth  model. 

It  is  necc8.sary  to  confirm  Lewantowicz  and  Keens’  claim  that  a  41-state  error 
model  gives  essentially  the  same  solutions  for  position,  velocity,  and  attitude  error 
as  the  full-order  truth  model  (10).  Their  methodology  for  arriving  at  41  states, 
based  on  their  68-state  truth  model  for  the  INS  which  had  already  eliminated  states 
69  through  93  from  the  93-state  truth  model,  was  to  eliminate  15  states  which 
contributed  one  or  more  orders  of  magnitude  less  error  contribution  to  the  attitude 
and  velocity  error  estimates  than  the  states  that  were  retained.  They  also  combined 
12  other  states  with  errors  which  exhibited  similar  behavior,  but  were  more  dominant 
(10).  This  thesis  takes  a  slightly  different  approach,  by  simply  eliminating  all  the 
states  they  recommend,  including  the  12  states  they  combined  with  other  states. 
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This  one  deviation  from  their  recommendations  doesn’t  appear  to  alfect  the  results 
adversely.  Tables  A. 7  and  A. 8  in  Appendix  A  contain  lists  of  the  41-state  INS  error 
model.  Appendix  C  gives  a  10-run  Monte  Carlo  comparison  between  the  93-state 
LN-93  truth  model  and  the  41-state  truth  model  discussed  above.  As  shown  by 
Figure  3.1  (and  Appendix  C),  the  41-state  model  does  give  position,  velocity,  and 
attitude  error-states  almost  exactly  matching  the  93-state  model  over  the  two-hour 
fighter  flight  profile.  At  this  point,  the  41  error-state  INS  model  is  considered  the 
truth  model  for  this  thesis. 

3.5  The  97-State  NRS  Error  Model 

By  merging  the  41-state  INS  model,  the  26-state  RRS  model,  and  the  30-state 
GPS  model,  a  new  97-state  NRS  model  is  formed.  This  NRS  model  performs  as  well 
as  Stacey’s  128  state  model  did,  but  does  so  3  times  faster  due  to  the  elimination  of 
31  INS  states.  Chapter  V  discusses  the  performance  of  Stacey’s  NRS  model  and  the 
97-state  NRS  model  and.  compares  their  performance  to  each  other  and  various  other 
filters  which  are  described  in  this  md  the  following  chapter.  Appendix  E  contains 
the  results  of  a  10-run  Monte  Carlo  ajialysis  performed  on  the  97-state  NRS  filter. 

3.5  The  69-State  NRS  Filter  Model 

Once  the  97-state  NRS  truth  and  full-order  filter  models  are  developed,  filter 
ordet  reduction  is  performed.  The  stab  -Iimension  for  the  ENRS  post-processing 
filter  is  limited  to  70  states,  this  number  was  set  as  the  target  for  a  NRS  reduced- 
order  filter.  Recognize  that  only  the  number  of  filter  states  is  to  be  reduced;  the 
truth  model  remdns  at  97  states. 
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Figure  3.1.  93  vs  41  State  INS  Models  (a)  Latitude  (b)  East  Velocity  and  (c)  East 
Tilt  Errors. 


93-State  INS  Model 

.... 

41-State  INS  Model 
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Using  research  accomplished  before  as  a  basis  (12,  5),  the  Kalman  filter’s  GPS 
error-model  weis  reduced  from  30  states  to  only  2  states  (10,  12).  The  two  remaining 
states  are  the  receiver  clock  bi<is  and  drift  states  (see  Appendix  A  or  GPS  error-state 
discussion),  leaving  a  69-state  NRS  filter.  Looking  at  the  plots  of  the  97-state  NRS 
filter  in  Appendix  E,  note  that  the  two  GPS  receiver  clock  error-states  are  many 
orders  of  magnitude  larger  thtin  the  eliminated  code  drift,  atmosphere,  SV  clock, 
and  SV  position  error-states.  This  huge  magnitude  difference  gives  an  indication 
as  to  why  these  states  can  be  eliminated  while  retaining  the  two  dominating  states 
(see  Chapter  IPs  filter  order  reduction  section).  However,  these  states  cannot  simply 
be  dropped  with  the  expectation  of  adequate  filter  performance.  The  filter  tuning 
techniques  of  Chapter  II  are  exploited,  noting  each  of  the  eliminated  states  is  part 
of  the  pseudorange  measurement  equation.  Thus,  the  measurement  equation  noise 
(see  Equation  (3.33))  is  increased  to  compensate  for  the  eliminated  states. 

As  stated  previously,  the  truth  model  measurement  noise  variance  was  9 
Chapter  II  recommends  2Klding  the  filter  covariance  values  of  the  eliminated  states  to 
the  measurement  noise  variance  to  compute  a  new  value  for  the  measurement  noise 
variance.  Performing  this  operation  gives  a  new  value  of  approximately  225  ft^  (9 
ft^  original  noise  -f  .09  code  loop  -I-  6.3  /<*•  atmosphere  -|- 16  /<*  SV  clock  -I-  64-3 
S V  position).  Remember,  this  new  value  of  the  measurement  noise  is  only  a  first  ■ 
best  guess,  further  tuning  may  be  needed  for  optimum  filter  performauice;  Chapter 
5  discusses  and  compares  the  69-state  NRS  filter  performance  to  the  other  filters 
analyzed  in  this  thesis.  Appendix  F  contains  the  plots  of  the  69-state  reduced-order 
NRS  filters  performance  against  the  97-state  truth  model. 

3. 7  Summary 

This  chapter  reviewed  the  LN-93  and  Stacey’s  NRS  error-state  models,  and 
then  developed  97-state  and  69-state  error  models  by  eliminating  31  INS  and  28 
GPS  states,  respectively.  These  three  NRS  filters  form  a  baseline  to  which  the 
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ENRS  filters  developed  in  the  next  chapter  are  compared.  Comparisons  of  filter 
performance  can  be  found  in  Chapter  V  while  the  state  variable  plots  are  located  in 
Appendices  D,  E,  and  F. 


IV.  The  Enhanced  Navigation  Reference  System  (ENRS) 

4.1  Overview 

This  chapter  give**  a  complete  description  of  the  ENRS.  First,  an  explanation 
of  DGPS  and  how  it  is  implemented  at  CIGTF  is  given.  Second,  a  DGPS  error 
model  is  developed.  Third,  the  DGPS  error  model  is  combined  with  the  41-state 
INS  and  26-state  RRS  models  from  Chapter  HI  to  form  an  89-state  ENRS  model  us¬ 
ing  only  DGPS  and  RRS  range  me2isurements  to  update  the  filter.  Fourth,  a  48-state 
reduced-order  ENRS  filter  is  developed.  Fifth,  the  DGPS  delta-range  measurement 
equation  is  derived  and  merged  into  a  22-state  reduced-order  DGPS  filter  (no  RRS) 
to  speed  the  tuning  of  delta-range  measurements.  Sixth,  the  RRS  range-rate  mea¬ 
surement  equation  is  overviewed  and  merged  into  a  46-state  CIRIS  filter  (no  DGPS) 
to  speed  tuning  of  the  range-rate  measurements.  Finally,  DGPS  delta-range  and 
RRS  range-rate  measurements  are  merged  into  the  48-state  reduced-order  ENRS 
filter  to  complete  the  development  of  a  post  processing  ENRS  filter  of  less  than  70 
states.  .  . 

4.2  Differential  Corrections  to  GPS  Measurements 

This  section  discusses  the  theory  behind  DGPS  as  it  is  is  being  implemented  at 
CIGTF.  Interinetrics,  Inc.  is  the  government  sponsored  contractor  responsible  for  the 
DGPS  reference  station  at  CIGTF  (1,  5).  Personal  interviews  were  conducted  with 
Mr  Darwin  Abbey  and  Mr  Scott  Dance  of  Intermetrics  to  determine  how  differential 
corrections  to  GPS  are  implemented  at  CIGTF.  The  following  discussion  of  DGPS 
comes  directly  from  these  interviews  and  the  DGPS  error  model  created  in  the  next 
section  is  i  combination  of  Intermetric’s  description  and  a  course  given  by  Navtech 
Seminars  on  DGPS  error  models  (1,5,  18). 

In  order  to  apply  differential  corrections  to  GPS  measurements,  a  ground  based 
reference  receiver  (GBR)  is  needed  as  well  as  the  airborne  GPS  receiver  (ABR)  used 


in  the  NRS.  Figure  4.1  shows  the  basic  DGPS  system  as  it  is  being  implemented 
at  CIGTF.  The  ground  based  receiver  must  be  capable  of  tracking  all  SVs  in  view, 
possibly  as  many  as  eleven  when  the  full  GPS  constellation  is  placed  in  orbit.  The 
GBR’s  antenna  position  has  been  permanently  fixed  and  surveyed  to  within  centime¬ 
ter  accuracy.  A  high  accuracy  rubidium  clock  is  used  in  place  of  the  GBR’s  normal 
clock,  greatly  decreasing  the  large  clock  errors  common  to  the  ABR  as  discussed 
in  Section  3.3.3.I.  Using  the  transmitted  SV  ephemeris  data,  the  GBR  computes 
its  position  (different  from  the  surveyed  position)  knbwn  as  the  ground  truth.  The 
data  (pseudorange  measurements,  SV  transmitted  corrections,  GBR  applied  correc¬ 
tions,  ground  truth,  clock  errors,  etc...)  is  fed  to  a  30386  personal  computer  which 
then  estimates  the  SV  position  errors,  SV  clock  errors,  and  atmospheric  delays  with 
great  precision.  Because  the  GBR’s  true  position  is  accurately  known  and  its  clock 
errors  are  much  smaller  than  the  ABR’s,  the  SV  position  and  clock  errors  and  atmo¬ 
spheric  delays  are  highly  observable.  This  is  unlike  the  NRS  GPS  model  where  the 
large  ABR’s  clock  errors  and  dynamics  cause  these  states  to  be  largely  unobserv¬ 
able.  These  errors,  now  called  differential  corrections,  are  time  tagged  and  stored  on 
magnetic  tape  or  disc.  Note  that  the  differential  corrections  could  be  iipmediately 
transmitted  via  a  data  link  fo  real-time  differential  corrections  if  the  need  ever  arises 

(1,  5). 

Remembering  that  the  NRS  is  a  post-processing  filter,  the  ro  Dairborne  pseudo¬ 
range  measurements  (which  are  also  stored  magnetically  and  time  tagged)  now  have 
the  differential  corrections  applied  before  they  are  analyzed  in  t  le  post-processing 
filter.  Of  course,  this  assumes  the  GBR  is  tracking  the  same  fmr  SVs  the  ABR 
was  tracking  (a  good  assumption  if  the  ABR  is  within  the  CIRI  i  test  range).  Us¬ 
ing  differential  corrections  in  this  manner,  the  SV  clock  error  Is  eliminated  from 
each  pseudorange  measurement  and  the  SV  position  errors  are  nearly  eliminated. 
Depending  on  the  distance  between  the  ABR  and  the  GBR  dur  ng  the  flight  pro¬ 
file,  the  atmospheric  propagation  errors  can  be  almost  totally  eliminated  for  close 
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Figxire  4.1,  GPS  Data  Analysis/Reference  Station 


trajectories  or  greatly  reduced  for  flights  within  200  miles  of  the  GBR.  Even  with 
long-range  DGPS  (flights  which  extend  more  than  200  miles  from  the  GBR),  the 
post-processing  in  the  30386  personal  computer  eliminates  some  of  the  atmospheric 
propagation  errors.  The  largest  remaining  errors  in  the  pseudorange  measurement 
are  the  ABR  clock  errors  (see  Section  4.2.2).  With  this  basic  knowledge  of  DGPS, 
a  new  DGPS  error  model  is  now  developed  assuming  that  differential  corrections 
have  previously  been  applied  tp  the  raw  pseudorange  measurements  from  the  ABR 
(1,5,18). 

4.2.1  The  22-State  DGPS  Model.  The  DGPS  error-model  is  composed  of  22 
states  (shown  in  Table  A.9,  Appendix  A).  Note  the  DGPS  states  are  almost  identical 
to  the  GPS  states  except  for  elimination  of  the  SV  clock  error  through  differential 
corrections,  and  removal  of  the  SV  code  loop  error.  Dance  states  that  compensation 
for  the  code  loop  error  (truth  and  filter  models)  is  provided  by  the  measurement 
noise  variance  in  the  pseudorange  measurement  as  implemented  in  this  thesis  (5), 
While  the  names  of  the  states  are  the  same,  only  the  receiver  dock  error  states  are 
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modeled  exactly  as  in  the  NRS.  These  two  error-states  are  again  modeled  by  (22): 

(4.1) 

(  iuclk^r  J  L  ®  ®  J  I  J 
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range  equivalent  of  ABR  clock  bias 
velocity  equivalent  of  ABR  clock  drift 


Again  the  initial  state  estimates  and  covariances  for  these  states  are  (22): 


2t/c/*fc(<o) 

XVclkiMo) 
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^Uclkt^UdkiMo)  = 


(4.2) 


(4.3) 
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As  noted  previously,  until  the  ABR  clock  error  is  determined,  it  is  the  single  largest 
source  of  error  in  DGPS  range  measurements.  While  the  two  states  discussed  above 
apply  to  all  DGPS  me2isurements,  The  troposphere,  ionosphere  and  SV  position 
errors  are  unique  to  each  SV.  Note  that  these  errors  have  a  much  smaller  contribution 
after  differential  corrections  than  in  the  GPS  model  (see  Section  3.3.3.1)  (18): 
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l.Oft^  0  0  0  0 

0  0  0  0 
0  0  .35/^2  0  0 

0  0  0,  .35/<2  0 

0  0  0  0  .35/<2 

and 

^{w^«Ps(0}  =  0  (4-6) 

0.001  0  0  0  0 

0  0.0004  0  0  0 

0  0  0  .0  0 

0  0  0  0  0 

0  0  0  0  0 

Once  again,  the  set  of  equations  above  apply  to  a  single  SV.  The  initial  covariance 
VcJues  in  Equation  (4.5)  and  atmospheric  error  dynamics  noise  variances  in  Equation 
(4.7)  were  obtained  from  (18)  and  are  modeled  slightly  larger  than  /\bbey  and  Dance 
(1,5)  recommend  since  a  conservative  filter  is  develcned.  There  are  four  such  sets  of 
matrix  equations  for  DGPS  SV  errors  modeled  in  this  thesis.  The  error-state  vector 
is  completely  specified  in  Appendix  A.  , 

4.2.2  DGFS  Pseudorange  Measuirment  Equation.  The  DGPS,  pseudorange 
measurement  equation  is  developed  identically  to  the  GPS  pseudorange  measurement 
equation  of  Section  3.3.3.2.  After  applying  differential  corrections,  the  measurement 
equation  is  modeled  as; 

^DGPS  ~  +  ^Rtrop  +  ^Rion  +  ^Rvdk  +  (4-8) 

where 


ift^lsec)  •  ^(r)  (4.7) 
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^DGPS  ~ 

DGPS  pseudorange  measurement,  from  SV  to  user 

Rt 

true  range,  from  SV  to  user 

S  Rtrop  ~ 

range  error  due  to  tropospheric  dela}' 

SRion  = 

range  error  due  to  ionospheric  delay 

SRuclk  = 

range  error  due  to  ABR  clock  error 

V  = 

zero-mean  white  Gaussian  measurement  noise 

As  in  the  comparable  equation  for  GPS  (Equation  (3.27)),  the  DGPS  pseudorange 
equation  above  includes  the  true  range  (which  can  never  be  known  exactly)  along 
with  terms  which  reflect  sources  of  error  and  uncertainty  inherent  to  DGPS  range 
measurements. 

4.2.2. 1  DGPS  Pseudorange  Calculation  From  INS  Data.  As  in  the  NRS 
subsystem,  it  is  desirable  to  formulate  a  difference  measurement  in  the  DGPS  model 
as  well.  Once  again,  two  sources  of  range  information  must  be  obtained.  Like  the 
NRS  case,  the  first  source  is  the  range  measurement  coming  from  the  DGPS  reference 
station  and  modeled  by  Equation  (4.8).  The  second  range  estimate  is  constructed 
by  differencing  INS-indicated  position  and  SV  (broadcast)  positions  to  calculate  the 
range.  Following  the  same  derivation  as  in  Section  3.3. 3.3,  the  DGPS  difference 
measurement  is  formed  as: 

Sz  =  Rjf/s  ■“  R-dgps 


■Sx^  - 

Vs-Vu" 

■hv 

i 

4  ^ 

1  \Rt.s\  J 

•Sx^  + 

’ys-y« 

•6y^  + 

.  l^/ATsl,  . 

-  -  (l]<5>Aon  -*  '•  ®  (4-9) 

The  pseudorange  me^lsu^ement  noise  variance  is  9  /t*  when  DGPS  pseudorange  mea¬ 
surements  occur  every  10  seconds.  As  nbted  in  the  NRS  case,  the  true  whole-valued 
range  (Rt)  formerly  present  in  the  individual  pseudorange  representations  {Rif,s  and 
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Rpops)  is  cancelled  in  the  differencing  operation.  The  brarket«»d  Coi'fficients  in  the 
equation  above  appear  in  the  H  matrix  u.sed  in  the  EKP'  update  txpiations  shown  in 
Section  2.2. 

4.3  The  89-Slate  E.\RS  Error  Model 

The  22-state  DGPS  error-model  developed  is  now  mergt“d  with  the  41 -state 
INS  model  and  26-state  RRS  model  to  form  an  Sfl-error-slate  K.\RS  model.  At  this 
point,  only  DGPS  pseudorange  and  RRS  range  measurements  are  available  to  the 
filter  for  position  updating.  A  lO-run  Monte  Carlo  analysis  is  performed  and  the 
plots  of  filter  j>erformance  are  located  in  Apf>endix  G.  This  filter  serves  a.s  a  ba.selinc 
for  DGPS  measurements,  and  a  discussion  of  its  |>erformance  relative  to  the  .N’RS 
filters  is  given  in  Chapter  V. 

4-4  The  43~State  E.\RS  Filter  (Position  Measurements  Only) 

Now  that  a  full-order  filter  has  been  di.scussed,  a  reduce<l-order  filter  of  70 , 
states  or  less  is  developed.  The  first  step  in  filler-order  reduction  uses  the  ntethod 
discus.sed  in  Sections  2.3  and  2.4.  Since  all  the  DGPS  error-states  are  small  except  for 
the  ADR’s  clock  erroM,  they  are  eliminated.  All  20  .states  eliminated  are  part  of  the 
pseudorange  equation,  so  again  the  measurement  noise  e  is  increased  to  compensate 
for  the  eliminated  states.  The  original  measurement  noise  variance  is  9fP,  and  the 
modified  no;  .e  variance  is. calculated  to  be  approximately  23JP.  Once  again,  this  is 
only  a  first  best  guc-is  and  additional  filter  tuning  may  l)e  niTessary  (s<’e  ('hapter 
V).' 

At  this  point  tlie  89- state  model  has  been  rerluced  to  69  states  by  the  elimina¬ 
tion  of  20  DGPS  stall's.  However,  DGPS  delta-range  and  RRS  range-rate  measure¬ 
ment  erjuations  are  to  l)e  implementt'dand  the  additional  computation  time  requirifl 
to  implement  4  DGPS  delta-range  and  6  RRS  range-rate  measurements  slows  the 
Monte  Carlo  analysi;  down  by  .'iO  percent  over  the  .NRS  69-state  filter  using  only 


position  aiding  measurements.  Further  order  reduction  is  necessary  to  ensure  the 
post- processing  filter  ran  complete  its  analysis  of  a  2-hour  flight  trajectory  within 
24  hours.  For  this  reason,  an  additional  21  INS  states  are  eliminated  to  reduce  the 
ENRS  filter  to  48  error-states.  Once  again  the  recommendations  of  f^wantowicz 
and  Keen  are  followed  for  this  INS  state  reduction  (10).  The  new  20-state  IN.S  filter 
is  formed  by  elimination  of  states  21  through  41  of  the  41-slate  INS  model,  and 
this  model  is  shown  in  Tables  A. 7  of  Appendix  A.  Zero-mean  white-gainssian  noise 
(lO-'^rad^/.src)  is  added  to  dynamics  ecjuations  of  the  first  two  l.NS  states  (related 
to  position  errors)  to  compensate  for  this  filter  order  reduction  (see  Section  2.4). 
The  plots  of  filter  performance  for  the  48  error-state  ENRS  filter  (no  velocity  aid¬ 
ing)  are  found  in  Appendix  H  with  a  discus.sion  of  filter  performance  relative  to  the 
previously  designed  filters  in  this  thesis  found  in  Chapter  V. 

DGPS  Dtlta-Rangt  Mfasunmfnts 

In  a  manner  similar  to  the  CPS  and  DCPS  pseudorange  measurement  equa- 
tion.s  development  in  Sections  3.3.3.2  an<I  4.2.2,  the  delta-range  measurement  equa¬ 
tion  ran  ahso  be  developed.  Typically,  the  velocity  between  the  us<»r  and  SV  is  not 
directly  measured  in  a  CPS  receiver.  Rather,  the  number  of  carrier-cycles  for  a 
period  of  lime  is  counted,  which  calculates  the  change  in  pseudorange  over  a  time 
interval,  not  the  time  rate  of  change  of  pseudorange  (12).  For  this  thesis,  since 
error-states  are  modeled  and  not  the  true  receiver  operation,  the  previous  pseudor¬ 
ange  measurement  (see  Fxpiation  (4.8))  is  subtracted  from  the  current  pseudorange 
mea.sureir"nt  and  divided  by  the  lime  interval  betw<*en  measurements  to  yield: 

+  (4.10) 


where 


^DGPs  ~  DGPS  delta-range  measurement,  from  SV  to  user 

Rt  —  true  delta-range,  from  SV'^  to  user 

SRvcik  =  delta-range  error  due  to  ABR  clock  drift  error 

V  =  zero- mean  white  Gaussian  measurement  noise 

Note  that  the  atmospheric  errors  in  Ek]uation  (4.8),  once  subtracted,  are  small  and 
are  ignored.  The  delta-range  equation  above  includes  true  delta-range  (which  is 
never  known)  along  with  an  ABR  clock  drift  error  term. 

4.5.1  DGPS  Delta-Range  Calculation  From  INS  Data.  As  in  the  GPS  subsys¬ 
tem,  it  is  desirable  to  formulate  a  difference  measurement  in  the  DGPS  delta-range 
model  as  well.  Once  again,  two  sources  of  delta-range  information  must  be  obtained. 
Like  the  pseudorange  case,  the  first  source  is  the  delta-range  measurement  which 
comes  from  the  DGPS  subsystem  and  which  is  modeled  by  Equation  (4.10).  The 
second  delta-range  estimate  is  constructed  by  taking  the  time  derivative  of  Equation 
(3.30).  As  a  reminder.  Equation  (3.30)  defined  Rj^^g  as: 

=  l(*v  -  +  (Vu  ~  Vs?  +  (4-11) 

and  taking  its  time  derivative  yields: 

=  .S,„r/Jiiss  (4.12) 

where: 

^iss  ~  ~  ^s)  iHv  ~  ys)(yi/  ~  y*)  4*  i^v  ~~  ^s)i^v  ~  *’s)I  (4-,13) 


\ 

\, 
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Remember  from  Section  3.3. 3. 3  that: 


Xy  =  [ii,  Vu  Za' 


(4.14) 


so  that: 


Xs  —  [^s  ys  ^sl 


Xy  =  [iy  ify 
Xs  =  (ij  i/s  ■*s] 


(4.15) 


(4.16) 


(4.17) 


Invoking  perturbation  theory  (2,  14),  Equation  (4.12)  is  written  as  a  first-order 
Taylor  series  to  approximate  the  (INS  derived)  calculation  of  user  to  SV  delta-range. 
The  truncated  (to  first-order)  series  is  expressed  as: 


=  R,  + 


dR,^,(X„Xy,XsXr)\ 


dR,„,(X„Xy,X„Xy)\ 


dR,^,(X„X^,X^.X,.)| 


\{XgXyXs<Xy)n 


dR,^g{Xg,X„,Xg,X,)\ 


\{Xg.XyXsXy)n. 


(4.18) 


Because  the  SV  fKwition  errors  are  modeled  as  bias  (constant)  terms,  their  time 
derivative  is  zero,  so  that  SX,  =  0.  When  Equation  (4.12)  is  substituted  into 
Equation  (4.18)  and  the  partial  derivatives  evaluated,  the  INS-derived  delta-range 
approximation  becomes: 

^INS  =  +  ((R;/y.v)('^l'  “  ~  (^/Ns*)('^'/ws)(*y  ~ 

+  ((«7N5)(i'e  -  i/s)  -  (^;^f)(^,ys)(ya  -  »5)]^ye 

+  ((«r,J5)(‘V  -  ^v)  -  -  ^s)JK 
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+  [i^lss)iyu  -  »s)]% 

+  [(-R^vsK^t^  - 

-  [iR7:,s)iiu  -  is)  -  iR;^i")iS,ss)i=^v  -  =^s)]6=^s 

-  [iR7L)iyu  -  iis)  -  (R7^i")iS,^,){y,  -  ys)]6ys 

-  Kss)(iv  -  is)  -  iR7^i^)iS.ss)(^a  -  ^s)]6^s  , 

'  (4.19) 

Now  the  delta-range  difference  measurement  is  formed  by  subtracting  Equation  (4.10)  from 
(4.19): 

Sz  =  R,^g  -  Rdgps 

+  ((^7^s)(yu  -  Vs)  -  iR7si')is>ss)iyv  -  ys)]^yc; 

+  [(«7;Js)(^c/  -  ^s)  -  iR7si")iS,ss)i^v  -  ^s)]6^v 

+  ((^7jis)(*«/ ■  *s)]^ii/ 

+  [(^7iis)(yi' -  ys)l% 

+  l(fi7:J5)("c-^s)15ic 

-  [(^7iy‘5K*v  “  *s)  -  iR7si^)iRiss)i^i;  "  *s)l^»s 

-  MsKyc.  -  y^)  -  («7^?)(^/^s)(yu  -  ys)]^ 

~  ((^/iJs)(^«/  ~  is)  “  iRi^s^)iRttfs)i^u  ~  ^s)]^*s 

-  [I]^fit/e4*  ~  ®  (4.20) 
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The  varicince  of  the  delta-r?nge  measurement  noise  v  used  is  0.16  when 
delta-range  measurements  occur  every  10  seconds  (12).  The  true  whole-valued  range- 
rate  (Rt)  formerly  present  in  the  individual  delta-range  representations  {R,f,s  and 
Rpcps)  is  cancelled  in  the  differencing  operation.  The  bracketed  coefficients  in  the 
equation  above  appear  in  the  H  matrix  of  the  EKF  update  equations  shown  in 
Section  2.2. 

4.5.2  The  22-State  DGPS  Filter  (No  RRS).  Now  that  the  DGPS  delta-range 
model  has  been  derived,  it  is  necessary  to  implement  this  algorithm  into  MSOFE 
to  analyze  its  performance.  To  speed  up  the  process  of  implementation,  a  22-state 
reduced-order  filter  is  modeled.  This  filter  contains  the  20  INS  states  discussed  in 
the  48-8tate  ENRS  filter  and  2  DGPS  clock  states.  ,  No  RRS  states  are  included  in  the 
filter  or  truth  model,  so  the  full-order  truth  model  for  this  system  is  63  states  (41  INS 
and  22  DGPS).  A  10-run  Monte  Carlo  analysis  using  the  2  hour  fighter  flight  profile 
on  this  filter  is  completed  in  5  hours,  allowing  two  or  three  analyses  per  day  until  the 
algorithm  is  properly  implemented  and  the  filter  well-tuned.  The  filter  performance 
plots  from  the  fighter  flight  profile  are  located  in  Section  I.l  of  Appendix  I.  For  a 
robustness  check,  the  filter  was  then  analyzed  using  a  10-run  Monte  Carlo  analysis  of 
a  2-hour  straight  flight  profile.  These  filter  performance  plots  are  located  in  Section 
1.2  of  Appendix  I.  A  discussion  of  the  results  of  this  filter  implementation  can  be 
found  in  Chapter  V. 

4.5  Transponder  Range- Rate  Measurements 

In  a  manner  similar  to  the  DGPS  delta-range  measurement  equation  develop- 
rhent  in  Section  4.5,  the  RRS  range-rate  measurement  equation  can  also  be  devel¬ 
oped.  From  Equation  (3.11),  the  RRS  range-rate  equation  is  modeled  by: 

=  ^«  +  +  V  (4.21) 
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where 


RRS  range-rate  measurement,  from  transponder  to  user 
true  range-rate,  from  transponder  to  user 
range-rate  error  due  to  calibration  bias  error 
zero-meau  white  Gaussian  measurement  noise 

Again  the  atmospheric  error  is  small  and  can  be  neglected.  The  delta-range  equation 
above  includes  true  range-rate  (which  is  never  known)  along  with  an  equipment 
calibration  velocity  error  term. 

4,6.1  RRS  Range-Rate  Calculation  From  INS  Data.  As  in  the  DGPS  sub¬ 
system,  it  is  desirable  to  formulate  a  difference  measurement  in  the  RRS  range-rate 
model  as  well.  Once  again,  two  sources  of  range-rate  information  must  be  obtained. 
Like  the  delta-range  case,  the  first  source  is  the  range-rate  measurement  which  comes 
from  the  RRS  subsystem  and  which  is  modeled  by  Equation  (4.21).  The  second  delta- 
range  estimate  is  constructed  by  t<iking  the  time  derivative  of  Equation  (3.15).  As 


a  reminder.  Equation  (3.15)  defined  Rif,s  as: 

Rtss  =  +  +  (4-22) 

and  taking  its  time  derivative  yields: 

Rins  ~  I'lHs  /  Rins  (4.23) 

where: 

I'jss  ~  [(®i;  “  *t)(*w)  "b  (ytf  —  yrK^v)  (^ir  ~ (^’24) 
Remember  from  Section  3.3.2.3  that: 

Xy  =  (X{,  jj,  (4.25) 

,Xy  =  [Xy  yj,  (4.26) 

4-13' 


^RRS  ~ 

Rt 

SRbv  — 

V  = 


The  transponders  do  not  move,  so  the  velocity  term  Xj.  =  0.  Invoking  perturbation 
theory  (2, 14),  Equation  (4.23)  is  written  as  a  first-order  Taylor  ^ries  to  approximate 
the  (INS  derived)  calculation  of  user  to  transponder  velocity.  The  truncated  (to  first- 
order)  series  is  expressed  as: 


When  Equation  (4.23)  is  substituted  into  Equation  (4.28)  and  the  partial  derivatives  eval¬ 
uated,  the  INS-derived  range-rate  approximation  becomes: 

=  Ri  +  [{R7ss)(iu)-iR7^i"){T,^s){^,-Xr)]6x^ 

+  [(iZ7^,)(j).,)-(/?r«TK^,;.s)(yv -itr)]^yv 

+  ((0(^v)  -  (^7v?)(^m5)(^v  - 

+  [(^7>JsKyv  “S'r)Kyv 

■b  [(^/;»5)(^tr  ~  ^t)]^'^V 

-  [(^7^.)(iv)  -  iR7^'s')iT:,s)(=^u  -  *r)]«ir 

-  P7;J5)(yv)  -  (fi7^T)(^/ivs)(yv  - 

.  -  ((0(^v)-(^7;vf  )(7’,^J(^,-z,)j^z^ 


(4.29) 

Now  the  range-rate  difference  measurement  is  formed  by  subtracting  Equation  (4.21)  from 
(4.29): 

dz  —  “  Rrrs 

,=  [(iZ7^,)(i  J 

+  [iR7ss)(yv)  -  (^7.?)(^ri.s)(2/.  -  Vryihu 

+  -  iR7^i^)iT:,s)i^u  - 

+  [(^7iisKj(v  - 

{{Rjns^^^v  ^r  )\^^u 

-  [(iZ7^,)(i„) -- )(T,^,)(i„ -- 

-  [iR7ss)iyu) - {R7s's){T,„styv - yr)¥yT 

-  -  V  (4.30) 


The  variance  of  the  range- rate  measurement  noise  v  is  0.09  /f*  when  range- rate 
measurements  occur  every  6  seconds  (19).  The  true  whole- valued  range  {RT)  formerly 
present  in  the  individual  delta-range  representations  {Rjffg  and  R„,)  is  cancelled  in 
the  differencing  operation.  The  br2w:keted  coefficients  in  the  equation  above  appear 
in  the  H  matrix  of  the  EKF  update  equations  shown  in  Section  2.2. 

4.6.2  The  46-State  CIRIS  Filter  (No  DGPS),  In  a  manner  similar  to  the  22- 
eiror-state  DGPS  filter,  a  46-error-state  CIRIS  filter  is  developed  to  implement  the 
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RRS  range-rate  meaisurement  algorithm.  The  filter  model  consists  of  the  20-state 
INS  model  and  26-state  RRS  model.  A  10-run  Monte  Carlo  analysis  using  the  2- 
hour  fighter  flight  profile  is  again  used  to  tune  the  filter,  and  the  filter  performance 
plots  axe  located  in  Section  J.l  of  Appendix  J.  As  a  robustness  check,  another  10- 
run  Monte  Carlo  analysis  using  a  2-hour  racetrack  flight  profile  is  performed.  The 
filter  performance  plots  of  this  analysis  are  located  in  Section  J.2  of  Appendix  J.  A 
discussion  of  this  filter’s  performance  is  found  in  Chapter  V. 

J^.l  The  48-State  ENRS  Filter  (Position  and  Velocity  Measurements) 

The  culmination  of  the  work  in  many  other  theses  is  fincdly  reached  when  the 
20-error-state  INS  modei,  the  26-state  RRS  model  with  range  and  range-rate  mea¬ 
surements,  and  a  2-state  DGPS  model  with  pseudorange  and  delta-range  measure¬ 
ments  are  assembled  into  the  48-error-state  ENRS  filter.  This  filter  is  tuned  using  a 
10-run  Monte  Carlo  analysis  using  the  fighter  flight  profile,  and  a  final  25-run  Monte 
Carlo  analysis  is  performed  to  demonstrate  this  filter’s  performance  relative  to  all 
the  other  filters  analyzed  in  this  thesis.  The  filter  performance  plots  are  found  in 
Appendix  K  while  a  complete  discussion  of  the  final  results  is  located  in  Chapter  5. 

4.8  Summary 

This  chapter  steps  through  the  development  leading  to  the  48-error-state  ENRS 
filter.  Differential  corrections  to  GPS  measurements  are  first  discussed,  followed  by 
development  of  a  DGPS  model.  Then  full-  and  reduced-order  models  and  delta-range 
and  range-rate  me^urement  equations  are  developed.  Finally,  all  this  data  is  merged 
into  a  48-state  post-processing  ENRS  filter  designed  to  determine  the  position  and 
velocity  of  the  test  article  accurately  so  its  performance  can  be  evaluated.  Remember, 
the  truth  model  throughout  this  chapter  is  composed  of  41  INS  error-states,  26  RRS 
error-states,  and  22  DGPS  error-states.  All  NRS  and  ENRS  filters’  performance  is 
discussed  in  detail  in  the  next  chapter. 
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V.  Results  of  Filter  Implementations 


5.1  Overview 

This  chapter  evaluates  the  performance  of  all  the  filters  implemented  in  this 
thesis.  Reduced-order  filters  are  compared  to  their  full-order  counterparts,  GPS 
filters  are  compared  to  DGPS  filters,  and  all  the  filters  are  analyzed  relative  to  the 
current  CIRIS  system.  As  a  final  presentation  of  results,  a  section  is  given  containing 
tables  of  the  filter  tuning  parameters.  Q(t)  and  R(t,)  if  different  from  their  truth 
model  values  for  all  the  filters  implemented  in  this  thesis.  Before  any  analysis  of 
filter  performance  can  begin,  an  important  assumption  pertaining  to  three  of  the 
truth/filter  models  in  this  thesis  must  be  explained  to  the  reader. 

5.1.1  The  Double  Precision  Factor.  The  truth/filter  models  implemented 
previously  by  Solomon,  Snodgrass,  and  Stacey  utilized  Fortran  on  VAX  mainframe 
and  workstation  computers  (20,  19,  22).  Their  computer  code  wais  normally  com¬ 
piled,  linked,  and  executed  using  “single  precision”  (the  7  most  significant  numbers 
are  stored).  For  Solomon  and  Snodgrass,  whose  work  dealt  mainly  with  implement¬ 
ing  an  updated  CIRIS  filter  without  GPS,  this  method  was  more  than  adequate. 
However,  when  Stacey  implemented  his  GPS  subsystem  and  brought  in  pseudor¬ 
ange  measurements,  all  his  filters  utilizing  GPS  diverged  and  could  not  be  tuned 
properly.  As  seen  in  Figure  5.1,  the  typical  truth  model  GPS  clock  bias  error  state 
attains  values  requiring  much  more  than  7  significant  figures  to  estimate  the  state 
adequately;  As  Stacey  had  mentioned  in  his  thesis,  “double  precision"  (14  most 
significant  figures)  is  required  for  optimal  filter  performance. 

Because  the  clock  bias  grows  to  such  a  large  number,  executing  in  single  pre¬ 
cision  mode  denies  the  filter  the  1-  or  2- foot  clock  bias  state  estimation  accuracy  re¬ 
quired  for  accurate  pseudorange  meas  urements.  However,  simply  switching  to  double 
precision  is  not  as  easy  as  it  seems.  When  Stacey’s  128-state  model  algorithms  were 
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Figure  5.1.  Typical  GPS  Clock  Bi£is  Error  State  for  Two  Hour  Flight  Profile 

modified  to  double  precision,  neither  the  VAX  mainframe  nor  workstation  computers 
had  enough  memory  available  to  link  the  code.  The  same  problem  arises  with  the 
97-State  NRS  and  89-State  ENRS  full-order  filter  models.  However,  all  the  reduced- 
order  filters  in  this  thesis  are  compiled,  linked,  and  executed  in  double  precision 
mode. 

At  this  point  an  assumption  is  made  so  that  the  full-order  filters  can  be  ana¬ 
lyzed,  and  the  reduced-order  filters  have  benchmark  filters  for  comparison  purposes. 
All  the  full-order  filters  truth  models  have  their  2  GPS  clock  states '  initial  covariance 
values  reduced  four  orders  of  magnitude  (see  Section  3.3.3  or  4.2.1).  By  reducing  their 
magnitudes,  the  filters  are  able  to  estimate  the  clock  bicis  error-state  adequately  in 
single  precision  mode,  and  these  2  clock  states  still  niaintain  errors  several  orders  of 
magnitude  larger  than  any  other  error  states.  Notice  in  Figures  D.8,  E.8,  and  G.8, 
that  ai  the  very  end  of  the  two-hour  ^.ight  profiles  the  clock  bias  state  approaches  the 
limit  of  single  precision  and  once  again  starts  to  diverge  slightly,  causing  the  clock 
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drift  state  to  diverge  also.  However,  overall,  these  full-order  filters  are  tuned  well 
enough  to  produce  accurate  navigation  solutions.  Note  that,  if  these  clock  states 
were  not  so  much  larger  in  magnitude  than  the  other  error  states,  this  approach 
v^ould  not  be  valid  and  the  results  could  not  be  accepted  without  serious  reserva¬ 
tions.  Also  note  that  the  CIRIS  filter  without  the  GPS  subsystem  does  not  need 
double  precision,  and  it  is  not  implemented  in  this  filter. 

In  summary,  the  128-  and  97-state  NRS  and  the  89-state  ENRS  truth  and 
full-order  filter  models  were  run  in  single  precision  mode  with  the  GPS  clock  error 
truth  states’  initial  covariance  values  reduced  4  orders  of  magnitude.  This  modifi¬ 
cation  in  no  way  affects  the  reliability  of  these  filters  when  used  as  benchmarks,  as 
the  2  GPS  clock  error-states  are  much  larger  than  any  other  error  states,  and  the 
GPS  clock  bias  state  continues  to  be  the  largest  source  of  errors  in  the  pseudorange 
measurement  by  several  orders  of  magnitude.  The  46-state  CIRIS  filter  is  also  run  in 
single  precision,  but  it  does  not  model  GPS  states,  so  double  precision  is  not  needed, 
All  of  the  reduced-order  NRS,  DGPS,  and  ENRS  filters  (evaluated  against  97-state 
NRS,  63-state  DGPS,  and  89-state  ENRS  full-order  double  precision  truth  models) 
are  programmed  in  double  precision  so  the  performance  of  these  filters  is  in  no  way 
compromised, 

5.2  The  NRS  Filters 

This  section  evaluates  the  performance  of  the  3  NRS  filters  analyzed  in  this 
thesis.  Remember  that  the  NRS  models  the  error  states  of  a  LN-93  INS  with  baro- 
altimeter,  RRS  transponders  with  range  measurernents,  and  GPS  receiver  with  pseu¬ 
dorange  measurements.  Stacey’s  128-State  NRS  filter  is  examined  first.  After  some 
preliminary  tuning  runs,  its  performance  is  outstanding.  Then  the  97-State  NRS 
filter  is  analyzed.  It  too  produces  a  navigation  solution  one  order  of  magnitude  bet¬ 
ter  than  CIRIS.  Finally,  a  69-state  reduced-order  NRS  filter  is  examined.  Although 
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not  as  impressive  as  the  first  two  NRS  filters,  its  performance  is  well  above  that  of 
CIRIS. 

5.2.1  Stacey’s  128^State  NRS  Filter.  As  the  largest  truth/filter  model  in  this 
thesis,  taking  over  8  days  to  complete  a  10-run  Monte  Carlo  simulation,  this  filter 
also  had  the  least  amount  of  tuning  performed  on  it  for  optimal  state  estimation.. 
Appendix  D  contains  the  plots  of  the  filters  performance  on  the  main  states  of  inter¬ 
est.  While  some  of  the  GPS  states  require  further  tuning,  overall,  this  filter  performs 
well,  with  the  possibility  of  a  significant  increase  in  precision  over  CIRIS  as  it  is  im¬ 
plemented  today.  Remembering  this  filter  was  implemented  in  single  precision,  Table 
5.1  contains  the  temporal  average  of  the  ensemble  average  of  true  filter  estimation 
errors  (Icr)  for  the  position,  velocity,  and  attitude  errors  over  the  two-hour  fighter 
flight  profile.  Ensemble  average  is  the  average  of  the  error-states  for  the  number  of 
Monte  Carlo  runs  performed  in  the  analysis  over  the  2-hour  trajectories.  This  vector 
is  the  (Time  in  the  filter  performance  plots  in  Appendices  D  through  K.  The  scalar 
temporal  average  is  simply  the  average  of  the  ensemble  average  for  the  2-hour  tra¬ 
jectories.  As  a  reminder  from  Chapter  1,  CIRIS  specified  accuracy  is  also  included 
in  this  and  all  tables  in  this  chapter  for  comparison  purposes.  Attitude  is  included 
in  this  and  future  tables  simply  to  round  out  the  navigation  solution  and  show  how 
INS  state  reduction  affects  these  error  states.  Remember,  no  attitude  measurements 
are  present  in  any  of  the  filters  implemented  in  this  thesis. 


Table  5.1.  Temporal  Average  of  the  Ensemble  Average  of  True  Errors  for  Stacey’s 
128  Erur-State,  Filter 


Filter 

Hal 

Longi-- 

tude 

(ft) 

Alti¬ 

tude 

(ft) 

East 

Vel 

(fps) 

Riil 

Up 

Vei 

(fps) 

East 

Tilt 

(arcs) 

North 

Tilt 

(arcs) 

Azi¬ 

muth 

(arcs) 

CIRIS 

14.00 

unia 

40.00 

0.100 

0.100 

0.400 

■Hil 

128-NRS 

1.59 

2.35 

5.58 

0.015 

0.011 

0.045 

1.65 

2.02 

12.83 

Stacey  never  resolved  the  precision  problem  for  this  filter  due  to  time  con¬ 
straints  and  wcis  only  able  to  perform  a  one-run  Monte  Carlo  simulation.  The  tun¬ 
ing  performed  in  this  thesis  on  his  filter  involved  only  the  random  bias  states  in  the 
RRS  and  GPS  subsystems.  RRS  range  bias  and  transponder  position  error  states 
and  the  GPS  SV  clock  and  SV  f>osition  error  states  have  a  dynamics  noise  strength 
Q{t)  addition  of  approximately  I0~^ftysec  fc.  the  performance  seen  in  Appendix 
D.  Although  not  yet  optimal,  the  filter’s  performance  at  this  point  speaks  for  itself. 
Using  only  RRS  range  and  GPS  pseudorange  measurements,  and  once  adequately 
tuned,  Stacey’s  filter  will  produce  a  navigation  solution  one  order  of  magnitude  more 
accurate  than  CIRIS. 

5.2.2  The  97-State  NRS  Filter.  This  filtei;  is  only  different  from  Stacey’s  filter 
by  use  of  the  41-state  INS  model  instead  of  Stacey’s  72-state  model  in  the  truth  and 
filter  models.  Since  the  4 1  -state  model  provides  true  error  states  which  almost  match 
those  for  the  93-state  LN-93  error  model  (see  Appendix  C),  this  filter’s  performance 
should  also  be  very  close  to  that  of  Stacey’s  128-state  filter.  This  simulation  was  also 
implemented  in  single  precision,  with  the'2  GPS  clock  error  states’  initial  covariance 
values  decreased  4  orders  of  magnitude.  Table  5.2  compares  the  128-  and  97-state 
NRS  filters  to  each  other  and  to  CIRIS. 


Table  5.2.  Temi>oral  Average  of  the  Ensemble  Average  of  True  Errors  for  128-  and 
97-State  NRS  Filters 


Filter 

Lati¬ 

tude 

(ft)  , 

lA>ngi- 

tude 

(ft) 

Alti¬ 

tude 

(ft) 

Ea.sl 

Vel 

(fps) 

.North 

Vel 

(fps) 

Up 

Vel 

(fps) 

East 

Tilt 

(arcs) 

North 

Tilt 

(arcs'; 

Azi¬ 

muth 

(arcs) 

CIRIS 

14.00 

14.00 

40.(K) 

0.100 

0.100 

0.400 

128-NHS 

1..59 

2.;r> 

.5..58 

0.015 

0.011 

0.015 

1.05 

2.02 

12.83 

j  97- NRS 

\.:\Ty 

2.71 

.5.28 

0.014 

0.010 

0.045 

1.07 

1.29 

9.71 

As  surmised,  the  j>erforn>ance  of  th'se  two  filters  is  almost  identical,  with  a 
large  increase  in  performance  over  the  (MKIS  sy  iem.  Any  di.screpancy  between 


Stacey’s  128-state  NRS  filter’s  performance  and  the  97-state  NRS  filter’s  perfor¬ 
mance  is  attributed  to  the  fact  that  more  time  was  spent  tuning  the  ‘jT-state  NRS 
filter  than  the  r28-state  filter.  Again,  the  only  filter  states  requiring  tuning  were 
the  random  bias  RRS  and  GPS  states  mentioned  in  the  previous-  section,  and  these 
states  have  increased  dynamics  equation  noise  similar  to  the  128-state  filter.  Ap¬ 
pendix  E  contains  the  filter  f>erformance  plots  of  the  main  states  of  interest.  Once 
again,  the  97-State  NRS  filter  provides  a  navigation  solution  one  order  of  magnitude 
more  accurate  than  the  current  CIRIS  system  while  using  only  RRS  range  and  GPS 
pseudorange  measurements. 

5.2.3  The  69-State  NRS  Filter.  This  first  attempt  at  a  rrrlucetl-order  filter 
shows  less  optima]  results  than  the  previous  two  NRS  full-order  filters.  Asa  reminder, 
the  69  error-states  implemente<i  in  this  filter  were  a  combination  of  41  INS  error- 
states,  26  RRS  error-states,  and  the  2  GPS  clock  error-states.  Also,  the  97-state  truth 
and  69-state  reduced-order  filter  models  ih  this  analysis  were  coded  and  executed  in 
double  precision,  so  the  GPS  clock  error  states’  initial  covariance  values  are  at  their 
true  magnitude.  The  filter  performance  plots  are  found  in  Appendix  F  while  'I’able 
5.3  compares  this  filter  with  the  other  NRS  filters  and  CIRIS. 


Table  5.3.  Temporal  Average  of  the  En.semble  Average  of  True  Errors  for  NRS 
Filters 
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l-CT 

5.28 
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o.ai5 

1.07 
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As  discussed  in  Section  3.6,  the  pseuflorange  measiirenient  no!s<’  variance  of 
this  filter  is  ailjusted  upward  to  comjwnsate  for  the  eliminated  GPS  error-states 
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(which  are  all  part  of  the  pseudorange  measurement  equation).  The  first  best  guess 
was  225  while  the  currently  implemented  value  of  300  ft^  is  still  not  enough 
to  ensure  adequate  performance  using  conservative  filter  tuning  methods  (see  the 
latitude  state  in  Figure  F.2).  Even  when  taking  into  account  the  filter's  nonoptirnal 
tuning,  it  appears  the  navigation  solution  from  this  filter  is  not  as  accurate  as  the  two 
full-order  filters.  However,  when  comparing  to  CIRIS  performance,  this  filter  is  at 
least  3  times  more  accurate  in  position  and  velocity  estimates.  This  is  a  substantial 
increase,  considering  only  RRS  range  and  GPS  ps<nidorange  equations  are  used  to 
update  the  filters  state  estimates.  It  can  be  deduced  from  this  filter  implementation 
that  the  eliminated  GPS  states,  even  though  they  have  relatively  small  magnitudes, 
do  affect  the  navigation  solution  accuracy  even  when  the  filter  is  optimally  tuned.  It 
is  with  this  thought  in  mind  that  one  is  lead  into  implementing  differential  corrections 
to  GPS  measurements  to  achieve  higher  accuracy  in  the  navigation  solution  when 
implementing  reduced-order  filters. 

5.3  The  ENRS  Filters 

Various  NRS  filters  have  l)een  implemented  and  analyzed  with  impressive  re¬ 
sults  compared  to  CIRIS.  The  incorporation  of  differential  corrections  to  GPS  pseu¬ 
dorange  mea.surements  takes  this  accuracy  a  step  further  with  the  ENRS.  The  89- 
State  ENRS  filter  is  analyzed  first,  Using  only  DGPS  p.seudorange  measurements, 
this  full-order  filter’s  performance  in  horizontal  position  error  standard  deviation  is 
on  the  order  of  I  ft.  Then  a  48-state  reduced-order  filter  is  examined.  Forty-one 
states  are  eliminated  with  only  a  slies  t  decrea.se,  in  velocity  and  attitude  error  esti¬ 
mation  performance.  Then  two  , iters  are  examined  which  implement  DGPS  delta- 
range  and  RRS  range-rat-  measurements  for  velocity  aiding.  Finally,  delta-range 
and  range-rate  mea.si  ements  an*  combined  into  the  48-state  ENRS  reducetl-order 
filter  as  the  conclusion  to  the  work  performed  in  this  thesis. 
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It  is  also  important  to  bring  out  at  this  time  that  when  an  ENRS  filter  is 
actually  implemented  at  CIGTF,  the  differential  corrections  should  ensure  that  the 
this  filter  performs  «dmost  as  well  in  real  life  as  it  does  in  simulation.  Because 
the  atmospheric,  SV  clock,  and  SV  position  GPS  error  sources  are  computed  very 
accurately  at  the  GPS  Reference  Station,  the  post-processing  filter  can  estimate  the 
true  position  errors  with  higher  accuracy  than  can  a  filter  utilizing  normal  GPS 
measurements.  The  added  error  sources  of  real  world  GPS  measurements  will  likely 
decrease  the  accuracy  of  the  NRS  filter  somewhat,  while  having  much  less  affect  on 
the  ENRS  filter  due  to  the  large  reduction  in  magnitude  of  the  error  sources  through 
differential  corrections  to  the  pseudorange  measurement. 

5.3.  t  The  89-State  NRS  Filter  (Position  Measurements  Only).  This  filter  is 
very  similar  to  the  97-state  NRS  filter,  except  for  the  assumption  that  differential 
corrections  have  been  applied  to  the  raw  pseudorange  measurements  before  filter 
processing.  It  includes  the  41-state  INS  model  and  the  26'.state  RRS  model  with 
the  GPS  clock  states’  initial  covariances  decreased  so  the  filter  executes  accurately 
in  single  precision  mode.  As  shown  by  Table  5.4,  the  inclusion  of  differential  correc¬ 
tions  improves  filter  performance  by  almost  50  percent  in  position  error  estimation 
accuracy  over  the  97-state  NRS  filler. 


Table  5.4.  Temporal  Average  of  the  Ensemble  Average  of  True  Errors  for  NRS  and 
89-State  ENRS  Filters 


Filter 

Lati¬ 

tude 

(ft) 

Longi¬ 

tude 

(ft) 

Alti¬ 

tude 

(ft) 

East 

Vel 

(fp.s) 

North 

Vel 

(fps) 

warn 

KBi 

East 

Tilt 

(arcs) 

North 

Tilt 

(arcs) 

Azi¬ 

muth 

(arcs) 

CIRIS 

14.00 

14.00 

40.00 

0.100 

0.100 

0.400 

128- NRS 

1.59 

mm 

5.58 

0.015 

0.011 

0.045 

1.65 

2.02 

12.83 

97- NRS 

1.35 

2.71 

5.28 

0.014 

0.010 

0.045 

1.07 

1.29 

9.74 

69- NRS 

3.28 

4.21 

9.04 

0.0.33 

0.026 

0.070 

1.37 

2.11 

18.42 

89- ENRS 

1.04 

3.80 

0.013 

0.010 

0.042 

0.94 

1.29 

6.46 
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The  filter  performance  plots  are  located  in  Appendix  G,  and  the  reader  will 
note  how  well  the  filter  is  tuned.  The  only  unusual  behavior  in  any  of  the  plots 
stems  from  the  fact  that  the  GPS  clock  bias  state  approaches  the  limit  of  accuracy 
of  single  precision  ne<u*  the  end  of  the  run  and  both  GPS  clock  states  diverge  slightly 
due  to  this  fact.  Other  than  this  small  anomaly,  this  filter  performs  better  than  any 
filter  examined  so  far. 

5.3.2  The  43-State  ENRS  Filter  (Position  Measurements  Only).  Since  the 
full-order  ENRS  filter  perfornis  so  well,  it  is  now  time  to  analyze  the  reduced-order 
filter  designed  to  fulfill  the  “less  than  70  state”  requirement  as  stated  in  Section 
1.2.  To  remind  the  reader,  48  states  were  chosen  over  the  69  states  of  the  NRS 
filter  because  the  addition  of  DGPS  delta-range  and  RRS  range-rate  metisurements 
increases  the  time  required  to  update  the  states  in  the  EKF  by  50  percent,  causing  the 
Monte  Carlo  analysis  to  take  almost  50  percent  longer.  Because  of  this  time  increase, 
21  additional  INS  states  relating  to  velocity  and  attitude  errors  were  eliminated 
along  with  the  20  upi>er  DGPS  error  states  as  discussed  in  Section  4.4.  Table  5.5 
once  again  compares  the  performance  of  this  48-state  reduced-order  ENRS  filter  to 
all  previous  filters  and  CIRIS.  As  can  be  seen,  the  position  states  suffer  almost  no 
degradation  due  to  the  order  reduction,  while  the  velocity  and  attitude  states  degrade 
significantly  from  the  eliminated  INS  states.  However,  the  addition  of  velocity  aiding 
measurements  should  increase  the  velocity  states’  accuracy  and  approach  the  full- 
order  filters’  performance. 

Appendix  H  contains  the  10-run  Monte  Carlo  analysis  plots  for  the  2-hour 
fighter  flight  profile.  Once  again,  using  the  techniques  of  Section  2.4,  this  filter  has 
been  conservatively  tuned  for  optimum  performance.  As  mentioned  in  Section  4.4, 
the  variance  of  the  pseudorange  measurement  equation  noise  v(/,)  was  increased 
from  its  full-order  value  of  9fP  to  30/t*  to  compensate  for  the  eliminated  DGPS 
states.  Also,  a  dynamics  equation  noise  of  \Q~^^ra<Pj sec  is  added  to  the  first  two 
INS  states  to  compensate  for  the  21  eliminated  states.  It  is  of  interest  to  note 
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Table  5.5.  Temporal  Average  of  the  Ensemble  Average  of  True  Errors  for  NRS, 
89-State  ENRS,  and  48-State  ENRS  Filters 


Filter 

Lati¬ 

tude 

(ft) 

Longi¬ 

tude 

(ft) 

nisi 

■Ql 

nsi 

Up 

Vel 

(fps) 

East 

Tilt 

(arcs) 

North 

Tilt 

(arcs) 

Azi¬ 

muth 

(arcs) 

CIRIS 

14.00 

14.00 

40.00 

0.100 

0.100 

liSiilil 

128-NRS 

1.59 

2.35 

5.58 

0.015 

0.011 

0.045 

1.65 

2.02 

12.83 

97-NRS 

1.35 

2.71 

5.28 

0.014 

0.010 

0.045 

1.07 

1.29 

9.74 

69-NRS 

3.28 

4.21 

9.04 

0.033 

0.026 

0.070 

1.37 

2.11 

18.42 

89-ENRS 

0.84 

1.04 

3.80 

0.013 

0.010 

0.042 

0.94 

1.29 

6.46 

48-ENRS-P 

0.90 

1.32 

3.05 

0.027 

0.020 

0.044 

3.83 

4.09 

19.23 

that  the  21  eliminated  INS  states  actually  feed  into  the  velocity  and  attitude  error 
states,  but  these  states  need  no  increase  in  noise,  while  the  position-related  states 
do  need  noise  aulded  for  acceptable  tuning.  Once  again,  this  48-state  reduced-order 
ENRS  filter  with  DGPS  pseudorange  and  RRS  range  measurements  increases  the 
navigation  solution  accuracy  in  position  by  better  than  one  order  of  magnitude 
compared  to  CIRIS.  Now  the  velocity  aiding  filters  are  explored  to  analyze  how 
velocity  measurements  can  increase  INS  velocity  error-state  estimation  accuracy. 

5.^  Velocity-Aiding  Filters 

The  filters  analyzed  in  this  section  implement  velocity-aiding  measurements 
into  the  DGPS  and  RRS  subsystems.  By  implementing  the  DGPS  delta-range  and 
RRS  range-rate  equations  separately  in  this  way,  each  measurement  is  separately 
tuned  for  optimal  performance  before  they  are  merged  together  into  the  final  48- 
state  ENRS  filter.  Also,  separating  the  DGPS  and  RRS  error  models  allows  faster 
execution  of  the  simulations  so  that  more  tuning  runs  can  be  performed  in  a  shorter 
period  of  time.  All  of  these  simulations  operate  in  double  precision  mode  using  a 
10-run  Monte  Carlo  analysis. 
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5.4-1  The  22-State  DGPS  Filter.  This  section  analyzes  the  results  of  the  22- 
state  DGPS  filter  used  for  DGPS  delta-range  measurement  equation  implementation. 
Table  5.6  contains  three  entries  for  the  22-state  DGPS  filter.  The  first  entry  “22- 
DGPS-P”  displays  the  average  errors  for  the  22-state  filter  using  only  pseudorange 
measurements  during  a  10-run  Monte  Carlo  fighter  flight  profile  analysis.  No  plots 
axe  included  in  the  appendixes  for  this  filter;  it  is  evaluated  simply  to  establish  a 
baseline  to  show  what  the  true  navigation  solution  errors  are  without  delta-range 
measurements.  The  second  entry  for  this  filter,  “22rDGPS-l”,  is  a  10-run  Monte 
Carlo  fighter  profile  analysis  incorporating  DGPS  delta-range  measurements.  The 
third  entry,  "22-DGPS-2’’  utilizes  a  straight  flight  profile  in  its  10-run  analysis. 


Table  5.6.  Temporeil  Average  of  the  Ensemble  Average  of  True  Errors  for  Previous 
and  22-State  DGPS  (Velocity  Aiding)  Filters 


Filter 

Lati¬ 

tude 

(ft) 

Longi¬ 

tude 

(ft) 

Alti¬ 

tude 

(ft) 

East 

Vel 

(fps) 

North 

Vel 

(fps) 

Up 

Vel 

(fps) 

East 

Tilt 

(arcs) 

North 

Tilt 

(arcs) 

Azi¬ 

muth 

(arcs) 

CIRIS 

14.00 

14.00 

40.00 

0.100 

0.100 

0.400 

],28-NRS 

1.59 

2.35 

5.58 

0.015 

0.011 

0.045 

1.65 

2.02 

12.83 

97-NRS 

1.35 

2.71 

5.28 

0.014 

0.010 

0.045 

1.07 

1.29 

9.74 

69-NRS 

3.28 

4.21 

9.04 

0.026 

0.070 

1.37 

2.11 

18.42 

89-ENRS 

0.84 

1.04 

3.80 

0.013 

0.010 

0.042 

0.94 

1.29 

6.46 

48-ENRS-P 

0.90 

1.32 

3.05 

0.027 

0.020 

0.044 

3.83 

4.09 

19.23 

22-DGFS-P 

1.95 

1.84 

3.14 

0.037 

0.044 

5.28 

3.94 

29.85 

22-DGPS.l 

2.10 

1.96 

3.25 

0.051 

5.46 

4.33 

45.02 

22-DGPS-2 

2.20 

2.14 

3.35 

0.049 

0.060 

5.56 

4.09 

63.24 

The  filter  performance  plots  for  the  delta-range  DGPS  filters  are  located  in 
Appendix  I  (22-DGPS-l  and  -2).  As  seen  from  the  plots  and  the  values  in  Table  5.6, 
the  DGPS  delta- range  measurement  is  corrupting  the  position  and  velocity  error- 
state  estimation  in  these  filters  a  small  amount  (remember  position  errors  are  part 
of  the  delta^range  measurement).  This  degradation  in  filter  performance  is  likely  to 
be  caused  by  cross-correlations  of  the  position  and  velocity  measurements  which  are 
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not  modeled  in  this  thesis,  and  it  is  currently  being  investigated  to  solve  this  problem. 
However,  even  with  this  anomaly,  this  reduced-order  filter  performs  as  well  or  better 
than  the  69-state  reduced-order  NRS  filter  (a  reduced-order  DGPS  only  filter  versus 
a  reduced-order  GPS/RRS  filter  issue).  Assuming  the  anomaly  with  the  DGPS 
delta-range  measurement  equation  is  remedied,  this  22-state  DGPS  filter  is  expected 
to  produce  position  and  velocity  error  estimates  50  percent  more  accurate  than  the 
69-state  NRS  filter  and  on  par  with  the  other  full-order  filters  discussed  previously. 
The  DGPS  delta-range  measurement  is  implemented  exactly  as  discussed  in  Section 
4.5.1,  with  no  changes  required  in  state  dynamics  noise  or  measurement  noise  from 
the  truth  model  values.  This  is  not  unusual  since  no  measurement  states  have  been 
eliminated,  essentially  a  “full-order”  measurement  equation.  Since  the  DGPS  delta- 
range  and  RRS  range-rate  measurements  are  very  similar,  it  is  no  surprise  if  the 
46-state  CIRIS  filter  analyzed  in  the  next  section  suffers  from  the  same  anomaly  as 
the  22-state  DGPS  filter. 


5.4-2  The  46-State  CIRIS  Filter.  Again,  this  filter  is  implemented  to  evaluate 
RRS  range-rate  measurements  optimally  before  they  are  incorporated  with  DGPS 
delta-range  measurements  in  the  48-state  ENRS  filter.  i  Table  5.7  contains  values 
for  three  CIRIS  filters  similar  to  the  DGPS  filters  discussed  in  the  previous  section. 
The  first  filter,  “46-CIRIS-P”,  is  a  46  state  filter  incorporating  only  RRS  range 
measurements.  It  is  used  as  a  baseline  to  evaluate  the  pejrformance  of  the  range-rate 
measurement  equations.  No  plots  are  contained  in  apper  dices  for  this  10-run  Monte 
Carlo  fighter  profile  analysis  since  filters  utilizing  only  RRS  range  measurements 
have  been  thoroughly  covered  in  previous  thesis  (20,  19,  22).  The  results  of  the 


second  CIRIS  filter,  “46-CIRIS-l”,  are  derived  from  a 


profile  analysis  using  RRS  range  and  range-rate  measurements.  The  third  filter. 


“46-CIRIS-2”,  was  analyzed  with  10-run  analysis  using  a 


0-run  Monte  Carlo  fighter 


racetrack  flight  profile  and 


also  incorporates  RRS  range  and  range-rate  measurements. 


Table  5.7.  Temporal  Average  of  the  Ensemble  Average  of  True  Errors  for  Previous 
and  46-State  CIRIS  (Velocity  Aiding)  Filters 


Filter 

Longi- 

I9B 

East 

North 

Up 

East 

North 

Azi- 

tude 

Vel 

Vel 

Vel 

Tilt 

Tilt 

muth 

(ft) 

(ft) 

1^ 

(fps) 

(fps) 

(fps) 

(arcs) 

(arcs) 

(arcs) 

CIRIS 

14.00 

14.00 

40.00 

0.100 

0.100 

0.400 

. . . . . I 

128-NRS 

1.59 

2.35 

5.58 

0.015 

0.011 

_ 

0.045 

1.65 

2.02 

12.83 

97-NRS 

1.35 

2.71 

5.28 

0.045 

1.07 

1.29 

9.74 

69-NRS 

3.28 

4.21 

9.04 

0.033 

0.026 

0.070 

1.37 

2.11 

18.42 

89-ENRS 

0.84 

1.04 

3.80 

0.013 

0.010 

0.042 

0.94 

1.29 

6.46 

48-ENRS-P 

0.90 

1.32 

3.05 

0.027 

0.020 

0.044 

3.83 

4.09 

19.23 

22-DGPS-P 

1.95 

1.84 

3.14 

isa 

0.044 

5.28, 

3.94 

29.85 

22- DGPS- 1 

2.10 

1.96 

3.25 

0.051 

InuKsi 

5.46 

4.33 

45.02 

22-DGPS-2 

2.20 

2.14 

3.33 

Kl 

0.049 

0.060 

5.56 

4.09 

63.24 

46-CIRIS-P 

3.12 

6.84 

18.09 

0.026 

0.100 

4.01 

4.39, 

47.48 

46-CIRIS-l 

2.80 

11.87 

B 

0.024 

0.046 

4.01 

4.34 

45.53 

46-CIRIS-2 

3.91 

6.20 

8.94 

■i 

0.026 

0.044 

4.08 

4.13 

12.92 

Appendix  J  contains  the  filter  performance  plots  for  the  filters  utilizing  range 
and  range-rate  measurements  (46-CIRIS-l  and  -2).  In  this  filter,  the  range-rate 
mecisurements  are  not  corrupting  the  position  states,  but  they  do  not  improve  the 
velocity  states’  estimation  except  possibly  in  the  vertical  channel.  In  this  filter,  the 
velocity  measurement  anornaly  does  not  have  a  pronounced  effect  as  in  the  DGPS 
filter,  possibly  because  the  transponders  are  stationary  and  the  GPS  satellites  circle 
the  earth.  Ih  any  case,  this  anomaly  still  hinders  the  evaluation  of  how  much  velocity 
aiding  will  enhance  the  navigation  solution.  The  RRS  range- rate  measurement  is 
implemented  exactly  as  discussed  in  Section  4.6.1,  with  no  changes  required  in  state 
dynamics  noise  or  measurement  noise  from  the  truth  model  values.  Again,  this  is  not 
unusual  since  no  measurements  states  have  been  eliminated,  essentially  producing  a 
“full-order”  measurement  equation.  For  comipleteness  and  in  the  hopes  the  anomaly 
is  eventually  corrected,  both  DGPS  delta-range  and  RRS  range-rate  measurement 
equations  are  merged  into  the  48-state  ENRS  filter  in  the  next  section. 


" 


5.5  The  Complete  48-State  ENRS  Filter 

Because  the  velocity  aiding  measurements  do  not  cause  the  DGPS  or  CIRIS 
filters  to  diverge  or  become  unstable,  a  decision  is  made  to  proceed  with  incorporating 
DGPS  delta-range  and  RRS  range-rate  measurements  into  the  48-state  ENRS  filter. 
This  way,  when  the  anomaly  is  discovered,  the  Fortran  code  can  be  quickly  fixed 
and  filter  analysis  preformed.  Table  5.8  contains  the  true  errors  for  the  navigation 
solution  for  all  filters  discussed  and  implemented  in  this  thesis.  Notice  that  even  with 
velocity  measurement  corruption,  the  48-state  ENRS  filter  performs  much  better 
than  CIRIS  or  even  the  69-statc  NRS  filter.  In  fact,  its  position  error-state  estimation 
is  one  of  the  best,  even  with  the  anomaly  present. 


Table  5.8.  Temporal  Average  of  the  Ensemble  Average  of  True  Errors  for  All  Filters 


Filter 

Lati- 

Ix)ngi- 

Alti- 

East 

Up 

East 

North 

Azi- 

tude 

tude 

tude 

Vel 

Vel 

Tilt 

Tilt 

nrtuth 

(ft) 

(ft) 

(ft) 

(fps) 

(fps) 

(fps) 

(arcs) 

(arcs) 

(arcs) 

CIRIS 

14.00 

14.00 

40.00 

0.100 

0.100 

0.400 

128-NRS 

1.59 

2.35 

5.58 

0.015 

0.011 

0.045 

1.65 

2.02 

12.83 

97-NRS 

1.35 

2.71 

5.28 

0.014 

0.010 

0.045 

1.07  , 

1.29 

9.74 

69-NRS 

3.28 

4.21 

9.04 

0.033 

0.026 

0.070 

1.37 

18.42 

89-ENRS 

0.84 

1.04 

3.80 

0.013 

0.010 

0.042 

0.94 

6.46 

48-ENRS-P 

0.90 

1.32 

3.05 

0.027 

0.020 

0.044 

3.83 

4.09 

19.23 

22-DGPS-P 

1.95 

1.84 

3.14 

0.044 

5.28 

3.94 

29.85 

22-DGPS-l 

2.10 

1.96 

3.25 

0.051 

5.46 

.4.33 

45.02 

22-DGPS-2 

2.20 

2.14 

3.33 

0.049 

0.060 

5.56 

4.09 

63.24 

46-CIRI3-P 

3.12 

6.84 

18.09 

0.026 

IQQI 

4,39 

47.48 

46-CIRIS-l 

2.80 

11.87 

0.024 

KfilB 

4.34 

45.53 

46-CIRIS-2 

3.91 

6.20 

8.94 

0.039 

0.026 

4.08 

4.13 

12.92 

48-ENRS 

1.19 

'  1.80 

5.46 

0.040 

0.029 

0.087 

5.14 

4:92 

31.21 

Plots  of  this  filter’s  state  estimation  performance  can  be  found  in  Appendix 
K.  These  plots  represent  a  25-run  Monte-Carlo  analysis  of  the  fighter  flight  profile. 
No  changes  were  made  to  the  dynamics  equation  or  measurement  noises  other  than 
those  discussed  in  previous  sections.  It  appears  that  when  the  velocity  measurement 


5-14 


equations  anomaly  is  remedied,  this  48-state  ENRS  filter  will  produce  navigation 
solutions  one  order  of  magnitude  more  accurate  than  current  CIRIS  system,  and 
will  rival  any  of  the  full-order  filters  implemented  or  analyzed  in  this  thesis. 


5.6  Filter  Tuning  Parameters 

This  section  contains  tables  of  the  filter  tuning  parameter’s  (Q(i)  and  R(<,)) 
final  values  for  the  filters  implemented  in  this  thesis  when  the  values  were  dif¬ 
ferent  from  the  truth  model  values.  Table  5.9  lists  the  Q{t)  values  from  among 
the  first  20  INS  error-states.  The  listed  states  were  the  only  states  requiring  tuning 
when  a  20-state  INS  filter  model  was  implemented,  and  the  tuning  values  were  the 
same  for  all  the  filters.  For  a  definition  of  the  listed  states,  consult  Table  A. 7  in 
Appendix  A. 


Table  5.9.  INS  States  Tuning  Parameters  Q(<) 


State 

State 

Truth  Model 

Filter 

Number 

Symbol 

Value 

Value 

1 

60, 

0.0 

lO~‘®ra</^/i(cc 

2 

60y 

0.0 

\Q~^^ra^  !  sec 

12 

6S3 

0.0 

lO'^fPfsec^ 

1? 

6Si 

0.0 

lO-^fPIsec 

14 

,  V,, 

0.0 

10~'^rad^/sec 

15 

Vv. 

0.0 

I0~^^rad^  fsec 

16 

0.0 

l0"*^ra(/^/sec 

17 

6gx 

0.0 

10~^^ug^ /sec 

18' 

69y 

0.0 

I0~^^ug^/sec 

19 

6gz 

0.0 

I0~^^ug^/sec 

Table  5.10  lists  the  26  RRS  error-state’s  tu^-iTi^  parameters.  All  of  the  RRS 
error-states  required  tuning  and  are  listed.  Remciriber  that  there  are  range  and 
velocity  bias  states  along  with  four  erre  states  (3  position  errors  and  1  atjuospheric 
error)  for  each  of  six  transponders  in  the  truth  and  filter  models  throughout  this 


thesis.  Table  A.5  in  Appendix  A  contains  the  definition  for  each  of  the  states  now 
listed. 


Table  5.10.  RRS  States  Tuning  Parameters  Q(<) 


State 

State 

Truth  Model 

Filter 

Number 

Symbol 

Value 

Value 

1 

SRt, 

0.0 

l0~^fPlsec 

2 

Svb 

0.0 

W^JPJs^ 

3,  7,11,15,19,23 

SPn. 

0.0 

.  2-lO-^fPlsec 

4,  8,12,16,20,24 

SPn^ 

0.0 

2  •  10~^fP/sec 

5,  9,13,17,21,25 

SPn, 

0.0 

2  •  fP ! sec 

6,10,14,18,22,26 

6.66  •  l0~^^fPjsec 

6.66  •  IQ'^fPIsec 

Table  5.11  lists  the  30  GPS  error-state’s  tuning  parameters.  Most  of  the  GPS 
error-states  required  tuning,  so  they  are  all  listed.  Listed  are  receiver  clock  bias 
and  drift  error-states  along  with  seven  error-states  (1  code-loop,  2  atmosphere,  1 
SV  clock,and  3  position  errors)  for  each  of  ,4  SVs  in  the  full-order  filter  models 
throughout  this  thesis.  The  reduced-order  GPS  filters  only  include  the  two  clock 
error-states  but  are  still  tuned  with  the  values  shown  below  for  these  two  states. 
Table  A. 6  in  Appendix  A  contains  the  definition  for  each  of  the  states  now  listed. 

Table  5.12  lists  the  22  DGPS  error-state’s  tuning  parameters.  All  of  the  DGPS 
error-states  required  tuning,  so  they  are  all  listed.  Listed  are  receiver  clock  bias  and 
drift  error-states  along  with  five  error-states  (2  atmosphere  and  3  position  errors) 
for  each  of  4  SVs  in  the  full-order  filter  models  throughout  this  thesis.  The  reduced- 
order  DGPS  filters  only  include  the  two  clock  error-states  but  are  still  tuned  with 
the  values  shown  below  for  these  two  states.  Table  A.9  in  Appendix  A  contains  the 
definition  for  each  of  the  states  now  listed. 
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Table  5.11.  GPS  States  Tuning  Parameters  Q(<) 


State 

'I'ruth  Model 

Filter 

Number 

Value 

Value 

1 

0.0 

mnwwmm 

2 

0.0 

I0~^ft^fsec 

3,10,17,24 

^Rcloopi 

0.5  ft^fsec 

O.Sft^fsec 

4,11,18,25 

4  •  10~^ft^fsec 

4  •  10~^ft^fs€c 

5,12,19,26 

4  •  10~^  ft^ / sec 

4  •  10~^  ft^lsec 

6,13,20,27 

0.0 

10-»/<V-sec 

7,14,21,28 

■m 

0.0 

10~'^  f  sec 

8,15,22,29 

^Vavi 

0.0 

10~^  f  sec 

9,16,23,30 

SZaui 

0.0 

I0~!^ft^fsec 

Table 

5.12.  DGPS  Stated  Tuning  Parameters  Q(t) 

State 

Numbe 

r 

, State 
Symbol 

Truth  Model 
Value 

Filter 

Value 

1 

0.0 

10~*/<^/sec 

2 

0.0 

I0‘^ft^/sec 

3i  8,13,] 

8 

SRtropi 

1  •  IQ-'^fi^jseC 

|||QjQ]jjy[QH22S||[| 

4,  9,14,] 

9 

^Riom 

4  •  10“‘‘/t^/sec 

15  •  \Q^*ft’^lsec 

iiaiPina 

SXavi 

0.0 

1  •  ft^  j  sec 

msm 

III 

wmm 

0.0 

1  •  10~'*/t^/sec 

7,12,17,; 

!2 

llliES9!illl 

0.0 

1  •  JO  !  sec 

The  last  table  in  this  section  lists  the  final  R(<,)  value  for  the  GPS/DGPS 
pseudorange  measurement  for  the  reduced-order  filters  implemented  in  this  thesis. 
The  full-order  GPS/DGPS  filters  measurement  noise  variance  is  always  9  The 
RRS  range  measurement  noise  variance  is  always  4  Only  the  69-state  NRS, 
48-state  ENRS,  ani  22-state  DGPS  reduced-order  filters  need  adjustment  of  their 
pseudorange  measurement  noise  variance,  so  only  values  for  these  filters  are  listed. 

Table  5.13.  GPS/DGPS  Measurement  Noise  Variance  Tuning  Parameters  R(<,) 


Filter 

Truth  Model 
Value 

Filter 

Value 

69-State  NRS 

9/i^ 

300/^2 

48-State  ENRS 

9/t" 

22-State  DGPS 

~loW~ 

5.7'  Summary 

This  chapter  discusses  the  results  of  all  the  filters  de.signed,  implemented,  and 
analyzed  in  this  thesis.  As  expected,  the  full-order  filters  perform  better  than  the 
reduced  order  filters,  while  the  filters  utilizing  DGPS  measurements  outperform  fil¬ 
ters  with,  only  GPS  measurements.  Again  as  a  reminder,  the  truth  model  for  th.e 
INS  is  41  states  (except  for  Stacey’s  72-state  model),  the, RRS  truth  model  is  26 
states,  the  GPS  truth  model  30  states,  and  the  DGPS  truth  model  is  22  states.  The 
anomaly  in  the  velocity-aiding  measurement  equations  hindered  the  analysis  of  tlie 
filters  implementing  this  type  of  measurement.  However,  every  filter  in  this  the^ 
sis  significantly  outperforms  CIRIS  and  show’s  the  benefit  of  augmenting  the  RRS 
transponder  measuieinents  with  GPS  and  DGPS  measurements. 


VI.  Conclusions  and  Peco  iimendations 

6. 1  Conclusions 

6.1.1  The  NR.^  Filters.  All  three  NRS  filters  performed  well  and  jneatly 
increzised  the  accuracy  of  the  navigation  Solution  over  CIRIS  as  it  is  implemented 
today.  The  following  sections  contain  specific  conclusions  drawn  from  the  NRS 
filters’  performance  while  recommendations  derived  from  their  implementation  arc 
found  later  in  this  chapter. 

6.1.I.1  ^1-State  INS  Model.  The  41-state  INS  error  model  proves  to 
be  an  adequate  substitute  for  the  93-state  LN-93  INS  error  model  in  both  the  truth 
and  filters  models  were  it  was  utilized.  As  shown  in  Appendix  C,  the  41-state  error 
model  followed  the  93-state  error  model  to  within  one  percent  over  the  entire  fighter 
flight  trajectory  with  respect  to  position,  velocity,  and  attitude  errors.  The  reader 
must  recognize  that  the  eliminated  states  from  the  93-state  model  feed  into  the  first 
ten  states  in  the  INS  error  model,  so  the  remaining  INS  states  (states  11  through  41 
in  Tables  A. 7  and  A.8  in  Appendix  A)  do  not  vary  from  the  93-state  error  model  and 
plots  from  these  states  are  not  included  in  Appendix  C.  This  reduction  in  states  in 
both  truth  and  filter  models  greatly  decrccised  the  time  required  for  filter  performance 
analysis.  . 

6.1.1.^  Single  versus  Double  Precision.  As  stated  in  Section  5.1.1,  dou¬ 
ble  precision  is  necessary  to  faithfully  inodel  GPS/DGPS  measurements  in  VAX  For¬ 
tran.  Recognizing  that  every  computer  system  implements  precision  differently,  the 
reader  is  warned  that  12-bit  accuracy  is  the  minimutn  accuracy  required  for  a  2-liour 
simulation  with  receiver  clock  errors  modeled  as  in  this  thesis.  However,  also  recog¬ 
nize  that  the  single  precision  assumption  used  in  the  full-order  filters  implemented 
in  this  thesis  only  affects  the  initial  filter  acquisition  (estimation)  of  the  receiver 
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clock  errors,  and  the  size  of  these  two  clock  errors  makes  the  decreased  magnitude 
assumption  a  viable  alternative  to  double  precision  simulation  analysis. 

6. 1.1. 3  NRS  Filter  Performance.  As  rtated  in  Chapter  V,  the  full-order 
128-  and  97-state  NRS  filters  perform  so  much  alike  that  any  discrepancy  in  perfor¬ 
mance  shown  in  the  tables  is  a  tuning  issue  and  noi  a  performance  issue.  However, 
the  69-state  NRS  filter  performance  degrades  significantly  when  compared  to  the  two 
full-order  filters  mentioned  above.  Even  though  the  69-state  filter  in  not  optimally 
tuned,  it  is  a  safe  assumption  that  its  position  and  velocity  error  estimation  accuracy 
will  be  worse  than  the  full-order  filters.  Recognize,  though,  that  the  69-state  filter  is 
at  least  3  times  more  accurate  than  the  current  CIRIS  system.  This  increase  in  per¬ 
formance  can  be  attributed  to  the  incorporation  of  GPS  pseudorange  measurerrtents 
to  augment  transponder  aiding  and  more  filter  states  used  for  the  analysis  (69  in  the 
NRS  versus  11  in  CIRIS).  How  much  each  of  these  attributes  individually  increases 
performance  is  not  in  the  scope  of  this  thesis;  however,  an  educated  guess  would  be 
that  increased  filter  states  give  the  most  increase  in  performance  in  evaluating  the 
69-state  NRS  filter’s  performance  increcise  over  CIRIS. 

6.1.2  The  ENRS  Filters.  Like  the  NRS  filters,  the  ENRS  filters  all  show 
a  great  increase  in  performance  over  CIRIS.  The  full-order  89-state  ENRS  filter 
outperforms  any  other  filter  implemented  in  this  thesis  by  producing  a  navigation 
solution  one  order  of  magnitude  more  accurate  than  CIRIS.  The  48-state  reduced- 
order  filter  IS  also  very  accurate,  and  produces  position  error  estimates  one  order 
of  magnitude  better  than  CIRIS  along  with  velocity  error  estimates  5  times  better 
(using  only  DGPS  pseudorange  and  RRS  range  measurements).  As  stated  in  Chapter 
V,  the  48-state  ENRS  filter  using  position-  and  velocity-aiding  measurements  has 
not  been  adequately  modeled  and  the  corruption  of  the  navigation  solution  by  an 
anomaly  eliminates  any  meaningful  conclusions  which  can  be  drawn  form  this  filter. 


However,  the  following  section  contains  a  conclusion  which  can  be  drawn  from  the 
ENRS  filters  successfully  implemented  in  this  thesis. 

6. 1.2.1  Utility  of  DGPS.  The  one  main  conclusion  which  can  be  drawn 
from  the  work  performed  in  this  thesis,  and  one  which  cannot  be  emphasized  enough, 
is  that  incorporation  of  DGPS  pseudorange  measurements  greatly  increases 
the  navigation  solution  accuracy  of  CIRIS.  When  either  the  89-sLate  full-order 
filter  or  the  48-state  reduced-order  filter  is  analyzed,  it  is  readily  apparent  that  DGPS 
provides  the  increase  in  navigation  solution  accuracy  required  for  CIRIS  to  be  used 
as  a  test  reference  against  the  GPS-aided  INSs  CIGTF  is  likely  to  be  testing  in  the 
next  decade.  This  one  conclusion  is  the  culmination  of  the  work  performed  in  thife 
thesis  and  the  3  other  theses  which  preceded  it  (20,  19,  22). 

6.2  Recommendations 

There  are  many  recommendations  which  could  be  derived  from  the  work  per¬ 
formed  in  this  thesis.  The  following  sections  contain  a  few  of  the  most  important  as 
determined  by  the  author  and  the  thesis  committee. 

6.2.1  Velocity-Aiding  Measurements  Require  Further  Study.  When  velocity- 
aiding  measurements  (DGPS  delta-range  and  RRS  range-rate)  were  implemented 
into  the  48-state  ENRS  filter,  corruption  of  both  position  and  velocity  error-state 
estimation  precluded  performance  analysis  of  the  filter.  At  first,  an  error  in  the 
algorithm  implementation  was  suspected;  however,  careful  analysis  proves  this  not  to 
be  the  case.  Then,  the  coordinate  transformation  from  the  navigation  frame  (NWU) 
to  Litton  ECEF  became  suspect.  When  a  different  transformation  was  implemented 
(results  not  shown  in  this  thesis),  the  corruption  of  position  and  velocity  states 
continued.  Finally,  during  the  thesis  defense,  the  correlation  between  position  and 
velocity  measurements  in  the  DGPS  subsystem  as  it  is  modeled  in  this  thesis  (and 
is  performed  by  the  ABR  in  real  life)  was  brought  to  the  surface.  Since  there  is  a 
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definite  correlation  between  position  and  velocity  measurements  in  the  truth  and  filter 
models  of  this  thesis,  it  is  necessary  to  model  them  as  being  correlated.  This  requires 
further  analysis  on  the  part  of  CIGTF  to  determine  to  true  correlation  values  so 
they  can  be  implemented  accurately  in  the  EKF.  Ihe  measurement  noise  covariance 
matrix  R(^)  is  modified  to  include  off-diagonal  terms  describing  the  correlation 
between  the  position  and  velocity  measurements.  Note  that  when  utilizing  MSOFE 
as  the  EKF  evaluation  tool,  scaler  measurements  are  assumed  (R(f,)  is  a  diagonal 
matrix)  so  that  for  correlated  measurements  (off-diagonal  terms),  this  matrix  raiist 
be  diagonalized  before  implementation.  Unfortunately,  Maybeck  has  stated  this  lack 
of  correlation  between  measurements  used  in  this  thesis  still  does  not  account  for  the 
corruption  of  the  navigation  solution,  since  the  truth  model  implements  uncorrelated 
measurements  as  well  (13).  This  anomaly,  which  has  plagued  the  three  previous 
theses  also  (20,  19,  22),  remains  undefined  and  a  subject  for  further  research. 

,  6.2.2  Return  to  a  93-State  INS  Truth  Model.  The  NRS  and  ENRS  truth  and 
filter  models  developed  or  analyzed  in  this  thesis  assume  that  either  a  72-state  INS 
model  (Stacey’s  NRS)  or  41-state  INS  model  (all  other  models  in  this  thesis)  perform 
close  enough  to  the  true  93-state  LN-93  error  model  so  that  any  real  difference  can 
be  ignored.  However,  before  real  meaisurements  are  used  to  test  any  of  the  filters, 
the  93-state  INS  model  should  be  included  in  the  truth  model  and  filter  performance 
against  a  full-order  truth  model  analyzed  in  case  further  tuning  is  required. 

6.2.3  Use  of  Real  Measurements.  Once  the  velocity-aiding  measurements’ 
anomaly  has  been  identified  and  corrected  and  the  filter  analyzed  against  the  true 
93-.state  INS  model,  actual  DGPS  and  RRS  mecisurements  should  be  used  to  test  the 
48-state  ENRS  filter  for  robustness.  It  is  a  foregone  conclusion  that  the  filter  will  need 
further  tuning  for  optimal  performance,  and  this  tuning  can  only  be  acconqdished 
using  real  measurements. 
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6.2.4  Reduction  of  the  RRS  error-states  in  the  Filter.  Another  worthwhile 
effort  could  be  the  reduction  of  the  26-.state  RRS  model  in  the  filter  to  a  2-state 
model.  The  RRS'  range  and  velocity  bias  states  would  be  retained  (see  Table  A. 5 
in  Appendix  A)  with  increased  variance  of  the  RRS  range  measurement  noise  to 
compensate  for  the  eliminated  states.  Recognize  this  reduction  does  not  affect  the 
RRS  range-rate  measurement  since  the  eliminated  states  don’t  appear  in  this  mea¬ 
surement.  Also,  since  DGPS  measurements  alone  provide  an  accurate  navigation 
solution,  only  a  minimal  degradation  of  the  navigation  solution  should  be  realized. 
Preliminary  work  using  this  technique  has  begun  utilizing  only  position-aiding  mea¬ 
surements  (not  included  in  this  thesis),  and  the  results  look  very  encouraging.  Note 
that  reduction  of  the  RRS  states  from  26  to  2  reduces  the  48-state  ENRS  filter  to 
24  states.  This  allows  for  inclusion  of  the  21  previously  eliminated  INS  states  (see 
Table  A.8  in  Appendix  A)  definitely  increasing  the  filter’s  ability  to  estimate  atti¬ 
tude  and  velocity  error-states  while  staying  under  the  70-state  post-processing  filter 
limit.  This  newly  proposed  filter  would  in  theory  have  41  INS  error-states,  2  RRS 
error-states,  and  2  DGPS  error-states  for  a  total  of  45  error-states  to  be  compared 
to  the  89-state  ENRS  truth  model  developed  in  Chapter  IV. 


Appendix  A.  Error  Model  State  Definitions 

This  appendix  contains  a  tabular  listing  of  the  93  Litton  INS,  26  RRS,  30  GPS, 
41  reduced  INS  truth  model,  and  the  22  DGPS  error-states. 

A.l  Litton  LN-93  Error-States 

Tables  A.l  through , A. 4  list  the  LN-93  error  model  (93  states)  as  defined  in 
the  Litton  CDRL  (11).  Note  that  this  document  cor  tains  several  errors,  which  have 
been  corrected  in  these  tables  (22). 

A. 2  RRS  Transponder  Error  States 

Table  A.5  lists  the  RRS  transponder  error  states  as  they  are  modeled  in  the 
NRS.  These  states  are  defined  in  and  extracted  from  (19).  A  total  of  26  states  are  in¬ 
cluded  to  model  the  error  characteristics  of  six  ground  transponders  plus  interrogator 
error  sources  (22). 

A. 3  GPS  Error  States 

Table  A. 6  lists  the  GPS  error  states  as  they  are  modeled  in  the  NRS.  These 
states  are  defined  in  and  extracted  from  (6).  The  definitions  are  believed  to  be 
ultimately  traceable  to  the  paper  by  D.B.  Cox  (4).  A  total  of  30  states  are  included  to 
model  the  error  characteristics  of  4  space  vehicles  plus  user  equipment  error  sources. 
(22) 

A. 4  Reduced  Order  INS  T'.'uth  Model  States 

Tables  A. 7  and  A. 8  list  the  41  INS  states  used  in  all  the  full-order  models  in 
this  thesis  except  for  Stacpy’s  128-state  model.  These  states  are  based  on  the  the 
recommendations  of  Lewantowicz  and  Keen  (10). 


A. 5  DGPS  Error  States 


Table  A. 9  lists  the  DGPS  error  states  as  they  are  modeled  in  the  ENRS.  A 
total  of  22  states  are  included  to  model  the  error  characteristics  of  4  space  vehicles 
plus  user  equipment  error  sources. 


State 

Sta 

Number 

Sym 

Table  A.l.  INS  System  Model:  INS  States  1 


Definition 


X-component  of  vector  angle  from  true  to  computer  frame 


Y-component  of  vector  angle  from  true  to  computer  frame 


Z-componcnt  of  vector  angle  from  true  to  computer  frame 


X-component  of  vector  angle  from  true  to  platform  frame 


Y-component  of  vector  angle  from  true  to  platform,  frame 


Z-compouent  of  vector  angle  from  true  to  platform  frame 


X-component  of  error  in  computed  velocity 


Y-component  cf  error  in  computed  velocity 


Z-component  of  error  in  computed  velocity  , 


Error  in  vehicle  altitude  above  reference  ellipsoid 


Error  in  lagged  inertial  altitude 


Error  in  vertical  channel  aiding  state 


Error  in  vertical  channel  aiding  state 


X-component  of  gyro  correlated  drift  rate 


Y-component  of  gyro  correlated  drift  rate 


Z-comnonent  of  gyro  correlated  drift  rate 


X-component  of  accelerometer  and 
velocity  quantizer  correlated  noise 


Y-  component  of  accelerometer  and 
velocity  quantizer  correlated  noise 


Z-component  of  accelerometer  and 
velocity  quantizer  correlated  noise 


X-component  of  gravity  vector  errors 


Y-component  of  gravity  vector  errors 


Z-component  of  gravity  vector  errors 


Total  baro- altimeter  correlated  error 


X-cdmponent  of  gyro  trend 


Y-component  of  gyro  trend 


Z-component  of  gyro  trend 


X-component  of  accelerometer  trend 


Y-component  of  accelerometer  trend 


Z-component  of  accelerometer  trend 


Table  A. 2.  INS  System  Model:  INS  States  30  — ♦  47 


State 

Number 

State 

Symbol 

Definition 

30 

6x 

X-component  of  gyro  drift  rate  repeatability 

31 

b. 

Y-component  of  gyro  drift  rate  repeatability 

32 

b. 

Z-component  of  gyro  drift  rate  repeatability 

33 

X-component  of  gyro  scale  factor  error 

34 

mmm 

Y-oomponent  of  gyro  scale  factor  error 

35 

Z-component  of  gyro  scale  factor  error 

36 

Xi 

X  gyro  misalignment  abput  Y-axis 

37 

X2 

Y  gyro  misalignment  about  X-axis 

38 

X3 

Z  gyro  misalignment  about  X-axis 

39 

‘'I 

X  gyro  misalignment  about  Z-axis 

40 

Vi 

Y  gyro  misalignment  about  Z-axis 

4,1 

Z  gyro  misalignment  about  Y-axis 

42 

Oxrx 

X  gyro  scale  factor  non-linearity 

43 

mam 

Y  gyro  scale  factor  ncfn-linearity 

44 

Z  gyro  scale  factor  non-linearity 

45 

mm 

X  gyre  scale  factor  asymmetry  error 

46 

mm 

Y'gyro  scale  factor  asymmetry  error 

47 

lamm 

Z  gyro  scale  factor  asymmetry  error 

Table  A. 3.  INS  System  Model:  NS  States  48  — >  69 


State 

Number 

State 

Symbol 

Definition 

48 

Vi. 

X-component  of  accelerometer  bias  repeatability 

49 

wEm 

Y-component  of  accelerometer  bias  repeatability 

50 

Vi. 

Z-component  of  accelerometer  bias  repeatability 

51 

X-component  of  accelerometer  and  velocity 
quantizer  scale  factor  error 

52 

5,. 

Y-component  of  accelerometer  and  velocity 
quantizer  scale  factor  error 

53 

5.4. 

Z-component  of  accelerometer  and  velocity 
quantizer  scale  factor  error 

54 

Sqa, 

X-component  of  accelerometer  and  velocity 
quantizer  scale  factor  asymmetry 

55 

Y-component  of  accelerometer  and  velocity 
quantizer  scale  factor  asymmetry 

56 

Sqa, 

Z-component  of  accelerometer  and  velocity 
quantizer  scale  factor  asymmetry 

57 

/xx 

Coefficient  of  er»‘or  proportional  to  square 
of  measured  acceleration 

58 

fyy 

Coefficient  of  error  oroportional  to  square 
of  measured  acceleration 

59 

/« 

Coefficient  of  error  proportional  to  square 
of  measured  acceleration 

60 

/xv 

Coefficient  of  error  proportional  to  products  of  acceleration 
along  and  orthogonal  to  accelerometer  sensitive  axis 

61 

H 

Coefficient  of  error  proporticbal  to  products  of  acceleration 
along  and  orthogonal  to  accelerometer  sensitive  axis 

62 

/vx 

Coefficient  of  error  proportional  to  products  of  acceleration 
along  and  orthogonal  to  acceleronicter  serisitive  axis 

63 

fy. 

Coefficient  of  error  proportional  to  products  of  acceleration 
along  and  orthogonal  to  accelerometer  sensitive  axis 

64 

U 

Coefficient  of  error  proportional  to  pioducts  of  acceleration 
along  and  oi  lhogonal  to  accelerometer  sensitive  axis 

m 

hy 

Coefficient  of  error  proportional  to  products  of  acceleration 
along  and  orthogonal  to  accelerometer  sensitive  axis 

66 

X  accelerometer  misalignment  about  Z-axis 

67 

Y  accelerometer  misalignment  about  Z-axis 

68 

Z  accelerometer  misalignnjc nt  about  Y-axis 

69 

(T^ 

Z-accclerometer  misalignment  about  X-axis 

Ttible  A.4.  INS  System  Model:  INS  States  70  -♦  93 


Definition 


State 

SI 

Number 

X-component  of  accelerometer  bias 
thermal  transient 


Y-component  of  accelerometer  bias 
therntal  transient 


Z-component  of  accelerometer  bias 
thermal  transient 


X-component  of  initial  gyro  drift  rate 
bias  thermal  transient 


Y-component  of  initial  gyro  drift  rate 
bias  thermal  transient 


Z-component  of  initial  gyro  drift  rat^ 
bias  thermal  transient 


X  gyro  compliance  term 


X  gyro  compliance  term 


X  gyro  compliance  term 


X  gyro  compliance  term 


X  gyro  compliance  term 


X  gyro  compliance  term 


Y  gyro  compliance  term 


Y  gyro  compliance  term 


Y  gyro  compliance  term 


Y  gyro  compliance  term 


Y  gyro  compliance  term 


Y  gyro  compliance  term 


Z  gyro  compliance  term 


Z  gyro  compliance  term 


Z  gyro  compliance  term 


Z  gyro  compliance  term 


Z  gyro  compliance  term 


Z  gyro  compliance  term 


93 


State  Sta 
Number  Sym 


Definition 

[NOTE:  Sims  =  Total  INS  States] 


Sins  +  1 

Range  error  due  to  equipment  bias 

Sins  +  2 

Sv\, 

Velocity  error  due  to  equipment  bias 

wsnm 

Transponder  1  x-component  of  position  error 

Sins  +  4 

Transponder  1  y-component  of  position  error 

Sins  +  5 

8Pt\, 

Transponder  1  z-romponent  of  position  error 

Sins  +  6 

Transponder  1  range  error  due  to  atm  propagation 

Sins  + 

Transponder  2  x-component  of  position  error 

Sins  +  8 

WSBM 

Transponder  2  y-component  of  position  error 

Sins  +  9 

8Pt2, 

Transponder  2  z-component  of  position  error 

Sins  +  10 

WESSSM 

msm 

MmSM 

Transponder  2  range  error  due  to  atm  propagation 

Sins  +  1 1 

Transponder  3  x-componeut  of  position  error 

Transponder  3  y-component  of  position  error 

Sins  +  13 

SPt3. 

Transponder  3  z-component  of  position  error 

SiNS  +  14 

102811 

13^91 

mssm 

Transponder  3  range  error  due  to  atm  propagation 

Sins  +  13 

Transponder  4  x-component  of  position  error 

Sins  +  16 

Transponder  4  y-component  of  position  error 

Sins  +  1 7 

8Pt4. 

Transponder  4  z-component  of  position  error 

Sins  -I- 18, 

IZ28B 

Transponder  4  range  error  due  to  atm  propagation 

Sins  +  19 

8Pts, 

Transponjler  5  x-component  of  ppsition  error 

Sins  +  20 

8Pts, 

Transponder  5  y-component  of  position  error 

Sins  +  21 

SPji. 

Transponder  5  z-component  of  position  error 

Sins  +  22 

msm 

Transponder  5  range  error  due  to  atm  propagation 

Sins  +  23 

8Pt^. 

Transponder  6  x-componeiit  of  position  error 

Sins  +  24 

1QS9I 

Transponder  6  y-component  of  position  erior 

Sms  +  25 

8Pt«. 

Transponder  6  z-component  of  position  error 

Sms  +  26 

I3S8I 

Transponder  6  range  error  due  to  atm  propagation 

Table  A.6.  GPS  Error  States 


State 

State 

Definition 

Number 

Symbol 

[NOTE:  Srrs  =  Total  RRS  States] 

\mmsmsm 

User  clock  bias 

EZSSHI 

User  clock  drift 

SV  1  code  loop  error 

Sins  +  Sbrs  +  4 

SV  1  tropospheric  error 

Sins  +  Srhs  +  5 

SV  1  ionospheric  error 

Sins  +  Srrs  +  6 

SV  1  clock  error 

Sins  +  Srrs  +  7 

SV  1  x-component  of  position  error 

Sins  +  Srrs  +  8 

SV  1  yrcomponent  of  position  error 

Sins  +  Srrs  +  9 

SV  1  zycomponent  of  position  error 

Sins  +  Srrs  +  10 

SV  2  code  loop  error 

Sins  +  Srrs  +  11 

SV  2  tropospheric  error 

Sins  +  Srrs  +  12 

SV  2  ionospheric  error 

Sins  +  Srrs  +  13 

SV  2  clock  error 

Sins  +  Srrs  +  14 

SV  2  x-component  of  position  error 

Sins  +  Srrs  +  15 

■m 

SV  2  y-component  of  position  error 

Sins  +  Srrs  +  16 

SV  2  z-component  of  position  error 

Sins  +  Srrs  +  17 

SV  3  code  loop  error 

Sins  +  Srrs  +  18 

^f^tropy 

SV  3  tropospheric  error 

Sins  +  Srrs  +  19 

mnm 

SV  3  ionospheric  error 

Sins  +  Srrs  +  20 

SV  3  clock  error 

.  Sins  +  Srrs  +  21 

SV  3  x-component  of  position  error 

Sins  +  Srrs  +  22 

^y»v3 

SV  3  y-comppnent  of  position  error 

Sins  +  Srrs  +  23 

^^SV3 

SV  3  z-component  of  position  error 

Sins  +  Srrs  +  24 

lOiil 

SV  4  code  loop  error 

Sins  +  Srrs  +  25 

SV  4  tropospheric  error 

Sins  +  Srrs  +  26 

mmm 

SV  4  ionospheric  error 

Sins  +  Srrs  +  27 

iZim 

SV  4  clock  error 

Sins  +  Srrs  +  28 

warn 

SV  4  x-component  of  position  error 

Sins  +  Srrs  +  29 

KSfflli 

S V  4  y-component  of  position  error 

Sins  +  Srrs  +  30 

SV  4  z-component  of  position  error 

Table  A. 7.  Reduced-Order  INS  System  Model:  INS  States  1  —*  20 


State 

Number 

Definition 

1 

S0^ 

X-component  of  vector  angle  from  true  to  computer  frame 

2 

K&iH 

Y-component  of  vector  angle  from  true  to  computer  frame 

3 

60, 

Z-component  of  vector  angle  from  true  to  computer  frame 

4 

<t>x 

X-component  of  vector  angle  from  true  to  platform  frame 

5 

4>y 

Y-component  of  vector  angle  from  true  to  platform  frame 

6 

4>Z 

Z-component  of  vector  angle  from  true  to  platform  frame 

7 

6V,, 

X-component  of  error  in  computed  velocity 

8 

mm 

Y-component  of  error  in  computed  velocity 

9 

6V, 

Z-compoiient  of  error  in  computed  velocity 

10 

6h 

Error  in  vehicle  altitude  above  reference  ellipsoid 

11 

Shi 

Error  in  lagged  inertial  altitude 

12 

6Sz 

Error  in  vertical  channel  aiding  state 

13 

6S^ 

Error  in  vertical  channel  aiding  state 

14 

ML 

X-ccmponent  of  2u:celerometer  and 
velocity  quantizer  correlated  noise 

15, 

^9 

Y-component  of  accelerometer  and 
velocity  quantizer  correlated  noise 

16 

V,. 

Z-component  of  accelerometer  and 
velocity  quantizer  correlated  noise 

17 

6gx 

X-component  of  gravity  vector  errors 

18 

Sgy 

Y-component  of  gravity  vector  errors 

19 

6g, 

Z-component  of  gravity  vector  errors 

20 

6hi 

Total  baro-altimeter  correlated  error 

A-9 


Table  A.8.  Reduced-Order  INS  System  Model:  INS  States  21  — +  41 


State 

Number 

jjjygs 

Definition 

21 

X-component  of  gyro  drift  rate  repeatability 

22 

■sH 

Y-component  of  gyro  drift  rate  repeatability 

23 

Z-component  of  gyro  drift  rate  repeatability 

24 

s,. 

X-component  of  gyro  scale  factor  error 

25  , 

Y-component  of  gyro  scale  factor  error 

26 

Z-component  of  gyro  scale  factor  error 

27 

V6, 

X-component  of  accelerometer  bias  repeatability 

28 

■an 

Y-component  of  accelerometer  bias  repeatability 

29 

Z-component  of  2u:celerometer  bias  repeatability 

30 

X-component  of  accelerometer  and  velocity 
quantizer  scale  factor  error 

31 

5^. 

Y-component  of  accelerometer  and  velocity 
quantizer  scale  factor  error 

32 

m 

Z-component  of  accelerometer  and  velocity 
quantizer  scale  factor  error 

33 

X-component  of  accelerometer  and  velocity 
quantizer  scale  factor  asymmetry 

34 

m 

Y-component  of  accelerometer  and  velocity 
quantizer  scale  factor  asymmetry 

35 

m 

Z-component  of  accelerometer  and  velocity 
quantizer  scale  factor  aisymmetry 

36 

A*! 

X  au:celerdnieter  misalignment  about  Z-axis 

37 

Y  awrcelerometer  misalignm'^nt  about  Z-axis 

38 

A<3 

Z  accelerometer  misalignment  about  Y-axis 

39 

<^1 

X-accelerometer  misalignment  about  Y-axis 

40 

<72 

Y-accelerometer  misalignment  about  Xraxis 

41 

<^3  . 

Z- accelerometer  misalignment  about  X-axis 

Table  A. 9.  DGPS  Error  States 


State 

Definition 

'  Number 

BHil 

[NOTE:  Srrs  =  Total  RRS  States] 

ABR  clock  bias 

WHBIKHM 

ABR  clock  drift 

Sms  +  Sfifis  +  3 

SV  1  tropospheric  error 

Sins  +  Srrs  +  4 

^  ^ion\ 

SV  1  ionospheric  error 

Sins  +  Srrs  +  5 

^Xjy, 

SV  1  x-component  of  position  error 

Sins  +  Srrs  -h  6 

liiS&IH 

SV  1  y-component  of  position  error 

Sins  +  Srrs  +  7 

SV  1  zycomponent  of  position  error 

Sins  +  Srrs  +  8 

SV  2  tropospheric  error 

Sins  +  Srrs  +  9 

SV  2  ionospheric  error 

Sins  +  Srrs  +  10 

SV  2  x-component  of  position  error 

Sins  +  Si^^is  +  11 

^y$vj 

SV  2  )'-component  of  position  error 

Sins  +  Srrs  +  12 

SV  2  z-compcnent  of  position  error 

Sins  +  +  13 

^fitrop3 

SV  3  tropospheric  error 

Sins  +  •S’/jas  +  14 

SV  3  ionospheric  error 

Sins  I-  Srrs  +  15 

,SV  3  x-component  of  position  error 

Sins  +  +  16 

^y>v^ 

SV  3  y-component  of  position  error 

Sins  +  Srrs  +  17 

SV  3  z-component  of  position  error 

Sins  f  •S'rr5  +  18 

^  f^trop\ 

SV  4  tropospheric  error 

Sins  +  5/ir5  +  19 

SV  4  ionospheric  error 

Sins  +  5rr5  +  20 

SV  4  x-component  of  position  error 

Sins  +  ■S’rr5  +  21 

^y$vi 

SV  4  y-component  of  position  error 

Sins  +  Srrs  +  32 

^ESHI 

SV  4  z-component  of  position  error 

All 


Appendix  B.  Litton  LN-93  Error-State  Model  Dynamics  Matrix 

(22) 

The  LN-93  error-state  dynamics  matrix  (F)  as  provided  by  Litton  is  a  93-by- 
93  Mray  that  contciins  a  large  number  of  elements  that  are  identically  zero.  Litton 
partitions  the  F  matrix  into  thirty-six  subarrays  (11)  reflecting  the  logical  divisions 
of  error  sources  discussed  in  Chapter  III. 

The  reader  should  note  that  only  the  NON-ZERO  elements  are  included  in  the 
tables  w^ich  follow,  and  should  further  note  that  the  revised  baro-altimeter  model 
states  are  NOT  included  in  this  set  of  ORIGINAL  F  matrix  elements  extracted  from 
the  Litton  document  (11).  ■ 

A  notational  convention  (22)  is  to  label  elements  of  the  C*  ,  sensor-to-true, 
matrix  as  where  i  is  the  row  and  j  is  the  column  in  the  transformation  matrix. 


Table  B.l.  Elements  of  the  Dynamics  Submatrix  Fn  (22) 


Element 

Term 

Element 

Term 

(1,3) 

-Py 

(1,8) 

-Cry 

(2,3) 

PT 

(2,7) 

Crx 

(3,1) 

py 

(3,2) 

-Px 

(4,2) 

—ill 

(4,3) 

ily 

(4,5) 

(4,6) 

(4,8) 

-Cry 

(5,1) 

n. 

(5,3) 

-il. 

(5,4) 

-Uit. 

(5,6) 

(5,7) 

Crx 

(6,1) 

-ily 

(6,2) 

fix 

(6,4) 

(6,5) 

-w.t. 

(7,1) 

-2Kiiy  ~  2V,ni 

(7.2) 

2v;n, 

^7,3) 

2Viay 

(7,5) 

-Ai 

(7,6) 

A, 

(7,7) 

-KCrx 

(7,8) 

2n, 

(7,9) 

—py  —  2ily 

(8,1) 

2v;fi, 

(8,2) 

-2Vxilr  -  2ViQ, 

(8,3) 

2Viily 

(8,4) 

Ai 

(8,6) 

-Ai 

(8,7) 

-2D, 

(8,8) 

—yiCRY 

(8,9) 

Px  +  2Dr 

(9,1) 

2VM, 

(9,2) 

2v;n, 

(9,3) 

-2Vyily  -  2l4n, 

(9,4) 

-Ay 

(9,5) 

-4r 

(9,7) 

Py  +  2Sly  +  VxCrX 

(9,8) 

~Pz  —  2Ux  +  VyCRY 

(9,10) 

2go/a 

(9,11) 

—k^ 

(9,12) 

-1 

(9,13) 

k2 

(10,9) 

1 

(10,11) 

-ki 

(10,13) 

(11,10) 

1 

(11,11) 

-1 

(12,11) 

k. 

(12,13) 

-^3 

(13,10) 

k4 

(13,11) 

(13,13) 

k4-l 

!■■■ 

■■■■■I 

Table  B.2.  Elements  of  the  Dynamics  Submatrix  F12  (22) 


Element  Term  Element  Term  Element  Term 


(4,14) 


(4,24) 


(5,14) 


(5,24) 


(6,14) 


(6,24) 


(7,17) 


(7,20) 


(7,29) 


(8,19) 


(8,28) 


(9,18) 


(9,23) 


(9,29) 


(13,23) 


C33< 


*4/600 


(4,15) 


(4,25) 


(5,15) 


(5,25) 


(6,15) 


(6,25) 


(7,18) 


(7,27) 


(8,17) 


(8,21) 


(8,29) 


(9,19) 


(9,27) 


(10,23) 


(4,16) 


(4,26) 


(5,16) 


(5,26) 


(6,16) 


(6,26) 


(7,19) 


(7,28) 


(8,18) 


(8,27) 


(9,17) 


(9,22) 


(9,28) 


(12,23) 


Table  B.3.  Elements  of  the  Dynamics  Submatrix  Fja 


Table  B.4.  Elements  of  the  Dynamics  Submatrix  Fh  (22) 


Element 


(7,48) 


(7,51) 


(7,54) 


(7,57) 


(7,60) 


(7,63) 


(7,66) 


(7,69) 


(8,50) 


(8,53) 


(8,56) 


(8,59) 


(8,62) 


(8,65) 


(8,88) 


(9,49) 


(9,52) 


(9.55) 


(9,58) 


(9,61) 


(9,64) 


(9,67) 


Term 


Element 


(7,49) 


(7,52) 


(7,55) 


(7,58) 


(7,61) 


Term 


Ci2 


c 


C^^A^A 


1 


n 

(7,67) 

—CnA 

n 

(7,68) 

C12A 

(8,48) 


(8,51) 


(8,54) 


(8,60) 


(8,63) 


(8,66) 


(8,69) 


(9,50) 


(9,53) 


(9,56) 


(9,59) 


(9,82) 


(9,65) 


(9,68) 


Cn\A^\ 


C^tAfAf 


CxsA 


Element 


(7,50) 


(7,53) 


(7,56) 


(7,59) 


(7,62) 


(7,65) 


Term 


C'ls 


C 


CxMT 


c 


(8,49) 


(8,52) 


(8,55) 


(8,58) 


(8,61) 


(8,64) 


(8,67) 


(9,48) 


(9,51) 


(9,54) 


(9,63; 


(9,66) 


(9,69) 


Table  B.5.  Elements  of  the  Dynamics  Submatrix  Fjs  (22) 


Table  B.6.  Elements  of  the  Dynamics  Submatrix  Fie  (22) 


Element 


Term  Element  Term  Element  Term 


(4,76) 


(4,79) 


(4,82) 


(4,85) 


(4,88) 


(5,76) 


(5,79) 


(5,82) 


(5,85) 


(5.88) 


(5,91) 


(6,76) 


(6,79) 


(6,82) 


(6,85) 


(6,88) 


(6,91) 


xmsmk\ 

iBspi 


IICaWjMli 

msmi\ 

\Msasmm 


(4,83) 


(4,86) 


(4,89) 


(4,92) 


(5,77) 


(5,80) 


(5,83) 


(5,86) 


(5,89) 


(5,92) 


(6,77) 


(6,80) 


(6,83) 


(6,86) 


(6,89) 


(6,92) 


C\iA^  Uib. 


(^22  A  f  Ufa  ^ 


Cai/lfwii, 


C32/lfu;,i,, 


(4,78) 


(4,81) 


(4,84) 


(4,87) 


(4,90) 


(4,93) 


(5,78) 


(5,81) 


(5,84) 


(5,87) 


(5,90) 


(5,93) 


(6,78) 


(6,81) 


(6,84) 


(6,87) 


(6,90) 


(6,93) 


Table  B.7,  Elements  of  the  Dynamics  Submatrix  F22  (22) 


Element  Term  Element  Term  Element  I  Term 


(H,14) 


(17,17) 


(20,20) 


(23,23) 


(18,18)  ■  -/? 


(21,21)  1 


(16,16)  -/?, 


(19,19) 


(22,22) 


Table  B.8.  Ehments  of  the  Dynamics  Submatrix  F55  (22) 


B.l  Elements  of  the  Process  Noise  Matrix 

The  Litton  document  (11)  includes  a  93-by-93  process  noise  matrix  (Q)  for  the 
LN-93  error  model.  Like  the  F  matrix,  the  Q  matrix  is  partitioned  into  subarrays 
which  correspond  to  the  error-state  subvectors  discussed  in  Chapter  HI.  The  vast 
majority  of  the  elements  in  Q  are  identically  zero.  Only  the  non-zero  elements  of  Q 
aie  shown  below. 


Table  B.9.  Non-zero  Elements  of  Process  Noise  Submatrix  Qn  (22) 


Element 

Term 

Kieliullt 

Term 

(4,4) 

(o,5) 

(6,6) 

(7,7) 

(8,8) 

(9,9) 

< 

Table  B.IO.  Non-zero  Elements  of  Process  Noise  Submatrix  (22) 


Element 

Term 

Element 

'I'erm 

(14,14) 

r('l5,15) 

(16,16) 

(17,17) 

(18,18) 

(1H,18) 

(20,20) 

(21,21) 

(22,22) 

(23,23) 

Appendix.  C.  Compariscn  of  93-state  and  41-State  INS  Models 

This  appendix  contains  plots  of  position,  velocity,  and  attitude  errors  of  the 
93-state  and  41-state  Litton  Ln-93  INS  truth  models  during  a  10-run  Monte  Carlo 
2-hour  flight  profile  (11).  No  alignment  was  performed  previous  to  the  start  of  the 
flight,  causing  the  large  magnitudes  in  all  the  errors.  These  plots  show  that  the  41- 
state  truth  model  almost  exactly  matches  the  93-state  truth  model,  and  it  is  suitable 
for  use  as  the  truth  model  for  all  the  filters  developed  in  this  thesis. 


Alttlud*  (ft)  Lx>ncltude  (decree*)  Latitude  (decree*) 


Figure  C.l,  2-Hour  Fighter  Flight  Profile  (a)  Latitude  (b)  Longitude  and  (c) 
Altitude 


C-2 


UtUlu*  bnr  (It) 


Appendix  D.  Stacey’s  12S- State  NRS  Filter  Performance  Plots 

All  plots  contained  in  this  and  subsequent  appendices  are  discussed  in  Chapter 
V  (Results).  All  plots  contained  in  this  and  subsequent  appendices  contain  five  traces. 
The  innermost  trace  (-  -  -)  on  each  data  plot  is  the  mean  error  time  history  for  the 
applicable  state.  Mean  Error  is  defined  as  being  the  difference  between  the  filters 
estimate  of  the  state  and  the  true  state  averaged  over  the  number  of  Monte  Carlo 
runs  performed.  The  equation  describing  this  relationship  is  defined  by  (14,  22): 

=  =  (D.l) 

where  is  the  filter-computed  estimate  of  variable  j  and  X(ruej(^)  is  the  truth 

model  value  of  the  same  variable,  at  time  t,,  for  sample  ;,  and  N  is  the  number  of 
time  histories  in  the  simulation  (10  in  this  thesis). 

In  addition  to  the  center  trace,  two  more  pairs  of  traces  are  plotted.  The  first 
pair  (represented  by  •  •  •)  is  symmetrically  displaced  about  the  mean  and  as  a  result 
follows  the  “undulations”  of  the  A/e(t,).  The  locus  of  these  traces  is  calculated  from 
(15,  22)  Me{ti)  ±  iJPJti)  ,  where  /’,(<<)  is  the  true  error  covariance  at  time  t,.  The 
true  standard  deviation  is  calculated  from  (14,  22): 

^ 

where  N  is  the  number  of  runs  in  the  Monte  Carlo  simulation  (10  in  this  thesis), 
and  Mj{ti)  is  the  mean-squared  value  of  the  Variable  at  each  time  of  interest  (such 
as  measurement  times). 

The  last  pair  of  traces  ( — )  represents  the  filter  computed  ±  (T/uter  values  for 
the  same  variables  of  interest  and  are  symmetrically  displaced  about  zero  because 
fhe  filter  “believes”  that  it  is  producing  zero-mean  errors  (16,  22).  These  quantities 


are  propagated  and  updated  in  the  MSOFE  (3,  22)  software  using  the  covariance 
propagation  equation  shown  in  Chapter  II.  These  traces  represent  the  filter’s  estimate 
of  its  own  error. 

D.l  Stacey’s  128-State  NRS  Filter  Performance  Plots 

The  plots  in  this  section  represent  results  of  a  10- run  Monte  Carlo  2- hour  fighter 
flight  profile  simulation.  In  these  runs  the  GPS  receiver  clock  states  initial  covariance 
values  have  been  decreased  4  orders  of  magnitude,  as  discussed  in  Chapter  V,  so  that 
single  precision  could  be  used  in  the  MSOFE  simulation.  This  filter  incorporater  both 
RRS  range  and  GPS  pseudorange  measurements. 


Figufe  D.3.  128-State  NRS  Filter  (a)  East  Velocity  (b)  North  Velocity  and  (c) 
Vertical  Velocity  Errors. 
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Figure  D.5.  128-State  NRS  Filter  (a)  Lagged  Altitude  (b)  Del  S3  and  (c)  Baro- 
Altimeter  Errors. 


Figure  D.6,  128-State  NRS  Filter  (a)  RRS  Range  Bias  (b)  RRS  Vel  Bias  and  (c) 
Tl  Atmosphere  Errors. 
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Figure  D.7.  128.State  NRS  Filter  (a)  Tl  X  Position  (b)  Tl  Y  Position  and  (c)  Tl 
Z  Position  Errors. 


! 


' 

1 : 

f 

: 

'  V  ..  . 

i 

» *  '  •  *  ■ » •  •  * 

“•  ( 

Figure  D.9. 

>  iboo 

128-State  NR 
SVl  Clock  E 

»boo  30<»0 

S  Filter  (a)  S’ 
rors. 

Bfl 

Mean  Er 

n 

1 

Mean  E. 

— 

0  i  ^ JxUrr 

i: 

Troposphere  (b)  SVl  Ionosphere  and  (c) 


ror  =  Mr  -  (Mr)trur 
^or  ± 


-11 


Appendix  E.  97-State  NRS  Filter  Performance  Plots 


The  plots  in  this  section  represent  results  of  a  10-run  Monte  Carlo  2-hour  fighter 
flight  profile  simulation.  In  these  runs  the  GPS  receiver  cloc'k  states  initial  covariance 
values  have  been  decreased  4  orders  of  magnitude,  as  discussed  in  Chapter  V,  so  that 
single  precision  could  be  used  in  the  MSOFE  simulation.  This  filter  incorporates  both 
RRS  range  and  GPS  pseudorange  measureiiients. 
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Figure  E.l.  2-Hour  Fighter  Flight  Profile  (a)  Latitude  (b)  Longitude  and  (c) 
Altitude 


97-State  NRS  Filter  (a)  Latitude  (b)  Longitude  and  (c)  Altitude  Errors. 
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Figure  E.7.  97-State  NRS  Filter  (a)  Tl  X  Position  (b)  Tl  Y  Position  and  (c)  Tl 
Z  Pof■'^Jn  Errors. 


Figure  E.8.  97-State  NRS  Filter  (a)  User  Clock  Bias  (b)  User  Clock  Drift  and  (c) 
SVl  Code  Loop  Errors. 
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Appendix  F.  69-Stnte  NRS  Filter  Performance  Plots 

The  plots  in  this  section  represent  results  of  a  10-run  Monte  Carlo  2-hour  fighter 
flight  profile  simulation.  This  filter  analysis  was  performed  in  double  precision,  so 
the  GPS  receiver  clock  errors  are  at  their  actual  values.  This  reduced-order  filter 
incorporates  both  RRS  range  and  GPS  pseudorange  measurements. 
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Figure  F-5.  69-State  NRS  Filter  (a)  Lagged  Altitude  (b)  Del  S3  and  (c)  Baro- 
Altimeter  Errors. 
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Appendix  G.  89-State  ENRS  Filter  Performance  Plots 


The  plots  in  this  section  represent  results  of  a  10-run  Monte  Carlo  2-hour  fighter 
flight  profile  simulation.  In  these  runs  the  GPS  receiver  clock  states  initial  covariance 
values  have  been  decreased  4  orders  of  magnitude,  as  discussed  in  Chapter  V,  so  that 
single  precision  could  be  used  in  the  MSOFE  simulation.  This  filter  incorporates  both 
RRS  rauige  and  DGPS  pseudorange  measurements. 
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Figure  G.l.  2-Hour  Fighter  Flight  Profile  (a)  Latitude  (b)  Longitude  and  (c) 
Altitude 


G-2 


Figure  G.2.  89-State  ENR5  Filter  (a)  Latitude  (b)  Longitude  and  (c)  All 
Errors. 
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Figure  G.3.  89-State  ENRS  Filter  (a)  East  Velocity  (b)  North  Velocity  and  (c) 
Vertical  Velocity  Errors. 
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Figure  G.9.  89-State  ENRS  Filter  (a)  SVl  Troposphere  and  (b)  SVl  Ionosphere 
Errors. 


Appendix  H.  4^-State  ENRS  Filter  Performance  Plots 

The  plots  in  this  section  represent  results  of  a  10-run  Monte  Carlo  2-hour  fighter 
flight  profile  simulation.  This  filter  analysis  was  performed  in  double  precision,  so 
the  GPS  receiver  clock  errors  are  at  their  actual  values.  This  reduced-order  filter 
incorporates  both  RRS  range  and  DGPS  pseudorange  measurements. 
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Figure  H.l. 


Figure  H.2.  48-State  ENRS  Filter  (a)  Latitude  (b)  Longitude  and  (c)  Altitude 
Errors. 
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Figure  H.5.  48-State  ENRS  Filter  (a)  Lagged  Altitude  (b)  Del  S3  and  (c)  Baro 
AUimeter  Errors, 
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Figure  H.8.  48-State  ENRS  Filter  (a)  Usei  Clock  Bias  and  (b)  User  Clock  Drift 
Errors, 
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Appendix  1.  22-State  DGPS  Filter  Performance  Plots 


The  plots  in  this  appendix  present  results  of  10-run  Monte  Carlo  simulations 
for  a  22-state  reduced-order  filter  used  to  implement  delta-range  measurements.  The 
first  section  contains  plots  from  the  2-hour  fighter  flight  profile.  The  second  section 
contains  plots  from  a  2-hour  straight  trajectory.  Both  filter  analyses  were  performed 
in  double  precision,  so  the  DGPS  receiver  clock  errors  are  at  their  actual  values. 
These  reduced-order  filters  incorporate  both  DGPS  pseudorange  and  delta-range 


measurements. 
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Figure  1. 1.  2-Hour  Fighter  Flight  Profile  (a)  Latitude  (b)  Longitude  and  (c) 
.Altitude 
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Figure  1.2.  22-State  DCIPS  Filter,  Fighter  Might  Profile  (a)  Latitiule  (h)  Longitude 
and  (c)  Altitude  F'rrors. 
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Figure  1.3.  22-State  DCPS  Filter,  Fighter  Flight  Profile  (a)  East  Velocity  (b)  North 
Velocity  and  (c)  Vertical  Velocity  Errors. 
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Figure  1.6.  22-State  DGPS  Filter,  Fighter  Flight  Profile  (a)  User  Clock  Bias  and 
(b)  User  Clock  Drift  Errors. 
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Figure  1.7.  2-Hour  Straight  Flight  Profile  (a)  Latitude  (b)  Longitude  and  (c), 
Altitude 
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Figure  1.8.  22-State  DGPS  Filter,  Straight  Flight  Profile  (a)  Latitude  (b)  Longitude 
and  (c)  Altitude  Errors. 
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Figure  1.9.  22-State  DGPS  Filter,  Straight  Flight  Profile  (a)  Ecist  Velocity  (b) 
North  Velocity  and  (c)  Vertical  Velocity  Errors. 
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Appendix  J.  4^-State  CIRIS  Filter  Performance  Plots 

The  plots  in  this  appendix  present  results  of  10-run  Monte  Carlo  simulations 
for  a  46-state  full-order  filter  used  to  implement  range-range  measurements.  The 
first  section  contains  plots  from  the  2-hour  fighter  flight  profile.  The  second  section 
contains  plots  from  a  2-hour  racetrack  trajectory.  Both  filter  analyses  were  performed 
in  single  precision  since  no  DGPS  states  are  modeled.  These  reduced-order  filters 
incorporate  both  RRS  range  and  range-rate  measurements. 
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Figure  J.4.  46-Stale  CHUS  Filter,  Fighter  F'light  Profile  Filter  (a)  F’ast  Tilt  (b) 
North  Tilt  and  (c)  Azimuth  F>roni. 
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Figure  J.5.  -le-State  CIRIS  Filter,  Fighter  Flight  Profile  Filter  (a)  flagged  Altitude 
(b)  Del  S3  and  (c)  Baro- Altimeter  Errors. 
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Figure  J.8. 
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Figure  J.ll.  46-State  CIRIS  Filter,  Racetrack  Flight  Profile  Filter  (a)  East  Tilt  (b) 
North  Tilt  and  (c)  Azimuth  Errors. 
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Figure  J.12.  46-State  CIRIS  Filter,  Racetrack  Flight  Profile  Filter  (a)  Lagged  Al¬ 
titude  (b)  Del  S3  and  (c)  Baroy Altimeter  Errors. 
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Appendix  K.  ENRS  Filter  Performance  Plots 


.  The  plots  in  this  section  represent  results  of  a  25-run  Monte  Carlo  2-hour 
fighter  flight  profile  simulation.  This  filter  analysis  was  performed  in  double  preci¬ 
sion,  so  the  DGPS  receiver  clock  errors  are  at  their  actual  values.  This  reduced-order 
filter  incorporates  RRS  range/range-rate  and  DGPS  pseudorange/delta- range  mea¬ 
surements. 
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Figure  K.l.  2-Hour  Fighter  Flight  Profile  (»)  Latitude  (b)  Ix>ngitude  and  (r), 
Altitude 


Figure  K.2.  48-SUte  ENRS  Filter  (a)  Latitude  (b)  l/ongitude  and  (c)  Altitude 
Errors. 
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Figure  K.5.  48-State  ENRS  Filter  (a)  Lagged  Altitude  (b)  Del  S3  and  (c)  Baro- 
Altimeter  Errors, 


Figure  K.6.  48-State  ENPi!  Filter  (a)  RRS  Range  Bias  (b)  RRS  Vel  Bias  and  (c) 
Tl  Atmosphere  Errors. 


Figure  K.7.  48-State  ENRS  Filter  (a)  Tl  X  Position  (b)  T1  Y  Position  and  (c)  Tl 
Z  Position  Errors. 
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